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Space Vehicle Power Systems 


EUGENE B. ZWICK’ and 
ROBERT L. ZIMMERMAN? 


Sundstrand Turbo Division 
Pacoima, Calif. 


Mr. Zwick is in charge of the Research Branch at Sundstrand Turbo. He received his B.S. in 
Physics from California Institute of Technology. Before joining Turbo in 1954 he was with 
Marquardt Aircraft Co. where he was engaged in the fields of aerodynamics and combustion. 


Mr. Zimmerman is in charge of the Systems Integration Branch at Sundstrand Turbo which 


handles all development projects. 


Mr. Zimmerman received his B.S. degree from University of 


California at Berkeley and before joining Turbo in 1953 had spent a number of years with Hughes 
Aircraft Co. as a specialist in electronic systems. 


pein es herein is a review of systems which can be 
utilized for the generation of secondary and/or propulsion 
power for advanced space vehicles. The power systems 
considered are: 

Chemically-fueled power supplies; open and closed cycle. 

Solar power supplies, closed cycle. 

Nuclear power supplies, closed cycle. 

Thermionic and thermoelectric power supplies. 

Solar photo-voltaic and solar recycling fuel cell power 

supplies. 

Included also in the review are the results of an initial pro- 
gram relating to the optimum temperature ratio, collector 
and radiator size for general power systems when used with 
nuclear, chemical or solar energy sources. The performance 
of the conversion devices per se as a function of the collector 
or radiator environment is not included. 

Descriptions are presented of a number of typical power sys- 
tems. An attempt is made to describe the general physical 
construction as well as to indicate the limitations of the sys- 
tem. Analyses and estimates have been prepared in order 
to determine the system weight as a function of power level 
and duration. 

Comparison of some of the systems studied is presented. 
Weights and sizes of the significant components are tabulated 
for |- and 100-kw systems. Variation of system weight with 
time and power level is plotted to determine the optimum 
system from a weight standpoint for various applications. 

System selection for various missions is summarized, and 
the role of turbomachinery in space vehicle power systems is 
noted. 


Description of Systems 


Many devices which are usable for conversion of energy to 
electrical power have been given serious consideration by 
various research groups and industry within this country and 
Europe. These devices, which are reviewed herein, can be 
categorized either as: Thermal engines including turbines, 
thermionic and thermoelectric generators; electrochemical 
engines including fuel cells and solar regenerated fuel cells; 
and quantum engines including photo-voltaic cells. 

The energy sources for these devices can be derived from a 
nuclear reactor, the combustion or decomposition of chemi- 


, Presented at the ARS Semi-Annual Meeting, June 8-11, 1959, 
San Diego, Calif. 
' Assistant Chief Engineer, Research Branch. Member ARS. 
Assistant Chief Engineer, Systems Integration Branch. 
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cals, the thermal energy from a solar source, or from a discreet 
quantum of energy at a particular frequency. 

Future space vehicles, depending on the application, can 
utilize the conversion devices or energy sources that have 
been mentioned. The determination of the most suitable de- 
vice for the advanced space vehicle depends on factors, such as 
watts per pound. duration, reliability of the device, effect of a 
malfunction on the society living within the regime of the de- 
vice (such as a malfunction of a nuclear reactor at launch 
when used as a space vehicle power supply), availability 
of the system and, of course, cost, environment and other 
usual considerations. 

The power requirements for space vehicles can be divided 
into two categories, the secondary and propulsion require- 
ments. The secondary power requirements are shown in 
Table 1. These range from 1 to 200 kw for periods of 4 hr to 1 
year. 

For electrical propulsion systems, however, depending on 
the length of the mission and the boost thrust available, the 
power will vary from 50 to 600 kw. The power as a function 
of specific impulse is shown in Fig. 1. The mass flow rate is 
inversely proportional to the specific impulse. The power, 
and therefore the prime mover weight, however, is propor- 
tional to the specific impulse. Therefore, the prime mover 
weight will increase as the total requirement for working fluid 
increases (extended missions) and the weight rate decreases. 
An optimum power level exists which depends on the duration 
of the mission. 

The important fact is that power ranging from one to hun- 
dreds of kilowatts for periods from a few hours to years will be 
required. 

Systems considered are: 

An open cycle chemically-fueled power supply wherein 
oxygen and hydrogen are used as an energy source. The 
combustion products are used to drive a rotary prime mover 
and thence to generate electrical power. 

A closed cycle chemically-fueled power supply wherein 
oxygen and hydrogen are used as an energy source to vapor- 
ize a working fluid (liquid sulphur). The recirculated work- 
ing fluid is used to drive a rotary prime mover and thence to 
generate electrical power. 

A closed cycle solar power supply wherein energy from the 
sun is collected and concentrated by a parabolic mirror. The 
concentrated energy is supplied to a boiler to heat the working 
fluid and then to generate power using the same mechanism 
as in the closed cycle chemically-fueled power supply. 
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Mission 


manned satellite 

cislunar vehicle (manned) 
lunar orbit probe (unmanned ) 
circumlunar vehicle (manned ) 


interplanetary probe to Mars 
interplanetary vehicle (manned) 


Table 1 Power requirements for space vehicles 


soft lunar landing and return (manned ) 


Time kw 
4-30 hr 1-15 
150 hr 3-30 
300 hr 1-2 
300 hr 3-30 
300 hr 10-50 
300 days 2-3 
300 days 10—200 


A closed cycle nuclear power supply wherein the energy 
from a nuclear reactor is transferred to the working fluid. 
The energy in the working fluid is used in the same manner 
as it is in the closed cycle chemically-fueled and solar power 
supply systems. 

A solar cell power supply wherein the incident photons from 
the sun are used to displace electron-hole pairs in a semi- 
conductor and thus generate power. The solar cells are 
connected in banks to provide the desired terminal voltage 
and current. 

A fuel cell-solarenergy power supply wherein incident quanta 
from the sun are used to dissociate the reaction products from 
a fuel cell. The reactants, after separation, are recycled to 
the fuel cell to produce electrical power. 

A thermoelectric power supply wherein a temperature differ- 
ence (Seebeck effect) is used to generate electrical power 
statically. The thermoelectric power for individual devices 
is extremely low; thus, hundreds of units operating at high 
temperature must be used to provide power at practical volt- 
ages. The heat source for these devices can be derived either 
from a nuclear reactor or solar collector system. 

A thermionic power supply wherein thermal energy is used 
to release electrons from a thermionic cathode. High cathode 
temperatures are utilized to obtain efficiencies comparable 
with other heat engines. Banks of thermionic generators are 
required for practical applications. The heat source for these 
devices can be derived either from a nuclear reactor or solar 
collector system. 


Open Cycle Chemically-Fueled Power Supply 


Open cycle power supplies use the hot exhaust gases from 
propellants to drive a prime mover, which provides shaft 
power for the generation of mechanical or electrical energy. 
Such units can use cryogenic propellants as well as the more 
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conventional mono- and bipropellant liquids and the solid j»ro- 
pellants. Fig. 2 is a block diagram of an open cycle sysiem 
which utilizes, for this example, cryogenic propellants. ‘he 
system consists of: Propellant tanks, combustion cham er, 
a turbine as a prime mover, an alternator and controls. 

A study of power supplies for long duration space vehi les 
indicates that, for an open cycle system, the minimum })ro- 
pellant weight is attained from the use of liquid hydrozen 
and liquid oxygen as the propellants. 

A simple open cycle system which utilizes cryogenic })ro- 
pellants may be pressurized by supplying heat to the pro- 
pellant tanks. The power into the heater is controlled as a 
function of the tank pressure with a pressure sensitive switch. 

The hydrogen leaving the storage tanks is used as a coolant 
to absorb the heat generated in the rotating equipment as 
a result of inefficiencies. The waste heat from the prime 
mover exhaust is used to raise the hydrogen temperature so 
that the hydrogen can be used to heat the oxygen. The hy- 
drogen and oxygen are heated to equal temperatures by pzss- 
ing these gases through a concurrent heat exchanger before 
entering the constant-pressure regulator; thus, gas is supplied 
to the reaction chamber at a fixed mixture ratio producing a 
2200 F hot gas to drive a turbine. The gas temperature is 
maintained below a stoichiometric temperature by utilizing a 
hydrogen-rich mixture of six times the stoichiometric require- 
ments. 

The gases, after being expanded through the nozzle, im- 
pinge on the buckets of a multistage single or multiple disk 
turbine. Although the general system described uses a 
turbine as a prime mover, substitution of a positive displace- 
ment device for a turbine would not materially affect the over- 
all concept or the performance. Table 2 represents a study 
of the characteristics of a 1-hp, 324-hr, open cycle prime 
mover with an optimum propellant. Studies indicate that, for 
short duration chemically-fueled open cycle systems, an 
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Fig. 2 Open cycle chemical power supply 
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optimum system would weigh approximately 600 lb, and the 
yolume would be approximately 60 ft’. The information in 
Table 2 has been scaled and presented in Fig. 15 to facilitate 
comparison with other systems. 


Closed Cycle Chemically-Fueled Power Supply 


The closed cyele system is very similar to the open cycle 
system described previously except that, in this case, the 
energy from the reaction is utilized to heat a working fluid 
whic! drives the prime mover. In the open cycle system, the 
hot gas products from the propellant were used to drive the 
prim: mover directly. For the closed cycle chemically- 
powered system, a cryogenic oxygen-hydrogen bipropellant 
isbe~t. In contrast to the open cycle system, the propellants 
are now burned stoichiometrically in order to obtain the maxi- 
mun: energy per pound from the reaction. 

Stress considerations limit the temperature of the prime 
mover to about 2000 F. The mass flow of the propellants 
to tie burner is controlled so that the turbine working 
fluid does not exceed this temperature. 

The closed cycle system which is represented in Fig. 3 
consists of: 

Oxvgen-hydrogen tanks which provide the chemical 
enerzy to the system. 

A bipropellant constant-pressure regulator to control the 
pressure to the gaseous propellants. 

A boiler which uses heat energy from the propellants to 
vaporize the working fluid. ; 

A turbine which transforms the energy available in the gas 
to mechanical power to drive an alternator. 

A condenser which radiates the unused heat in the cycle. 

A pump which returns the working fluid from the condenser 
to the boiler for recirculation in the system. 

‘The oxygen and hydrogen are supplied from their respective 
tanks by thermal pressurization. Tank pressure is 
controlled by varying the power to the tank heaters. The 
hydrogen gas is utilized as a coolant for the rotating machine 
so that its temperature can be maintained at a prescribed 
level. Both the gaseous oxygen and hydrogen after passing 
through concurrent heat exchangers are supplied to the bi- 
propellant pressure regulator and then to the reaction cham- 
ber portion of the boiler. The working fluid which is sup- 
plied to the boiler from the condenser is vaporized by the heat 
transferred from the reaction. The high pressure high tem- 
perature working fluid is expanded through the nozzle and 
impinges on the turbine. The superheat in the exhaust gas 
and the latent heat of condensation are radiated to the back- 
ground temperature of space. The resulting liquid working 
fluid at the outlet from the condenser is passed through an 
alternator cooling jacket and then to the boiler. 

The closed cycle system, in addition to the basic elements 
which have been considered, also requires controls for the 
liquid content of the boiler, the boiler pressure and the tur- 
bine speed. 

The selection of a working fluid for a Rankine engine de- 
pends on a number of factors. The most important of these 
are the temperatures to be utilized in the cycle. In order to 
achieve a high thermal efficiency, it is important that the 
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Fig. 4 Enthalpy-entropy diagram for sulphur 


critical temperature of the working fluid be at or above the 
maximum temperature to be utilized in the cycle. Superheat- 
ing above the critical temperature results in little gain in ther- 
mal efficiency but a very considerable increase in material 
problems. It is, of course, important that the fluid remain 
liquid down to the minimum temperature to be utilized in the 
cycle. An additional fluid property which is of considerable 
significance in system design and durability is the shape of the 
Mollier diagram. It is quite desirable to have a fluid on 
whose Mollier diagram the saturated vapor line shows a posi- 
tive slope. Comparison of a number of fluids, including 
water and mercury, led to the selection of sulphur as a closed 
cycle working fluid. It has a relatively low melting point, 
240 F; its critical temperature is high, approximately 1900 F, 
and its Mollier diagram indicates a positive slope of the 
saturated vapor curve up to approximately 80 atm as indi- 
cated in Fig. 4. 

A study program for a 324-hr, 1-hp, closed cycle system 


Table 2 Chemical open cycle system performance, O. and H: power supply 
Total energy, 324 hp hr 


Shaft power, 1.0 hp 


-——— Prime mover _ Propellant 
Type Weight, Il S.F.C., 
Type Ib/hp hr Weight, lb Volume, ft* Tank weight, Ib Total weight, lb 
‘our stage turbine 20 H2-O2 1.8 556 62.0 50 626 
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24,000 1.68 545 
64,000 1.33 433 


Table 3 Chemical closed cycle system performance, O. and H: power supply 


Rotational speed, rpm _ S.F.C., lb/hp hr Fuel weight, lb Fuel volume, ft* Tank weight, lb 


Prime mover Total system 
weight, lb weight, lb 
22.4 24 60 640 
17.8 20 37 508 


which utilized sulphur as the working fluid and hydrogen 
and oxygen as the propellants has been completed. Table 3 
indicates the weight and volume characteristics of this system. 

The table indicates that for the closed cycle compared with 
the open cycle system, the machinery weight is higher, but the 
superior performance results in a total system weight approxi- 
mately 10 per cent less. The volume, for this oxygen- 
hydrogen system operating stoichiometrically, is a factor 
of three less than the open cycle requirements. 


Closed Cycle Solar Power Supply 


The closed cycle solar power supply will be identical to 
the chemical closed cycle system except for the energy source. 
For the solar power system, a parabolic mirror is used to col- 
lect the incident solar energy and focus this energy into a 
boiler. The latent heat of condensation and superheat, in 
a manner identical to the chemical closed cycle system, is radi- 
ated to space (Fig. 5). 

The incident energy available to a collector in the vicinity 
of the Earth’s surface is 440 Btu/hr-ft?. This corresponds to 
approximately 0.13 kw/ft?. The energy available at the 
collector is less than 440 Btu/hr-ft? by the amount that is re- 
radiated by the boiler. The boiler access area, which in the 
limit is determined by the diameter of the sun’s image at the 
focus, in practical applications will be the area required to 
accommodate the various aberrations in the mirror system. 

The mirror size can be determined from the Carnot and 
device efficiency. The Carnot efficiency 


_Tu-Te 
TH 
is the maximum thermal efficiency obtainable from a heat 
engine working between a high temperature 7, and a low 
temperature Tp. Device efficiency is defined here as the ratio 
of the actual thermal efficiency of the device in question to 
that achievable by a Carnot engine. This ratio f is a useful 
means of comparing various heat engines. The value of f 
tends to remain constant for a given cycle over a wide range 
of operating points. Maximum thermal efficiency (product 
of Carnot and device efficiency) results in a minimum size 
collector. The Carnot efficiency is dependent on the boiler 
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I Pump 


Fig. 5 Closed cycle solar power supply 


temperature, which is limited by material considerations, «ind 
the radiator temperature. Low radiator temperatures, how- 
ever, result in large radiator area requirements, since radiant 
heat rejection is proportional to the fourth power of the r:di- 
ator temperature. 

A preliminary study was made to determine the optimum 
temperatures for the boiler and radiator. Expressions for the 
weight of the collector and radiator (as a function of tem) er- 
ature ratio, area, specific weights, device efficiency and ti-m- 
perature) were optimized for minimum weight. The resiilts 
of the study are presented in Figs. 6, 7 and 8. 

Fig. 6 permits the computation of a parameter R which  le- 
pends on 7'y, the high temperature in the cycle; W; and W, 
the specific weights of the collector and radiator; nx the ffi- 
ciency of the collector, and € the emissivity of the radiator. 
Inefficient collection yields high values of R; inefficient racia- 
tion yields low values of R. 

In Fig. 7 the optimum radiator temperature is presented. 
a, the ratio of Tp to Tz is given as function of R with the device 
efficiency f as a parameter. For high values of R (collecior 
size problems) the radiator optimizes at low temperatures, 
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Fig. 6 R’ vs. maximum cycle temperature for a power system 


with a solar collector, prime mover and radiator 
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Fig. 7 Optimum temperature ratio for a power system with a solar collector, prime mover and radiator 


which yields relatively high thermal efficiency. For low 
values of R, the radiator optimizes at high temperatures 
which tends to minimize radiator area. 

Fig. 8 presents the combined collector and radiator weight 
as a function of R and device efficiency. The weight is nor- 
malized here by the specific weight of the collector and the 
collector efficiency. 

An example of the use of the curves for a 100-kw system 


follows. Fig. 6 
Tu = 2240R 
R’ = 100 
W,/W. = 0.4 
€/nx = 1 
thus 
R = 40 
Fig. 7 
R = 40 
f = 0.68 
thus 
a = 0.465 
Tr = 0.465 X 2240 
Tr = 1030 
Fig. 8 
R =40 
=068 
= 28 
NK = 0.9 
= 0.2 


At cust 1959 


thus 


W = 6.5 lb/kw (collector + radiator) 
nru = f(1 — @) = 36.5 per cent 
Pix = 100kw/0.36 = 275 kw 
275kw 

= = 2370 ft2 
W. = 2370ft? x 0.2 = 474 Ib 
Wer 650 — 474 = 176 lb 


The results from the curves do not consider variations in 
the performance of the prime mover in the systems. An 
actual system would have to take into account the H-S 
curve for the working fluid as well as the relationship be- 
tween device efficiency and pressure ratio. 

The weights of closed cycle power systems, using the data 
of the collector-radiator study and including the weight of 
other equipment, such as turbines, alternators and boilers 
established in prior studies, are presented in Fig. 16 as a func- 
tion of power for time invariant systems. The same infor- 
mation is presented in Fig. 15 at powers of 1 and 100 kw, so 
that it can be compared with time variant systems. 

Practical considerations will limit the application of solar 
energy. Before a 100-kw device can be made feasible, it will 
be necessary to devise erecting schemes and mechanical con- 
figurations so that large, low weight parabolic mirrors can 
be unfolded in space. The mirror, for the 100-kw system 
for example, must be 60 ft in diameter. ; 


Closed Cycle Nuclear Power Supply 


The closed cycle nuclear power supply shown in Fig. 9 is 
similar to the closed cycle solar and closed cycle chemical 
power supplies in that it uses a Rankine cycle heat engine to 
convert thermal energy into electrical power. Heat is sup- 
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Fig. 8 Optimum radiator and collector per unit specific weight for a power system with a solar collector, prime 
mover and radiator 


plied to the working fluid from a nuclear reactor and is re- 
jected to outer space by radiation. The reactor for this ap- 
plication will be a small nuclear reactor, with or without 
shielding, depending on the mission for which the power supply 
is to be utilized. 

Current technology indicates that small nuclear reactors 
using sulphur as the working fluid can achieve a thermal output 
of approximately 2 megawatts (2000 kw) per ft? with a 
weight of approximately 300 lb per ft*, exclusive of shielding. 
The minimum reactor size which may be made critical is ap- 
proximately 3 ft?. Thus, for any space vehicle power supply 
using nuclear energy there will be a minimum of 1 megawatt 
(1000 kw) of thermal energy available from a }3-ft* reactor 
weighing 150 Ib. Using a cycle efficiency of 35 per cent ob- 
tainable from closed cycle sulphur systems, up to 350 kw can 
be provided by the use of this minimum 1 megawatt heat 
source. 
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Fig. 9 Closed cycle nuclear power supply 


For large power outputs, the weight of the radiator is the 
dominant factor in total system weight. A study was made 
to determine the criteria for minimizing the radiator area as a 
function of the maximum cycle temperature and the device 
efficiency. The thermal efficiency is given by the product of 
the Carnot and device efficiency. 

The results of this study show that for a device efficiency 
of 1.0, the minimum radiator area is achieved when the ratio 
of radiator to maximum temperature is 0.75 (see Figs. 10 and 
11). The thermal efficiency for this condition is 0.25. For 
device efficiencies approaching zero, the optimum temper- 
ature ratio approaches 0.8 and the thermal efficiency ap- 


proaches 0.2f. An example of the use of the curves for a 400- 
kw system is, from Fig. 10 
Py = 400 kw 
= 0.8 
7.26 
24 
4.22 
25 
.18 
H 79 16 
TH 4.14 
TH- TR 
¢= DEVICE EFFICIENCY 10 
Tr=RADIATOR TEMP (°R) We 7.08 
Ty=MAX CYCLE TEMP (°R) dos 
76+ 4.04 
4.02 
0 2 4 6 8 10 


DEVICE EFFICIENCY -f 


Fig. 10 Optimum temperature ratio and thermal efficiency for a 
power system with an internal heat source, prime mover and 
radiator 
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Fig. 16 shows the power vs. weight for the unshielded 
nuclear closed cycle system. For power levels less than 350 
kw, the weight increases due to increases in the prime mover 
and radiator weight. For power levels in excess of 350 kw, 
the reactor size as well as the prime mover increases propor- 
tions lly to the output. The comparison with time variant 
syst-ms is shown in Fig. 15. 

Specific weights for the radiator calculations were assumed 
to |» 2.0 lb per ft? for a power level of 1 kw, and 0.5 Ib per 
ft? r power levels in excess of 100 kw. 


Solar Cell Power Supply 


Solar cells have recently been used as direct converters of 
solar energy for power applications. These cells have been 
introduced commercially by a number of companies. 

‘The semiconductor solar cells consist of P-type and N-type 
silicon. The N-type silicon has an excess electron due to the 
doping of the silicon with a material which has five electrons 
in the valenee orbit. The P silicon has a lack of electron 
or : hole due to the doping of this part of the cell with a ma- 
terial with three electrons in the valence orbit. When photons 
from the sun impinge upon the surface of the solar converter 
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Fig. 11 Optimum radiator area for a power system with an in- 
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they are absorbed with a resultant displacement of electron- 
hole pairs. The displacement of an electron or a hole near 
the junction will result in a continuous shift of electrons and 
holes within the complete lattice and through the load, thus 
providing electrical current. 

Solar converters presently provide output with a conversion 
efficiency of approximately 10 per cent. Thus for a 100-w 
solar conversion system, approximately 8 ft? of silicon cells 
will be necessary. The figure of merit for solar cells has been 
reported to vary from 0.1 to 0.25 lb per w. The lower 
value is used for large multikilowatt converters. Current 
research is directed at increasing the conversion efficiency of 
solar cells by variation of composition and stacking. This 
latter technique involves the interception of light passing 
through one cell by a second cell placed behind it which is 
sensitive to wave lengths to which the surface cells are trans- 
parent. These gains in efficiency will, of course, be accom- 
panied by an increase in weight per unit area and a reduction 
in area requirements. Series parallel connections are neces- 
sary to provide the required voltage and current from solar 
cells. For a 3-kw system, approximately 100,000 cells with 
a total weight of 600 lb will be required. 

Solar cells have an advantage over the systems discussed 
previously, in that they can be assembled into systems of 
very low power level with little fixed weight penalty. Thus 
a 10-w system will require approximately two times the 
weight of a 5-w system. In contrast, the weight of turbo- 
machinery at low power levels tends to approach an irreducible 
minimum independent of power level, so that a 10-w solar 
closed cycle system would probably weigh substantially 
the same as a 50- or a 5-w system. 

For systems which operate without incident solar energy, 
it is necessary to utilize electrochemical batteries in con- 
junction with the solar cell. The batteries for this type of a 
configuration will be charged during the period when the solar 
collector can be oriented toward the sun and discharged dur- 
ing the dark part of the orbit. 

For alternating current requirements it will be necessary 
to use either a rotary-type inverter or a static transistorized 
inverter. 

The regime wherein the solar cell configuration is optimum 
is shown in Figs. 16 and 18. 


Fuel Cell-Solar Energy Power Supplies 


The fuel cell-solar energy system makes use of two sources 
of energy. These are the exothermic energy available from 
the reaction of the stored reactants of the fuel cell system and 
the energy available as radiant energy to reconstitute the 
product of the reaction to the original form. 

The fuel cell does not differ from a conventional battery in 
its theory of operation. The difference is physical; in a con- 
ventional battery the chemicals are stored on the plate ma- 
terials, whereas in the fuel cell the chemicals are stored ex- 
ternally in separate containers. 

A typical space vehicle, one possible application for a fuel 
cell, might require power for only a small part of each day 
spent in space. Thus, in a conventional cell, the reac- 
tants would have to be available for the total expected 
power-time product delivered to the end item (a transmitter, 
television camera, space surveillance equipment or other). 

However, if solar energy is utilized, the large power pulse 
can be furnished by the reaction of the stored reactants. 
The dissociation of the product, charging the fuel cell system, 
is accomplished over a longer period of time with a com- 
paratively small solar collector. The block diagram shown 
in Fig. 12 represents a typical solar energy charging system. 
Proper mechanization of the system will result in the overall 
conversion of solar energy into electricity with an electro- 
chemical process as an intermediate step. 

The solar conversion scheme is in the research stage, and, 
as such, no information is available pertaining to the weight, 
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Fig. 12 Diagram of recycling fuel cell using solar energy to 
photo-dissociate product of fuel cell 


volume or, indeed, the practical feasibility of such a device. 

In order for the fuel cell-solar conversion scheme to be prac- 
tical and competitive, it will be necessary for this device to 
make maximum use of the available energy from the sun. 
The silicon solar cell presently provides 100 w per m?, and 
the closed cycle solar-turbines provide approximately 350 w 
per m? of collector area. The fuel cell-solar system must 
provide more than 100 w per m? in order to compete with the 
simple solar cell. 

Conventional photolysis techniques for dissociating the 
reaction products are effective for wave lengths less than 
3100 A. The incident energy, however, to 3100 A is only 
27 w perm. If the photolysis range is extended to 4100 A and 
thence to 5100 A by the use of sensitizers, the incident power 
increases to 134 and 330 w per m?, respectively. 

A type of fuel cell which has received considerable research 
attention for the past 100 years uses the hydrogen-oxygen reac- 
tion. It should be emphasized that hydrogen and oxygen 
reactants provide an electrochemical device with minimum 
weight and volume from a stored chemical viewpoint. A 
number of other possibilities exist, such as a nitric oxide- 
chlorine reaction or a lithium-hydrogen reaction which, from 
a system point of view, might provide a better solar regener- 
ated fuel cell end item. 

In particular, the nitric oxide-chlorine reaction has a much 
wider spectral sensitivity for photochemical dissociation than 
hydrogen-oxygen and appears very attractive. The system 
would consist of a cell utilizing the reaction product nitrosyl 
chloride as a solvent for a suitable electrolyte. Overall con- 
version efficiency of the nitric oxide-chlorine solar regenerative 
fuel cell will be comparable to silicon solar cells. Advantages 
of this system include the ease with which large systems can 
be fabricated together with the inherent storage capabilities 
of the system. 

The fuel cell differs in one major respect from the other 
closed or open cycle systems. Each of the other systems is 
limited by the Carnot efficiency. The fuel cell per se is not 
Carnot-limited and theoretically can achieve 100 per cent 
efficiency. Fuel cells using hydrogen-oxygen have been re- 
ported with efficiencies of 60 to 70 per cent. The major 
disadvantage is the weight. Present estimates range from 
130 to 400 lb per kw for 600 psi, 240 C and 15 psi, 50 C 
batteries. 


Thermoelectric Power Supply 


Recent advances in the technology of materials have made 
feasible the generation of power by thermoelectric methods. 
P- and N-type materials compounded from high atomic 
weight substances, such as lead, bismuth, tellurium and mer- 
cury have been investigated. Research studies indicate 
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that theoretical device efficiencies of 5 to 7 per cent for bis- 
muth telluride-bismuth and P- and N-type bismuth telluride 
are obtainable. 

Research activity has been directed toward materials which 
provide high thermoelectric power (mv/deg C) with minimum 
resistivity and thermal conductivity. It is desired to opti- 
mize these parameters in order to maximize the generated 
voltage and to minimize the internal power (J?R) and therinal 
conductivity losses. 

Two schemes have been considered which use a theriio- 
electric generator for space application. One method uses a 
nuclear reactor to provide the thermal energy and a radiator 
to maintain the temperature of the cold junction. ‘lhe 
second method uses a solar mirror to focus the incident energy 
on to the junction. These systems are presented in Figs. 13 
and 14. 

The efficiency for these heat engines, like the closed cycle 
turbine or reciprocator, can be expressed as the product of 
the device and Carnot efficiencies 


Tn — Te 


TH = 


Expressions for f have been derived which involve the th«r- 
moelectric power of the materials used S, the product of the 
resistivity p and the thermal conductivity k 


M + (Te/Tx) 


v(TatTr 
5/ Tr 


Sit 
(V + kop»)? 


Static thermoelectric engines are presently limited in tem- 
perature difference as well as maximum temperature. The 
reported thermal efficiency for laboratory devices operating 
at 2000 F with a AT of 990 deg is approximately 2 per cent. 
The output voltage for these devices is approximately 0.2 mv 
per deg C. One hundred and eighty-seven couples will be 
required for a 20-v bank with a AT of 990 F. The devices, in 
addition, may require pumps and motors to circulate the work- 
ing fluid and coolant, and thus approach a rotating closed 
cycle heat engine in complexity without providing the per- 
formance. The reliability of static devices may be low due 
to the requirements for large numbers of high temperature 
connections. 

The parameter R (Fig. 6) for a hot junction temperature 
of 2460 R, a radiator weight of 2.0 lb per ft”, an emissivity of 0.9 
and a collector at 0.2 lb per ft? with 90 per cent efficiency is 
equal to 15. An examination of Fig. 8 for a device efficiency 
of 5 per cent indicates W’ = 600, which yields, for the col- 
lector and radiator, a specific weight of 133 lb per kw. The 
specific weight for the generator is reported to be approxi- 
mately 1000 lb per kw; accordingly, the system would weigh 
1133 lb to produce 1 kw. 

A somewhat different set of operating conditions is re- 
quired if the static thermoelectric generator is used with a 
nuclear heat source. There will, in general, be a temperature 
limitation to the reactor, say for example 1500 F. For this 
temperature and an efficiency equal to 5 per cent of the Carnot 
efficiency (f = 0.05), a radiator area of 36 ft? is required at 
an emissivity of 1.0 (see Fig. 11). Assuming e = 0.9, 
40 ft? will be needed. At 2 lb per ft?, the radiator weight will 
be 80 lb. The minimum size nuclear heat source is about a 
3-ft® device and will weigh about 150 lb. The total system 
weight including auxiliary equipment, will be approximately 
1230 lb for a 1-kw system. 
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Fig. 13 Thermoelectric solar power supply 
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Fig. 14 Thermoelectric nuclear power supply 


Thermionic Power Supply 


The thermionic generator is similar to the thermoelectric 
device except that a thermionic cathode is used instead 
of P-N junctions. Physically, the thermionic device 
consists of an anode and cathode with finite spacing. The 
potential from the device is dependent on the thermochemical 
potential, the thermoelectric power and the voltage drop due 
to lead resistance. The current is a function of the temper- 
ature, the work function for the materials and the space 
charge. 

Efficiency expressions have been derived for these devices 
which are identical to the expressions for the thermoelectric 
engine for the limiting case wherein AT’ approaches zero. 
The thermal efficiency may again be expressed as a product 
of the thermal efficiency of a Carnot engine and a device 
efficiency 


nrH = Nef 
The device efficiency is in turn a function of a parameter Z 
(VZ+1-1) 


where Z is a function of the open circuit voltage, the temper- 
ature and the Boltzmann constant. 

Laboratory evaluation indicates that thermal efficiencies 
of 13 per cent can be attained with cathode temperatures of 
2300 F and anode temperatures of 1000 F. The device effi- 
ciency is approximately 28 per cent. The terminal voltage 
at this optimum efficiency (for this configuration) was 0.5 v. 
Che current density was approximately 1 amp per cm?. 
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A 1-kw unit, considering the present state of the art, would 
thus require approximately 54 series units with 50 units in 
parallel for a total of 2700 units (assuming 1-cm? cathodes and 
24 v). Space applications would be similar to the thermo- 
electric generator (Figs. 13 and 14). No realistic weight 
estimates for a practical thermionic device are available. 


System Comparisons 


The characteristics of some of the various systems de- 
scribed in the preceding sections are tabulated in Table 4. 
In this table 1- and 100-kw systems have been considered 
with estimates for the weight and the size of the power 
supplies and their significant components. All the signifi- 
cant assumptions are tabulated in the notes listed. Sev- 
eral systems which have been discussed in the text have been 
omitted from the tabulation and the subsequent comparisons 
in Figs. 15 and 16. The thermionic systems and the solar 
regenerative fuel cell systems were omitted due to the un- 
availability of weight estimates for practical systems at the 
present time. The thermoelectric systems have been omitted 
because the available weight estimate of 1000 lb per kw fora 
thermocouple device seems unreasonably high, and no more 
realistic weight estimates were available. 

Several important features of Table 4 are discussed below. 
It will be observed, for example, that the turbine weight does 
not change in going from a 1-kw system to a 100-kw system. 
Virtually the same turbine would be used with the difference 
in power level resulting from a combination of increased are 
of admission and higher back pressure. In the closed cycle 
systems, this higher back pressure corresponds to a higher 
condenser temperature, which is indicated in the notes. Thus 
an increase in exhaust density by a factor of 10 and an in- 
crease in the are of admission from 10 to 100 per cent yields 
roughly a 100 to 1 increase in the power generated. 

Another item noted in Table 4 which needs a word of ex- 
planation is the alternator weight assumed. Although the 
alternator never represents a large fraction of the total sys- 
tem weight, it should be noted that a decrease in specific 
weight from 2 to 0.3 lb per kw has been assumed as the power 
increases from 1 to 100 kw. Recent experimental work 
demonstrated a 22-kw alternator weighing 11 Ib. This was 
not a fully developed device, but it is certainly indicative of 
the future trends in rotating electromagnetic machinery. 
Electrostatic generators have been reported in the literature 
with similar weight estimates. 

The radiator temperatures, sizes and weights do not corre- 
spond to the optimum temperatures indicated in the analyti- 
cal studies presented in Figs. 6, 7, 8, 10 and 11. The differ- 
ences arise from practical considerations of the vapor pressure 
characteristics of the working fluid and the effect of back 
pressure on turbine performance. 

The information presented in Table 4 together with suit- 
able interpolations and extrapolations are presented in Figs. 
15 and 16. In Fig. 15 a comparison of system weights as a 
function of time is presented. The question at issue here is 
whether or not there is any application for the chemically- 
fueled space vehicle power supply. We find that the chemical 
open cycle system is lighter than all other systems at the 1-kw 
level up to about 10-hr duration, at which point the solar 
closed cycle powerplant takes over. For durations longer 
than 10 hr, this solar powerplant is optimum from a weight 
viewpoint. At the 100-kw level, the chemical open cycle 
powerplant maintains a weight superiority up to approxi- 
mately 2 hr where the unshielded nuclear closed cycle device 
becomes the optimum powerplant for longer durations. It 
should be noted that the chemical closed cycle machine is 
superior to the chemical open cycle machine for durations in 
excess of 35 hr at the 1-kw level, but for this duration the solar 
powerplant is lighter than either one of them. 


3 At Sundstrand Turbo Division. 
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Table 4 Power supply characteristics 


TYPE CHEMICAL CHEMICAL SOLAR NUCLEAR 
OPEN CYCLE CLOSED CYCLE CLOSED CYCLE | _ CLOSED CYCLE 
; Energy Source Liq. Hz - Liq. 02 Liq. Hz - Liq. 02 Solar Unshielded 
Nuclear Reactor 
Power Level 1 Kw 100 Kw 1 Kw 100 Kw 1 Kw 100 Kw, 1 Kw 100 Kw 
WEIGHT 
¢ Energy Source 2.8 #/Hr 170 #/Hr 2.1 #/Hr 160 #/Hr 6# 480# 150# 150# 
Prime Mover 
Turbine 208 208 204 20# 20# 208 208 20# 
Ex -- -- 5# 20% 5# 100# -- 
Alternator 24 30# 2# 308 2¢ 308 30# 
Radiator -- -- 18# 170# 18# 170# 24# 170# 
: TOTAL WEIGHT 22# + 2.8 #/Hr| 50# + 170 #/Hry 45# + 2.1 #/Hr 240#+ 160#/H 51# 800# 196# 370# 
Energy Source 0.32 ft?/Hr | 19.2 £t°?/Hr 0.092 5.65 30 2400 ft 0.5 0.5 
Prime Mover 
Turbine 200 in.? 200 in.? 200 in.? 200 in.? 200 in.f 200 200 200 in.? 
H Ex 250 875 in.? 500 in.? 1 15 in.7] 520 in. = = -- -- 
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h, = 48% = 802 = 60% nN, = 80% = h, = 802 M, = 602 = 802 | 4 watts 10 vatte 
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NOTES Kw Kw Kw Kw Kw Kw Kw ne Q = 10% 
Neie = 85% Ty = 1780°F | Ty = 1780°F | Ty=1780°R| Ty = [Ty=1600F |Ty=1600°F 
To = 400°F | T, = 540°F 400°F| 540°F |7,,=540°F 
4OL Ry - 60% 07 by weight c Cc Cc Cc Cc 
a rad rad rad a fz 
Pressure Pressure 
= = 85 = 
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Fig. 15 Comparison of space vehicle power systems 


There are many considerations which might alter the exact 
crossover points for these various powerplants; for example, 
no provision has been made in the solar powerplant for heat 
storage on the dark side of an orbit. Similarly, no provision 
for shielding has been made for the nuclear powerplants. 
Such considerations will undoubtedly vary the crossover 
time and perhaps the powerplant selection for particular 
applications. The significance to be drawn from Fig. 15 is 
that, at least from a weight point of view, the open cycle 
chemical power supply has a place in the sun, but only for 
relatively short duration missions. For those missions where 
erecting a mirror or carrying a nuclear powerplant are not per- 
mitted, the choice between the open and closed cycle power- 
plants will depend on the relative importance of powerplant 
weight and size. There is roughly a 3 to 1 propellant-volume 
ratio in favor of the closed cycle machine, which may be 
significant for the longer durations. 

The variation of system weight with power level for those 
systems whose weight is essentially independent of time is 
presented in Fig. 16. Here again, the merit in such a plot 
lies in the estimation of approximate power levels at which 
the weight crossovers occur. Allowing for the uncertainty 
in the present state of the art in solar ceils together with 
some optimism in the closed cycle solar powerplant weight, 
we would estimate a crossover between these two at about 
100 w. The crossover between the closed cycle solar power- 


plant and the closed cycle nuclear powerplant occurs at ap- 
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Fig. 16 Space vehicle power systems weight vs. power 
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proximately 20 kw. In addition to these crossover points, it 
is rather interesting to observe the specific weight values 
which these powerplants approach at high power levels. The 
solar cells, as indicated in Fig. 16, should achieve a specific 
weight of about 100 lb per kw at large power levels. Some of 
the more optimistic investigators in this area have estimated 
that they will achieve a 50-lb-per-kw specific weight for large 
power levels, but this state of the art advance is, to the best 
of our knowledge, not possible at the present time. The 
closed cycle solar powerplants should achieve a specific weight 
of approximately 10 lb per kw at large power levels; the closed 
cycle unshielded nuclear powerplants should achieve approxi- 
mately 2 lb per kw at about a 1-megawatt rating. 


Conclusions 


The powerplant selection indicated by Figs. 15 and 16 is 
presented graphically in Fig. 17 together with an estimate as 
to when these powerplants will be available. An estimate of 
the crossover point between batteries and solar cells at low 
power levels has also been included for the sake of complete- 
ness. We see that below a power level of 100 w, batteries 
are optimum up to about 20-hr duration, at which point solar 
cells become the lightest available system. Both of these 
devices are currently available. From 100 w to 20 kw, chem- 
ical powerplants are optimum for relatively short durations, 
say up to 10 hr, whereas for the longer durations a solar 
closed cycle machine is optimum. As indicated in the graph, 
however, the solar powerplant is not presently available, 
although we estimate that a working machine would be 
available by approximately 1961. Chemical powerplants 
are the proper choice for interim applications. 

For power levels in excess of 20 kw, chemical powerplants 
are suitable for short durations on the order of 5 hr or less, 
depending on power level, whereas nuclear powerplants, which 
should be available by 1964, are optimum for the longer dur- 
ation application. Missions requiring 20 kw or more will use 
chemical powerplants until more advanced systems become 
available. 

An important conclusion which we have reached as a result 
of the study presented here is indicated by Fig. 18. Here 
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we show a closed cycle powerplant which can derive its heat 
energy from either chemical, nuclear or solar means. Each 
of these powerplants is well suited to practical missions in 
space. The superiority of these powerplants over other 
means of power generation which have been considered stems 
largely from the high thermal efficiency which is achievable 
from Rankine cycle engines. In particular, at high power 
levels where total system weight is directly related to thermal 
efficiency, these closed cycle systems are clearly superior. 
Turbomachinery will find application in space vehicle power 
systems for a long time to come. 
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Lockheed X-17 Rocket Test Vehicle 


RONALD SMELT' 


and Its Applications 


Lockheed Aircraft Corp. 
Sunnyvale, Calif. 


The Lockheed X-17 test vehicle was originally designed, in 1955, to permit experiments to be made 
in flight at the Mach numbers and Reynolds numbers appropriate to re-entry of ballistic missile 
nose cones. It consists of three stages of solid propellant rocket motors, and carries an instru- 
mented model re-entry body. Although the model is much smaller in scale than the simulated 
re-entry body, full-scale Reynolds numbers at Mach numbers up to 15 are obtained by a flight pro- 
gram in which the second and third stages are fired on the downward part of the flight trajectory. 
Details of design and performance, and a description of the telemetering equipment developed for 
the vehicle are given in the report. No guidance equipment is carried; the X-17 depends upon aero- 
dynamic stabilization for correct orientation, both during the ascent and in the re-entry phase. 
The resulting oscillations and dispersion have been analyzed in considerable detail, and the results 
of the analysis are summarized. Although most of the flights of the X-17 have been made to obtain 
re-entry data, some recent flights have utilized the vehicle in nuclear experiments at altitudes ex- 
ceeding 400 miles. The vehicle then operates as a conventional sounding rocket, all three stages 
firing in rapid succession during ascent. The performance of the X-17 vehicle operating in this 
manner is included in the report, together with some discussion of the techniques adopted to sta- 


bilize at high altitude by spinning. 


HEN the Air Force began the development of the present 

generation of long range and intermediate range ballistic 
missiles, it was immediately apparent that one of the major 
problems to be solved was the design of the vehicle to with- 
stand the conditions of re-entry into the Earth’s atmosphere 
near the end of its flight. The aerodynamic characteristics 
of this phase of the flight were unprecedented: Mach number 
in the range from 15 to 20, air stagnation temperatures around 
20,000 F, and excessive rates of heat transfer into the vehicle, 
indicated by theoretical work to rise to several thousand Btu 
per sec per ft? of surface. None of the conventional experi- 
mental aerodynamic facilities in the country could approach 
the simulation of these conditions, and although new types of 
facilities with extremely short flow durations, the shock tube 
and are-heated wind tunnels, showed some promise of being 
able to duplicate re-entry conditions, they were in a very early 
stage of development. They had not yet passed the acid test 
of aerodynamic facilities, the attainment of some measure of 
agreement with aerodynamic data from free flight. The re- 
entry problem presented a particular difficulty here, since the 
critical quantity, the heat transfer rate, is greatly dependent 
upon the point of transition from laminar to turbulent flow, 
and even at lower speeds in wind tunnels factors, such as 
stream turbulence and model surface roughness make it very 
difficult to simulate free flight transition conditions with cer- 
tainty. 

The Lockheed X-17 hypersonic test vehicle was developed, 
at the request of the Air Force Ballistic Missile Division, to 
give design information on this critical problem. Propelled 
by solid fuel rockets, and much smaller in scale than the 
ballistic missiles whose re-entry it simulated, it could be put 
into operation much more cheaply and rapidly than the full- 
scale vehicle. Design of the X-17 vehicle began in February 
1955, and in July 1956 it gave the first flight data on heat 
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transfer rate at the Mach numbers and Reynolds numbers of 
typical re-entry. (The technique used to attain full-scale 
Reynolds numbers, in spite of the smaller scale of the X-17 
vehicle, will be described.) Since that time, 26 X-17 vehicles 
have been launched on the Air Force re-entry program, and 
20 of these have been successful in providing aerodynamic 
data. 

Because of its unique features of hypersonic performance 
and excellent reliability, the X-17 test vehicle has been utilized 
on other programs subsequent to the Air Force re-entry in- 
vestigation for which it was first developed. Some of these 
later vehicles have been operated as re-entry vehicles, follow- 
ing a flight path similar to those of the original program. The 
X-17 has also been utilized in the Argus program for nuclear 
research at extremely high altitude, taking advantage of the 
fact that a high speed vehicle of this type is necessarily also 
capable of flight to very great heights. More details of the 
performance of the X-17 in these alternative applications are 
given later. 


Description of the X-17 Test Vehicle 


Vehicle Configuration 


The X-17, shown in Figs. 1, 2 and 3, is a three-stage tan- 
dem configuration, powered by solid propellant rocket mo- 
tors. Separation of the second stage is achieved by an ex- 
plosive charge usually fired by a pressure switch. Thrust 
from the third-stage rocket motor separates the third stage 
approximately 3 sec after second-stage separation. 

The rocket motors for all three stages are manufactured 
by the Thiokol Chemical Corp. An XM-20 (Sergeant) 
motor is used for first-stage boost, a cluster of three XM-19 
(Recruit) motors for the second stage, and an XM-19F1 
motor for third stage. Powerplant characteristics are listed 
in Table 1, and typical static test thrust curves are presented 
in Figs. 4 and 5. The variation of weight during rocket 
motor burning is given for each stage in Fig. 6. 
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Fig. la Lockheed X-17 on launching pad 


Fig. lb Launching of Lockheed X-17 


Two different methods of aerodynamic stabilization «are 
used on this vehicle. The first stage is stabilized by lifting 
surfaces consisting of four straight, tapered fins, with a 4 per 
cent thick biconvex section, arranged in cruciform. The 
second and third stages are stabilized by conical skirts with 
approximately 10-deg semivertex angles. 

In addition to aerodynamic stabilization, dispersion is 
minimized by spinning the first stage up to about 2 »ps 
immediately upon leaving the ground, using two small 
spin rockets clamped to the outside of the case of ‘he 
XM-20 rocket. These are jettisoned after operation, «nd 
in configurations designed for re-entry testing, the large fi: st- 
stage fins then take over, gradually reducing the spin to 21 
These first-stage fins are made large enough to provide raid 
damping of the large amplitude oscillations in angle of att: ck 
which occur when the vehicle re-enters the atmosphere, so as 
to ensure that the second stage separates and fires w! en 
aligned along the re-entry flight path. 


Table 1 Rocket motor characteristics 
XM-20 XM-19 XM-19E1 
Performance, sea level 
web burning time, 
70 F, see 23.8 1.53 1.53 
average thrust, 70 F, 
lb 50,000 34,861 35,950 
total impulse, 70 F, 
Ib-see 1,344,000 57,024 58,800 
total fimpulse/ 
initial weight, sec 164 164 169 
specific impulse, sec 191 216 223 
Weight, lb 
768 66.5 66.5 
nozzle 295 15.6 15.6 
liner 89 2.6 2.6 
propellant 7033 264 264 
total weight 8185 348.7 348.7 
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Fig. 2 Lockheed X-17 typical configuration 
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Fig. 3 Lockheed X-17 equipment location 


In configurations intended for high altitude applications, 
dispersion due to misalignment of second- and third-stage 
motor thrust, or to upper atmosphere winds, is minimized by 
spinning the combined second and third stages immediately 
after separation from the first. This is achieved by the addi- 
tion of small tabs at suitable points around the perimeter of 
the second-stage rocket nozzles, dimensioned to impart a 
spin rate of about 8 rps. A more detailed consideration of the 
whole problem of stabilization and dispersion is given in the 
Appendix. 


Equipment 


The majority of the flights of the X-17 vehicle have taken 
place from the Air Force missile range at Cape Canaveral, 
and the equipment normally provided in the first stage, as in- 
dicated in Fig. 3, is defined by the requirements of that range 
for missile tracking and destruction in the event of malfunc- 
tion. 

By far the most critical equipment item is the telemetering 
system carried in the final stage to transmit test data to 
ground. Evidently, since it is a major part of the third stage, 
its weight and bulk are of extreme importance. Further, it is 
required to withstand extremely severe acceleration condi- 
tions. These will be evident from the data given in Figs. 5 
and 6, where it will be observed that during third-stage burn- 
ing, a thrust of about 35,000 lb is acting upon a total mass 
which is decreasing from 500 to 230 lb, an acceleration of the 
order of 100 g. During re-entry experiments, conditions are 
still worse; as can be seen from a typical trajectory, such as 
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Fig. 4 Thrust characteristics of XM-20 rocket 
motor at 70 deg, sea level 
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Fig. 5 Thrust characteristics of XM-19 rocket 
motor at 70 deg, sea level 
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Fig. 6 Weight variation during burning 


in Fig. 11, the high drag associated with re-entry appears 
toward the end of third-stage burning and imposes a de- 
celeration of about 100 g immediately after an acceleration of 
the same order. When it is remembered that the thrust of 
the rocket motor is not completely steady, so that its varia- 
tions represent considerable shaking superposed on the 
high steady acceleration, the severity of the environment will 
be apparent. 
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Fig. 7 Lockheed X-17 telemetering transmit- 
ter assembly, original form 


Fig. 9 Lockheed transistorized FM-FM tele- 
metering system, typical component 


To meet these conditions, Lockheed Missiles and Space 
Division found it necessary to develop its own special tele- 
metering equipment. Fig. 7 shows the FM-FM system as 
used originally in the vehicle. It provides 16 continuous 
channels of information, coded in accordance with the 
standards set by the Inter-Range Instrumentation Group; 
one or two of the 16 channels are subcommutated, giving up 
to 48 channels of sampled data. Its total weight, including a 
10-w RF section, is 44 lb. Before use, it is tested on the 
ground in a centrifuge, to assure that it can withstand the 
accelerations and decelerations discussed, and it is also s\b- 
jected to shaking at accelerations up to 10 g, over a wide ‘re- 
quency range up to 2,000 eps. 

In contrast to Fig. 7, a section of the nose cone recovered 
from the desert at Holloman Air Force Base is shown in Fig. 8. 
Evidently the telemetering equipment is not yet able to w::h- 
stand the final decelerations upon impact after such a flig/it! 

More recently, the steady improvement in transistor te: h- 
niques has permitted the development of telemetering equ ip- 
ment which can meet even more severe environmental 
specifications than the original X-17 equipment, and is capa’ le 
of being packaged into much smaller bulk and weight. Fig. 9 
shows a typical component of the fully-transistorized tc\e- 
metering equipment; it permits the same amount of inforn.a- 
tion to be transmitted as the original system, with about 40 
per cent of the total weight and a considerable saving in 
battery power. 

It will be appreciated that all the instruments contained in 
the third stage of the vehicle—measuring pressures, tempera- 
tures, heat transfer rates and other quantities required to de- 
termine the aerodynamic characteristics—are subject to the 
same difficulties due to high accelerations. In addition, the 
instruments on the surface near the nose are subject to the ex- 
cessive heating rates during re-entry; temperatures up to tlie 
melting point of copper have been recorded. Again, it was 
necessary to develop special end instruments for many of 
these measurements. 


Performance of the X-17 Vehicle 


Performance as a Re-Entry Vehicle 


In the origina] application of the X-17 test vehicle to the 
aerodynamic problems of re-entry, the basic objective was to 
furnish information at Reynolds numbers appropriate to a 
full-scale re-entry body, although using a smaller scale model 
on the nose of the X-17 vehicle. Evidently this requirement 
can be met if the smaller scale of the X-17 model is compen- 
sated by making measurements at a lower altitude. The 
major difficulty is then to impart to the model, at this lower 
altitude, the hypersonic velocities anticipated in the full-scale 
re-entry body at its higher critical altitude; to achieve this 
in free fall would imply that the drag-to-weight ratio of the 
model must be much less than the full-scale object. Since this 
is extremely difficult to achieve in practice, the alternative 
expedient is employed of canceling the model drag by an ap- 
plied thrust during re-entry. 

To achieve this objective, the X-17 vehicle utilizes only 
the first stage to reach its maximum altitude, usually some- 
what in excess of 500,000 ft when launched at an angle of 85 
deg. The first-stage motor burns out at about 90,000 ft, 
when the velocity of the vehicle is slightly over 5,000 fps, and 
the vehicle coasts to apogee and into the descent. Typical 
trajectories during this phase of the flight are shown in Fig. 
10. Separation of the spent first stage from the remainder of 
the vehicle is intentionally not initiated during this part of the 
flight so as to make use of the high stability margin of the 
complete configuration with its large fins to align the vehicle 
along the descending flight path before the air density is 
sufficient to produce high heating rates. 

At a point on the descent defined by the setting of a timer, 
the first stage is separated, and the second stage fires, followed 
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Fig. 10 X-17 trajectories, first stage 


in about 3 see by the third stage. This drives the velocity 
rapidly up to hypersonic values. The maximum Mach num- 
ber and maximum Reynolds number attained during this 
portion of the descent can be varied over quite a wide range 
by selection of the altitude at which second-stage ignition is 
initiated. Evidently, if ignition is deferred until a low alti- 
tule is reached, the resulting maximum Reynolds number can 
be very high, but the maximum Mach number is not so great 
as that given by ignition at high altitudes. Fig. 11 compares 
the performances which can be obtained by typical high al- 
titude and low altitude second-stage ignition; the Reynolds 
number in this figure is a typical value utilizing the diameter 
of the nose of the vehicle as a characteristic length. 

This technique, appropriately referred to as the ‘nail- 
driving’ method, has proved to be very successful in provid- 
ing data at extremely high Reynolds numbers. It will be 
observed from Fig. 11 that the maximum Reynolds number 
obtained can be as high as 20 million for a flight with delayed 
second-stage ignition, and is around 8 to 10 million at the 
highest Mach numbers in flights aimed specifically at achiev- 
ing the highest speeds. Of course, the major importance of 
these high Reynolds numbers lies in the fact that they permit 
the test to cover the complete range of Reynolds number in 
which transition phenomena are normally expected to occur, 
so that the aerodynamics of both laminar and turbulent con- 
ditions can be studied. 


Utilization of the X-17 Vehicle for Maximum Altitude 


An alternative type of operation permits the X-17 test 
vehicle to be used as a sounding rocket attaining very high 
altitude, and a number of flights have been made in this 
manner. To obtain maximum altitude, the vehicle is again 
launched into an almost vertical flight path, but the second 
stage is now fired after burnout of the first stage, with the 
minimum possible delay to permit satisfactory separation. 
The third stage again fires as soon as possible after burnout 
of the second stage. The performance of the vehicle in this 
mode of operation, defined by the maximum altitude which 
can be reached, is of course entirely dependent upon the 
weight of instrumentation and equipment to be carried to the 
extreme altitude. Fig. 12 shows the calculated maximum 
height as a function of the total weight to be carried with the 
third stage, but excluding the weight of the third-stage rocket 
motor itself. The flights already made with the X-17 vehicle 
in this condition conform quite closely to the calculated curve 
on Fig. 12. 


Applications 


The large number of flights made in the X-17 re-entry re- 
search program have given aerodynamic information for a 
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Fig. 11 Variations in Mach number and Reynolds number 


very wide variety of re-entry vehicle nose shapes, ranging from 
simple hemispheres to the more complex shapes selected in 
the ballistic missile programs. Pressure distributions meas- 
ured on these shapes during X-17 flights, at Mach numbers up 
to 14 or 15, have shown excellent agreement with current 
hypersonic blunt body theories. Distribution of heat transfer 
coefficient over the surface of the same shapes, in both laminar 
and turbulent flow, has been obtained over the same Mach 
number range. This has permitted a selection to be made 
among the several alternative theoretical methods of cal- 
culating heating rates at re-entry conditions. 

Probably the most valuable contribution of the X-17 re- 
entry flight data has been in the experimental determination 
of the criteria for transition from laminar to turbulent flow 
on a re-entering vehicle. As described earlier, a typica) X-17 
experiment covered the Reynolds number range at which 
transition is likely to occur, and the transition point could be 
quite accurately located from the measured distribution of 
surface temperature. The position of the transition point is of 
vital importance in the determination of re-entry heating on 
ballistic missiles, and much of the X-17 test program was 
therefore aimed at investigation of factors influencing transi- 
tion. In addition to variations in nose cone shape, the tests 
included several alternative surface roughness conditions, 
with both distributed and local roughness elements. 
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The maximum Mach number attained in these flights, about 
15, is adequate for direct application to intermediate range 
ballistic missiles, but is less than the Mach number at which 
high heating rates may be expected when re-entering from a 
satellite orbit or long range ballistic missile trajectory. In 
spite of this lower Mach number, heating rates encountered 
during the X-17 experiments were well in excess of the heating 
rates to be expected on these higher speed re-entering vehicles. 
In consequence, the X-17 served as a useful test vehicle for 
the nose cone materials to be used to withstand re-entry 
heating conditions. On several flights, the measured surface 
temperatures in critical areas reached the melting point of the 
material, but the nose cone generally maintained its struc- 
tural integrity under these extreme conditions and gave early 
assurance that the re-entry problem could be solved. 

To appreciate how the X-17 flights could produce heating 
rates exceeding ICBM or satellite re-entry conditions, al- 
though its Mach number was limited to a maximum of 15, it 
should be noted that the full-scale re-entry vehicle is gen- 
erally several times larger than the model carried on the X-17. 
Since the special trajectory and firing program of the X-17 
were designed to give full-scale Reynolds numbers at the 
critical high heat rate conditions, these occurred at lower alti- 
tude and at several times the atmospheric density appropriate 
to full-scale re-entry. The rate of convective heat transfe 
H to the surface is defined very generally by an expression 


H = kyCyppVAT per unit area 


where ky depends primarily upon Reynolds number and is 
therefore roughly the same on the X-17 model and the full- 
scale vehicle, and the density p is several times larger on the 
X-17 model; Cp, V, AT, the specific heat, velocity and air-to- 
surface temperature difference, are dependent upon Mach 
number, the product varying roughly as (Mach number)*. 

In general, the Mach number of the X-17 test was suffi- 
ciently close to the corresponding full-scale value so that 
the influence of the small model scale predominated in this 
expression; the heating rate on the model, per unit area, was 
usually greater than the corresponding full-scale condition. 

The X-17 re-entry test program has made an interesting 
contribution to the problem of radio transmission through 
the ionized boundary layer on a hypersonic vehicle. At the 
highest speeds, around Mach number 15, the transmitted tele- 
metering signal from it invariably showed a dip in signal 
strength which could be correlated quite well with the 
actual speed from flight to flight. Evidently this resulted 
from the commencement of ionization of the air near the 
vehicle, and some work was done to identify the mechanism 
of signal loss from ionization, whether from antenna mis- 
match, antenna pattern change or direct absorption. 

The most recent applications of the X-17, as a high altitude 
test vehicle in almost vertical flight, have contributed con- 
siderably to the study of the influence of the Earth’s magnetic 
field on ionized gas flow at extreme altitudes. The general 
nature of these effects has recently been discussed in a number 
of articles in the technical press, and will therefore not be re- 
viewed in this paper. 


Appendix Stability and Dispersion Problems 


One of the major factors contributing to the reliability of 
the X-17 vehicle is the absence of guidance equipment. 
Operation without guidance implies that the various errors 
which can modify the flight path—the errors in launch, the 
asymmetries of the missile itself, the misalignments in the 
thrust and the atmospheric disturbances—remain uncorrected 
during the flight and contribute to the final dispersion. To 
insure that the changes in trajectory due to these factors do 
not prevent attainment of the objectives of a particular 
flight, it has been necessary to undertake a comprehensive 
program of calculation of the actual flight paths, taking into 
account probable values of the major disturbances. Since it 
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is necessary to include the changes in mass and moment of 
inertia of the vehicle as the fuel is expended, the aeroelastic 
distortions of the vehicle structure under the aerodynamic 
loads and the effects of roll on the missile dynamics, the ca!- 
culation of the flight trajectory is fairly complex. The equ:- 
tions of motion including all these variables have been solvi«| 
directly by the use of high speed digital computers; in addi- 
tion, a method of separation of the overall flight trajectory 
from the local oscillation at a particular point has been utiliz 
to permit graphical solution of the oscillatory motion on :n 
analog computer. 

In the present paper, it is not feasible to present the resu!'s 
of these calculations in detail. Instead, the following par.- 
graphs will outline the most important general features whi: 
appear from the results. Evidently the critical question s 
the response of the vehicle to any disturbance; i.e., how quick 
the disturbance will be damped out as a result of the aerod: - 
namic damping. There are damping components due to t! e 
acceleration of the missile under the motor thrust, and al-» 
damping from the change in angular momentum with time : s 
the rocket grain burns; these are added to the normal aerod: - 
namic damping, which arises mainly from the large fins of t! e 
first stage. It appears from the results of our calculations, 
however, that the major effect contributing to the dampii.z 
is the influence of varying air density, which changes tlc 
natural frequency of pitch oscillation as the vehicle tran-- 
verses the atmosphere. It has been shown by Friedrich an | 
Dore (1)? in a recent paper that the damping of a vehic'e 
passing through a region of varying air density will depend 
upon the change of aerodynamic stability in accordance wit! 
the approximate equation 


2 
1 


Ay Wn2 
where 

A = amplitude of the oscillation 

wn = frequency of pitch oscillation 

¢ = damping factor 
CnagS T-D ly 

mV mV Ty 


using conventional aerodynamic symbols. 

The natural frequency of pitch oscillation w, varies roughly 
in accordance with the dynamic pressure. At the beginning 
of the trajectory it increases steadily as speed is gained, up to 
an altitude of about 35,000 ft. It then drops rapidly due to 
the increase in air density. In accordance with the equation 
given above, the oscillations at launch and disturbances on 
the initial part of the trajectory are very rapidly damped, but 
above a height of 35,000 ft the damping ceases to be as effec- 
tive, until at extremely high altitudes the motion becomes 
completely divergent. This is clearly illustrated in Figs. 13 
and 14, which show the oscillations of a typical] trajectory of 
the X-17 utilized in a re-entry test. It will be seen that the 
large angle of attack induced at launch dies out within about 
7 sec, and that the subsequent variations in angle of attack 
are extremely small until a period of about 40 sec has elapsed. 
The oscillations then begin to build up to a complete diver- 
gence at the maximum altitude. On the return path to 
ground, the rapid increase in dynamic pressure again induces 
an extremely large effective damping, so that the large error 
in angle of attack which has been built up at high altitude is 
very rapidly reduced to zero. 

The experimental data collected on flights of the X-17 ve- 
hicle to the present have generally confirmed these calculations 
on the nature of the disturbances in trajectory. Examining 
first the ascending part of the trajectory, it is evident from 
Fig. 13 that the disttrrbances occurring in the few seconds im- 
mediately after launch are large in magnitude, and althoug!: 
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they are heavily damped, they can be expected to contribute 
considerably to the final dispersion of the vehicle. Experi- 
mentally it has been found that the greatest single cause of 
dispersion is the effect of prevailing winds on the vehicle dur- 
ing the initial part of the flight, and that this effect is felt for 
approximately 6 sec of flight time (to an altitude of about 
1000 ft). For this reason, the prevailing winds are included 
in predictions of the flight path, and a final prediction is made 
on the basis of the wind determination approximately 15 min 
in advance of the flight. It is of interest to compare this pre- 
diction with the flight path deduced from the wind strength 3 
min after the flight. As a typical example, the dispersion in 
a horizontal plane at an altitude of 200,000 ft has been deter- 
mined, both from the calculations on the basis of these two 
wind determinations and also from the actual flight trajectory 
as determined by radar. Calculations based upon the wind 
strength 15 min before flight gave a dispersion at 200,000-ft 
altitude which averaged 17,400 ft from the actual determined 
position of the vehicle. When the wind determination 3 min 
after the flight was utilized, the mean error in the same 9 flights 
w:s only 13,400 ft. 

n operation as a re-entry test vehicle, there is no further 
m:. jor disturbance on the ascending trajectory, and separa- 
tion of the first stage followed by ignition of the second-stage 
ro-kets occurs in the damped part of the descending trajec- 
tory. Fairly large disturbances are therefore permissible 
without significant effect upon the resulting motion. Calcu- 
lations show that a disturbance inducing an angle of attack of 
20 deg in the second-stage rocket motor does not -produce 
critical conditions. This is fortunate, since the second stage 
consists of three motors in parallel, and there is a possibility of 
delay of ignition of one of the motors relative to the other two. 
Extension of the calculations to include this delay shows that a 
delay of 0.02 sec does not produce seriously high angle of at- 
tack or critical loads. However, if one rocket should fail to 
ignite completely, the effect would be disastrous, since the 
moment produced is equivalent to a thrust misalignment of 
about 3 deg. There is no case yet during the flights of the 
X-17 where such an accident has occurred. 

When the X-17 is used as a high altitude sounding rocket, 
conditions at separation of the first and second stage are 
much more critical than in the re-entry mode of operation, 
since these separation disturbances and thrust misalignments 
then occur during the phase of flight where the decreasing dy- 
namie pressure is effectively reducing the vehicle damping. 
It is then necessary to resort to spinning of the vehicle to 
minimize the influence of misalignments. Calculations show 
that with a spin rate of 2.5 eps the second stage will develop 
a maximum angle of attack of approximately 1.9 deg, with an 
associated maximum change in flight path angle of 1.15 deg, 
if a thrust misalignment angle of 0.2 deg is assumed. The 
problem of the third stage during high altitude sounding rocket 
operation is even worse; the same thrust misalignment of 0.2 
deg will produce a maximum angle of attack of 20 deg with the 
same spin rate of 2.5 eps and a corresponding change in flight 
path angle of about 6 deg at third-stage burnout. To reduce 
these large disturbances, it is necessary to resort to much 
higher rates of spin, approaching 10 rps. Fig. 15 shows the 
effect of the third-stage spin rate upon the maximum angle of 
attack at burnout, and illustrates the very powerful effects of 
high spin rate. 

It is appropriate at this point to discuss generally the in- 
fluence of spin upon the flight path as utilized in the X-17 
vehicle. It is important to emphasize that in comparison 
with the spinning shell, the rates of spin are low, so that 
gyroscopic torques and Magnus effects are usually negligible. 
At extremely high altitudes the gyrostopic torques may be 
significant, when the aerodynamic moments are negligible 
and the rate of change of angle of attack is low, as is usually 
the case. The gyroscopic effects then convert a pure diver- 
gence into a precision motion, so that the deviations can be 
l-ss than in nonspinning flights as shown in Fig. 14. At lower 
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Fig. 13 Angle of attack oscillation in typical re-entry trajectory, 
ascent 
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Fig. 14 Angle of attack oscillation in typical re-entry trajectory 
(continued), descending phase 


altitudes, the spin is merely utilized to ensure that asym- 
metries in missile construction and in motor thrust do not 
build up a large change of angle of attack during the motion. 
Disturbances which do not rotate with the missile, such as 
wind shear or gusts, will have almost identical effects whether 
the missile is spinning or not. However, at high altitude, the 
wind effects are evidently small, so that the moments created 
by the missile itself due to fin or thrust misalignments are 
much more significant; in this case, spinning can be ex- 
tremely valuable in minimizing the resulting disturbances, as 
is evident from Fig. 15. To illustrate its major effects, it is 
instructive to examine the trim characteristics of a missile as 
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the rate of roll is increased from zero. Fig. 16 shows both the 
magnitude of the resultant trim and the phase angle between 
the roll and no-roll conditions, as a function of the ratio of the 
roll frequency to the pitch frequency, the individual curves on 
each figure referring to changes in the damping factor as de- 
fined above. It will be seen from this figure that to produce a 
favorable effect upon the motion, it is necessary that the spin 
rate should be much more rapid than the pitch frequency. 
When the ratio of the rolling frequency to the pitch fre- 
quency approaches unity, the effect known as “‘roll resonance’ 
will occur. The disturbances then build up continuously wit |: 
time, so that extremely high effective trims can be obtaine:| 
if the damping is not high. At high altitudes, as discussed, 
the damping is usually extremely low, so that coincidence of 
rolling frequency and pitch frequency can be catastrophic in 
its effect upon the flight path. To avoid these effects, an! 
also to take maximum advantage of the influence of high spi: 
rates in reducing the effect of disturbances, the third stage o/ 
the X-17 is provided with jet tabs which spin it to the hig): 
rate of about 10 rps for high altitude sounding work. Th» 
missile passes through the resonance condition in reaching this 
spin rate, but the acceleration through resonance is sui- 
ficiently rapid to prohibit a significant buildup in oscillation. 
The spinning of a large missile at such high rates presents 
some structural problems, and, indeed, the upper limit to 
which the missile can be spun is determined entirely by the 
structural characteristics. The structural problem involve: 
is essentially the classical problem of shaft whirl, but there are 
a number of complicating features. The center of mass is not 
exactly along the axis of the missile, but has a distributed ec- 
centricity. In addition, the spin rate is time varying, so that 
the problem of the whirling shaft must be generalized. A 
further effect which must be taken into account arises from 
the thrust and drag applied to the vehicle, making the prob- 
lem analogous to that of a whirling shaft in compression (2). 
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Pressure, temperature, chemical composition and electrical conductivity data are presented for 
gas mixtures behind strong air shocks and detonations of hydrogen-air and acetylene-air mixtures 
up to Mach 20. A digital computer program has been developed which solves the simultaneous 
aerodynamic and chemical equilibrium equations. The latest thermodynamic and collision cross- 
section data are used. Data are also presented wherein the initial gas composition was varied to 

- obtain the maximum electrical conductivity for a given shock or detonation strength. 


»CIENTIFIC and engineering interest have been focused 
J upon the properties of high temperature gas mixtures by 
the increasing effort being devoted to high speed gas dy- 
namies. Dissociation and ionization affect the properties 
of the gas mixture and open up new opportunities for study 
and application to the chemist, physicist and engineer. For 
example, magnetohydrodynamics is critically dependent 
upon the concentration of ions present in a gas mixture. 

One-dimensional steady flow equations of fluid mechanics, 
together with the concept of thermodynamic equilibrium and 
the basic thermochemical data are the foundations upon which 
such studies rest. The basic equations are presented in the 
next section, and the sources of thermochemical data which we 
used are listed in the Appendix. The problem is reduced to 
the simultaneous solution of a large number of nonlinear alge- 
braic equations. In our particular case, where we consider 
chemical systems with hydrogen, oxygen, nitrogen, carbon 
and argon atoms, the number of equations is 15. A digital 
electronic computer (IBM 650) was used to solve this system 
of equations. 

Knowing the thermodynamic state of a gas mixture, i.e., 
pressure, temperature and chemical composition, permits the 
calculation of the electrical conductivity of the mixture if the 
electron-atom collision cross sections are known. Such cal- 
culation has been carried out here for several important cases 
of air shocks and strong hydrogen- and acetylene-air detona- 
tions. Finally, we have explored the effect of initial chemical 
composition to find what mixtures will yield maximum elec- 
trical conductivity in the product gases of a given strength 
shock or detonation. 


Fundamental Equations 


Consider the one-dimensional steady flow of a homogene- 
ous, chemically metastable, gaseous mixture. The thermo- 
dynamic and physical properties and composition of this re- 
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actant mixture are assumed known. We seek to determine the 
thermodynamic and physical properties and equilibrium com- 
position of the product mixture. Specifically we treat the 
class of reactions 


aoCrHy + + a2N2 + asA B,CO2 + + 
B.0: + B;CO + + + + BoH + + 
BuHe + + + BuO; + BisNO2 + BisNO* [1] 


where z and y determine the particular hydrocarbon fuel enter- 
ing the reaction, and ao, a, @, a3, Vi, T; and p, are known. 
The last three quantities are the velocity, temperature and pres- 
sure of the reactant gas, respectively. The final four species 
treated in Equation [1] are present only in relatively small 
amounts, their concentrations being about 10-4 times those 
of the major product gas constituents. As a result, these 
species are treated as a perturbation to the main problem, 
to determine the 15 quantities: = 1,2,...,12, V2, Ts and 
p2 of the product gas. From the solution, other properties of 
the product gas can be calculated. 

To determine the 15 unknowns, we require 15 relations be- 
tween them. These are provided by the various conserva- 
tion laws and the chemical equilibrium equations. 

For one-dimensional steady flow, conservation of momen- 
tum is expressed by 


Pi + = pot p2V2? [2] 
Conservation of mass is expressed by 
= prVe [3] 


The fact that the total number of atoms of each atomic 
species entering the reaction given by Equation [1] remains 
constant is expressed by 

zx Bs 
= 282+ Bs + 2Bu 


1 1 1 1 1 
a = B+ + By + + 5 Be 


1 1 
a: = Bs + + 


where we have set ao = 1. 
For one-dimensional steady flow of a homogeneous mixture 
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of gases, the conservation of energy is expressed by 
1 1 
+ = Vs" + TBihi(T2) [5] 


where h; is the specific enthalpy of the ith species. 

For the reaction expressed by Equation [1], chemical 
equilibrium of the products requires that the following set of 
equations be satisfied 


= 2 2N 
Nol 2) B:2p: 
Bx 
K T:) = H.O @2H+0 
#,0( 7'2) 
= O. = 20 
oh Bs*ps 
Kon(T2) OH @O+H 
BsBop, 
zB; 
Ku(T2) = H, = 2H 
Kyxo(T:) = NO2N+0 
GTA: ., 
Keo, T:) = 2CO. @ 2CO + O. 


In Equations [5 and 6] Ah,(7) and K;(T) are known func- 
tions of the temperature (see Appendix A). Equations [2 
through 6] then constitute 15 nonlinear algebraic equations 
in the 15 unknowns: £;, 7 = 1, 2,..., 12 and Ve, 72 and po. 
This set of simultaneous equations has been solved by means 
of the Newton-Raphson method using a digital computer 
(IBM 650). The details of the solution are given in (1).4 

The concentrations of the four rare species C, O3, NO» and 
NO? are then computed from the equilibrium equations 


Bc = coO2C+0 
T2) 
Bo, = O; = 30 
= [ @ Nz + 20; 
Bxo+ = NO = NO+ + 


§ The conductivity of the product is computed by using the 
equation proposed by Chapman and Cowling (2) for slightly 
ionized gases consisting of rigid elastic spherical molecules 


kT kT V8n 27 myme 
We introduce the collision cross section Q defined by 
Q = tr? 


and assume that the mass of the electron is negligible com- 
pared to that of the molecule. We obtain finally 


(mkT)”2 Q 
or 


yxo+ 1 mhos 


VT Q cm 


‘ Numbers in parentheses indicate References at end of paper. 


o = 5.16 X 10-” [9] 


where 7 = temperature, R, and Q = collision cross section, 

The value of the collision cross section Q to be used in 
Equation [9] is the sum of the Maxwell-averaged collision 
cross sections of the neutral molecules and the ions, each 
weighted with the number density. That is 


+ 
Q=- 


where 
Nn = number density of neutrals 
n; = number density of ions 


Qn = electron collision cross section of neutrals 
= electron collision cross section of ions 


Q; 
Dn j = total number density of particles 
j 


In the temperature range considered here, the number density 
of ions is much less than that of the neutrals. Therefore, the 
contribution of the ions to the total collision cross section 
may be neglected. Thus 
n 


af Din; 
Jj 


For the gas mixtures considered here, the major product gas 
constituents are Ne, N, O2,O, CO and H. Therefore, Equa- 
tion [10] may be written 


nN NOs no - nco nH 
+ = On + Qo. + = Q + — Qo + 
wah an ah 


or 
Q = + + yo.Qo. + yoQo + ycoQco + [11] 


The values of the collision cross sections for atomic oxygen 
and nitrogen are not known with certainty. Reference (3) 
presents values of this parameter for O and N. The cross 
section for H in this temperature range is unknown. A value of 
1.5 X 107 em? has been assumed, however, since this is of 
the same order of magnitude as the cross section for H at 
higher temperatures (4). In addition, cross-section values for 
the atomic species O and N are of this magnitude in the tem- 
perature range of interest. Fig. 1 shows the electron collision 
cross sections for several interesting species. The values for 
the molecular species were obtained experimentally by the 
Townsend method and are taken from (4). 


Results 


The basic equations were solved for the case of air shocks. 
The results of the computations are shown in Figs. 2, 3 and 4. 
The temperature vs. the shock Mach number curve is shown 
for three different initial pressures: 10~*, 10-* and 1 atm. 
The lower temperatures obtained at the lower initial pres- 
sures result from the increase in dissociation as the pressure 
is decreased. In addition, the waviness in the curves results 
from the various dissociation processes which are occurring. 
The dissociation which occurs in the product gas as the Mach 
number is increased is shown in Fig. 3. Of particular interest 
is the appearance of ionized nitric oxide in appreciable quan- 
tities at the higher Mach numbers. The product gas pres- 
sure as a function of Mach number is shown in Fig. 4 for an 
initial pressure of 10~3 atm. 

Figs. 5 and 6 show the results of solving the basic equations 
for the fuels hydrogen and acetylene, respectively, burning in 
air. On the curve of temperature vs. Mach number is indi- 
cated the subsonic and supersonic Chapman-Jouguet points. 
The steady-state system of equations is not solvable in the 
Mach number range between these two points. The curve of 
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Fig. 1 Collision cross sections for slow electrons 
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Fig. 2. Effect of initial pressure on the product gas temperature 
behind air shocks 


pressure vs. Mach number for both fuels was found to be 
virtually coincident with the same curve for air shown in Fig. 
4. The effect on the dissociation of the combustion products 
of acetylene of increasing the initial Mach number (hence the 
final temperature) is shown in Fig. 6. Again, the rise in the 
quantity of ionized nitric oxide at the higher Mach numbers 
should be noted. 
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Fig. 3 Constituents of the product gas behind air shocks vs. 
shock Mach number 
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Fig. 4 Product gas pressure vs. shock Mach number behind 
air shocks 


In performing magnetohydrodynamic experiments in a 
shock tube, it is desirable to use reactant gases which will 
yield the highest quantity of ionized constituents. This is a 
consequence of Equation [9] which indicates that the elec- 
trical conductivity of a gas is directly proportional to the 
mole fraction of ionized species. In the range of tempera- 
tures considered here, only ionized nitric oxide (NO*) ap- 
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Fig. 5 Product gas temperature vs. Mach number for fuel-air 
detonations 
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Fig. 6 Constituents of the product gas for a stoichiometric 
acetylene-air detonation vs. 
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Fig. 7 Effect of nitrogen concentration on the quantity of ionized 
nitric oxide in the product gas of mixtures initially O, and N: 
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Fig. 8 Effect of nitrogen concentration on the quantity of 
ionized nitric oxide in the product gas of mixtures initially H,, 
oO. and N, 
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Fig. 9 Effect of nitrogen concentration on the electrical conduc- 
tivity of mixtures initially O. and N, 


pears in appreciable quantities due to its relatively low ioniza- 
tion potential (9.5 ev). 

Figs. 7 and 8 indicate the effect of varying the quantity of 
nitrogen in O.-N»2 and H.-O.-N2 mixtures such that the largest 
amount of NO* will be yielded for Mach numbers ranging 
up to 20. The corresponding curves of electrical conductivity 
vs. the quantity of nitrogen in the reactant gas are shown in 
Figs. 9 and 10. Only at a Mach number of about 17.5 for 
the O.-N, mixtures and a Mach number greater than 20 for 
the H»-O--N2 mixtures does the optimum composition coin- 
cide with that of air. At the lower Mach numbers, a larger 
amount of nitrogen than is present in air is required to obtain 
the maximum conductivity. 

Fig. 11 is a compilation of electrical conductivity computa- 
tions for acetylene-air and hydrogen-air detonations and air 
shocks. It is seen that in the Mach number range con- 
sidered, the electrical conductivity behind air shocks exceeds 
the conductivity of the product gases resulting from H:-air or 
C.H.-air detonations. 

In an effort to attain higher electrical conductivities than 
are possible by using air or optimizing the quantity of nitrogen, 
the technique exists of “seeding” a gas mixture with a sub- 
stance that ionizes easily. The alkali elements have low 
ionization potentials and might be used to increase the elec- 
trical conductivity. The computation of the equilibrium 
constants for three of the alkali metals (cesium, potassium 
and sodium), using the Saha equation, is shown in the 
Appendix. Fig. 11 indicates the large increase in electrical 
conductivity which may be achieved using cesium to ‘“‘seed”’ 
air. The quantities of cesium shown are 0.001 and 0.01 per 
cent by volume, that is, one mole of cesium to 10° and 104 
moles of mixture, respectively. 


Aveust 1959 


p, 
T,=537° R 
| MOLE OF 

MOLE OF Op 


| 


= 


ELECTRICAL CONDUCTIVITY (MHOS/CM) 


\ 
5 + 
| \ 
MAXIMUM 
T 
| 
mya 
/\2.5 / 
| 
M=10 
| | 
04 06 OF 08 09 10 


INITIAL MOLE FRACTION OF No 


Fig. 10 Effect of nitrogen concentration on the electrical conduc- 
tivity of mixtures initially H., O. and N, 


100 
AIR SHOCK +0.0!1% Cs 
AIR SHOCK +000! % Cs 
Fj 
= 
5 
AL 
> AIR SHOCK —7/ 
= 
4 STOICH. STOICH. 
-| AIR—F H2-AIR 
7 7 
77 7] 
f/ 
WwW 
Tf 7 
p, =10°ATM. 
T,=537°R —— 


8 i0 i2 14 6 20 
SHOCK OR DETONATION MACH NUMBER 


Fig. 11 Electrical conductivity vs. shock or detonation Mach 
number 


577 


| 
=| 
| 
| | } aa 
| 
| 
=| 
4 | 
4 
| 
| 
] 
| 
| 
i 
f 
‘ 


Appendix 
A Sources of Thermochemical Data 


The major sources of thermochemical data were tables 
compiled by the National Bureau of Standards (6, 7). In 
some cases, where NBS data were not available, sources which 
were believed to contain the most accurate, modern data were 
employed. Table 1 lists the heats of formation used in the 
computation and the various data sources. 

For purposes of computation, the composition of air at 
standard conditions was taken as 


Molecule Mole fraction 
0.780881 
0.209795 
A 0.009324 


We have compared the results obtained from our computer 
program with the data in (9 and 11) for the case of air shocks. 


Table 1 Thermochemical data 


Heat of formation, 


AHy, at zero deg K Source 

Gas Keal/mole (Ref. no.) 
O- 0 (6) 
O 58.586 (6) 
H, 0 (6) 
H 51.620 (6) 
No 0 (9) 
N 112.535 (9) 
OH 10.0 (6) 
NO 21.477 (6) 
—57.107 (6) 
CO, —93 . 9686 (6) 
CO —27.2019 (6) 
C (solid) 0 (6) 
C (gas) 170.39 (6) 

Free-energy function or 

Equilibrium Constant Enthalpy 
Gas (Ref. no.) (Ref. no.) 
N2 (7) (7) 
N (7) (7) 
(7) (7) 
O (7) (7) 
Hz (7) (7) 
H (7) (7) 
OH (10) (8) 
NO (9) (8) 
H.O (10) (7) 
CO, (7) (7) 
co (7) (7) 
Cc (9) (8) 


Table 2 Equilibrium constants for the alkali metals 


Metal Na K Cs 
Ionization poten- 
tial, ev 5.14 4.33 3.87 
Temperature, K logioK 
1000 —24.901 —20.818 —18.499 
2000 —11.193 —9.151 — 7.992 
3000 —6. 434 —5.073 —4.300 
4000 —3.963 —2.942 —2.362 
—2.425 —1.608 —1.145 
—1.363 —0.682 —0.296 
—0.578 0.004 0.336 
0.029 0.539 0.829 
0.517 0.970 1.228 


Equilibrium composition of the product gas was found to 
agree very well with both references. 

A comparison of our electrical conductivity data with the 
data in (5) for the case of air shocks indicates excellent agree- 
ment over the entire Mach number range considered. 

In (1), the heat of formation of nitrogen was taken as 85.120 
kcal/mole. The value presently accepted is 112.535 kcal/mole 
which, as indicated in Table 1, was used in the present report. 


B Computation of the Equilibrium Constants for 
Ionization Reactions Involving the Alkali Elements 


To compute the thermal ionization of a monatomic gas, 
Saha has applied the Nernst equation for the equilibrium 
constant. By neglecting the contribution of the upper elec- 
tronic states of the neutral atom to the specific heat, and con- 
sidering the atom, ionized atom and electron to be ideal gases, 
he derived the following equation for the equilibrium con- 
stant 


5 


log K(T) = — 


For the electron, g. = 2. In the case of the alkali metals, ga = 
2andga+= 1. Inaddition, U = 23,0537. Substitution into 
Equation [B-1] gives 


5 
log — 6.491 |B-2] 


log K(T) = —5041 


where J is in electron volts, and T is in deg Kelvin. 

Equation [B-2] has been used to solve for the equilibrium 
constants of the three alkali metals, sodium, potassium and 
cesium, in the temperature range from 1000 to 9000 K. Table 
2 indicates the results. 

It should be noted that these results in Table 2 assume the 
occurrence of singular ionization only, according to the reac- 
tion equation 


A = At+t+e- 


where A is the neutral atom, A+ the ionized atom, and e~ the 
electron. The equilibrium constant is defined by 


log K(T) = »| 


Ya 
where 
Yn = mole fraction of the nth constituent in the equilibrium 
mixture 
p = pressure of the mixture 
Nomenclature 
Dy» = diffusion coefficient 
g = statistical weight 
h = specific enthalpy 
I = ionization potential, ev 
k = Boltzmann constant 
K(T) = equilibrium constant 
m = molecular weight (see Eq. 5); mass 
Me = electronic mass 
Ne = number density of electrons 
nj = number density of particles 
Pp = pressure 
Q = electron-particle collision cross section 
re = distance between particle centers at collision 
= temperature 
U = ionization potential, cal per mole 
V = velocity 
x = number of carbon atoms 
7] = number of hydrogen atoms; mole fraction 
a = moles of reactants 
B = moles of products 
€ = electronic charge 
p = density 
o = electrical conductivity 
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Meteoric Ablation Through 


Interfacial Instability 


SIN-I CHENG’ 


Princeton University 
Princeton, N. J. 


The interfacial instability of the Taylor-type of the molten layer on the surface of a meteoric body 
is analyzed when the representative situation on the surface of a blunt nosed body is taken into 
consideration. The presence of a component of the effective gravity field tangent to the gas-liquid 
interface considerably increases the seriousness of the instability, especially in the region of the 
transonic shoulder. In the transonic region, the logarithmic increment per second of the most 


critical component of initial disturbance is likely to be significantly larger than that in the stagna- 
tion region where the normal component alone acts. In view of this result, the recovered meteor- 
ite, with smooth nose but ablated downstream, as illustrated in Fig. 2B, can be understood from the 
ablation mechanism through interfacial instability. The value of the logarithmic increment per 
second at a given location is inversely proportional to the local value of the kinematic viscosity of 
the liquid at the interface. The contribution from the normal component of the gravity field is 
proportional to R'? (where R is the radius of curvature of the body, a linear scale representing geo- 
metrically similar bodies), but is rather insensitive to the thickness h of the liquid film. The con- 
tribution from the tangential component is insensitive to the body seale R, but is linearly propor- 
tional to the local thickness h of the liquid film. The total logarithmic increment per second of the 
combined excitation over the blunt nose region will, therefore, increase with the linear scale of the 


body and with the thickness of the molten layer, but not so seriously as ~R!/? and ~h. 


"P.HE VELOCITY of meteoric objects entering the atmos- 

phere is usually so great that the impact of the atmos- 
pheric particles will “melt” the surface material of the body. 
The surface of the meteoric object will be covered by a layer 
of the molten material when the collisions are sufficiently fre- 
quent. The stability of the molten layer, therefore, becomes 
of interest. An energy approach (1)* has been developed to 
treat such stability problems. It is shown that the presence 
of an effective gravitational field normal to the liquid-gas inter- 
face and directed toward the gas side is important (1). In ad- 
dition, it has been shown that the interface in the stagnation 
region is always unstable, and the instability is predomin- 
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antly of the Taylor type, possibly leading to mass ablation 
from the body (2). A method of determining the amplifica- 
tion rate based upon the energy approach is also presented 
in (2) and illustrated with sample results. It is found that 
the interfacial instability can be quite vigorous, and it is. 
therefore, considered necessary to determine the amplification 
rate more carefully by taking into account several factors 
which were ignored in (1): The effect of the tangential com- 
ponenent of the gravitational field along the local liquid-gas 
interface, and the effect of the variation of the viscosity co- 
efficient of the liquid, etc. For this purpose, the energy 
equation must be modified slightly. The analysis of the insta- 
bility condition of the molten layer on different portions of the 
meteoric surface, from the stagnation point downstream, leads 
to several conclusions of practical importance. 

It is found that the tangential component of the gravity 
field is always destabilizing over the nose region, and by its 
action alone, the transonic region over the blunt nose may be 
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Fig. 1 Schematic diagram of coordinate system 


more vigorously unstable than the stagnation point region. 
The destabilizing action of the normal component of the 
gravity field increases with the linear scale of the nose di- 
mension. Simple formulas are given for estimating the 
logarithmic increment per second (L.I./sec) over this subsonic 
and transonic nose region. The conditions in the region down- 
stream of the transonic region cannot be determined without 
detailed numerical work for specific instances. 


Energy Equation for Disturbance Motion 


When a small segment of the body surface is considered, 
it will be assumed that the mean or the undisturbed flow of 
the liquid layer and the gas stream is essentially a locally 
parallel flow. The undisturbed position of the gas-liquid 
interface will be chosen as the z-axis (y = 0) with the solid- 
liquid interface located at y = —1. Following (1), the linear 
scale is here nondimensionalized by using the local thickness 
of the liquid film h(cm) as the reference length. The mean 
flow properties are necessarily the local values. Since the 
film thickness h is small compared with other linear dimen- 
sions of the body, the variations of the mean flow properties 
along the surface, within a length equal to A are small, and 
the assumption of nearly parallel flow is then justifiable. 

The object under consideration is moving at a very great 
velocity with respect to the surrounding atmosphere. The 
deceleration, large as it may be, will not change the relative 
velocity percentagewise to any appreciable extent during a 
period of, say, 1 sec. Nevertheless, the actual displacement 
of the object during 1 see may be of the order of 10 or more 
miles with consequent variation of the atmospheric environ- 
ment due to altitude change. In order that a quasi-steady 
mean flow may be assumed, the results of (1) clearly indicate 
that the two characteristic times h/U» and v,/Gh should be 
sufficiently small. Here Uo is the mass velocity of the liquid- 
gas interface, G the magnitude of the deceleration and y; the 
kinematic viscosity of the liquid in consistent dimensional 
units. Representative conditions are such that this quasi- 
steady-state assumption is well justified. 

The mean pressure p will be taken as constant across the 
liquid layer, but the mean temperature 7’ must vary signifi- 
cantly. At the gas-liquid interface, the liquid temperature 
will be rather close to the stagnation temperature of the gas; 
while at the liquid-solid interface, the liquid must be at the 


u? + v? f dU 1 f ( or f du \ du 
~ dV = - —dV—- — iV—2 — »— }— dV 
ol p 2 dy or or oy] dy 


A 


ov 


melting temperature, which is much lower than that of the 
gas-liquid interface. The absence of a sharp melting tem- 
perature for an amorphous solid is not of great significance as 
later results indicate. The large temperature gradient pro- 
duces a high heat transfer rate Q into the solid to melt further 
the solid surface. The density of the liquid may be taken as 
constant if the thermal expansion is neglected—an assump- 
tion easily justified. The large temperature variation across 
the liquid layer will, however, lead to a large variation of the 
viscosity coefficient u(7). It is not expected that the vari- 
ation of viscosity coefficient will result in any change in the 
qualitative stability behavior of the liquid layer, but it may 
equally well be asserted that any variation will change the 
amplification rate quantitatively. 

The deceleration G (cm/sec?) of the body, and also the 
effective gravity field when referred to body frame of refer- 
ence will, in general, be inclined to the normal of the loca! 
liquid-gas interface (the y-axis) at an angle 6. The compon- 
ent of the gravity field normal to the interface will then be 
G, = G cos 6, and that tangential to the interface will be 
G. = —G@sin 0. The gravitational potential may then be 
written as Q = yG, + 2G, where both G, and G, will be taken 
as positive if directed in the positive y- and positive x-direc- 
tions, respectively. The positive z-direction is defined as the 
direction of the flow velocity U» (cm/sec) of the interfacial 
mass. This velocity will be taken as the reference velocity 
so that U(y = 0) = 1 (Fig. 1). The tangential component of 
the gravity field G, tends to modify the boundary layer flow 
so that the mean velocity profile U(y) will be different accord- 
ing to whether the G, field is present or not. The tangential 
field G, will exchange energy with the disturbance motion if 
it involves net mass displacement in the direction of the field 
G,. The tangential field should, therefore, appear directly 
in the equation of energy balance of the disturbance motion. 

Let us consider only the two-dimensional disturbances, 
propagating in the direction of the mean flow (z-direction) 
defined by the disturbance stream function 


U(r, y, = gly) exp fia(x — ct)}} 
with 
u = vy = exp [ia(x — ct))}} 
v = = —iag(y) exp [ia(x — ct)]}} 
c= + 1¢; 


and the instantaneous position of the disturbed interface 


yo = exp — ct)}} [2] 
where 
G = real part of the terms following 
@ = angular wave number = 27 times the number of com- 
plete waves in a distance equal to h in the z-direction 
Cr = wave velocity in the units of U» 
ac; = amplification rate in the unit time h/U) 


The stream function is used because the density of the media 
(liquid as well as gas in the vicinity of the liquid-gas inter- 
face) isassumed to be constant. Take the scalar product of the 
disturbance velocity vector (u, v) with the equation of motion 
and integrate over the control volume V defined by —1 < y < 
o and0<2< A = 2x/a, ie., a strip bounded by the solid 
and open on the gas end, and extending over a wave length 
in the z-direction, asin (1). The following expression for the 
rate of change of the disturbance energy of all the fluid ele- 
ments occupying the control volume V at any instant t is ob- 
tained after neglecting terms which are of third order or 
higher in disturbance quantities 


[3] 
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where A indicates the discontinuity of the following quantity 
across the liquid-gas interface, i.e., the difference of the quan- 
tity evaluated at y = O- lessthat aty = 04. Fy = Uo?/G,h, 
F, = U,?/G.h represent the component Froude numbers 
with proper sign associated with G, and G,. The first term 
on the right-hand side is the Reynolds stress excitation. The 
following group of terms is the direct viscous contribution, 
which vanishes as the Reynolds number R, = p,Uoh/w ap- 
proaches infinity. Here yu is nondimensional, expressed in 
terms of the liquid viscosity coefficient uo just below the 
liquid-gas interface, i.e., u(y = O-) = 1. The first integral 
in this group is the ordinary viscous dissipation term, pro- 
portional to the square of the vorticity. The second integral 
is due to the variation of the viscosity coefficient in each 
layer, and the third term is the contribution from interfacial 
discontinuities. The last integral over the wave length 
is due entirely to the interfacial discontinuities of fluid density 
and of static pressure in the presence of surface tension. 
Only two terms in Equation [3] involving du/dy and F, are 
new when compared with the corresponding equation in 

Following exactly the procedure as described in (1), replace 
uv and v from Equation [1], yo from Equation [2] and Ap by 
(1/W)(0%yo/Ox?) with W designating the Weber number, 
I" = p,Uo*h/r (7 stands for surface tension, dynes/cm), and 
carry out the integration with respect to x in Equation [3]. 
Then divide the whole expression by 


Nono*c; exp = gogo* exp — + 
a 


The last equality is due to the assumption that the liquid-gas 
interface is a streamline at all times. Here (*) indicates the 
complex conjugate of the quantity. The final expression is 
then 


M : 
at[(1 — = A’ + [(1 + (E D) 


By applying the timewise average over the period 22/ac,, 
the time independent terms M, E and D remain unchanged 
as given by Equation [2], but the interfacial excitation A’ be- 
comes 


1 a? 1 — (1 — cjg’ 


[7] 


and the disturbance energy balance over a wave length A = 
2r/a and over a period 27/ac, is 


— + = At (1 6)? + (E 


This equation is identical in form to the corresponding equa- 
tion in (1), but the definitions of A and D are slightly dif- 
ferent. The presence of the tangential component gravity 
field G, or 1/F, in dimensionless form greatly increases the 
complexities of the problem, even qualitatively. This is due 
primarily to the fact that A the interfacial excitation term 
can no longer be reduced to a form independent of the dis- 
turbance motion. This dependence is extremely important. 
The contribution to the disturbance energy of the tangential 
component gravity field may add to, or subtract from, that 
of the normal component, depending upon the mean flow 
condition. The simultaneous presence of a tangential com- 
ponent of the gravity field and a normal component of the 
gravity field directed toward the less dense medium, will 
not necessarily result in the Taylor-type instability of the 
interface. Moreover, it is no longer clear that the waves, 
nearly stationary with respect to the mean interfacial mass 
motion, are the most unstable. Although the normal gravity 
field does not participate in determining the stability bound- 
ary of the waves not stationary with respect to the mean mass 
motion at the interface (c, ¥ 1), the tangential component 
does. That the contribution of the tangential component 
of the gravity field may or may not add to that of the normal 


[4} 
; component, depending on the fluid mechanical situation, 
where can be visualized directly from the following discussion. 
as The normal component of the gravity field is operative in 
f 2 + energy exchange with the y-component mass velocity of os- 
2 = a Qyyr* cillations, whereas the tangential component is operative 
with the z-component. The velocity oscillation of the inter- 
B= 
a gy dy 
2 i - 4 fe du | 


2F, 


Comparison with the corresponding results in (1) shows that 
both the apparent mass M and the Reynolds stress excitation 
remain unchanged in form. The viscous dissipation term 
D is augmented by a new term proportional to du/dy. The 
interfacial excitation term given as Equation [6] is greatly 
complicated especially because part of the contributions of 
G, ~ 1/F, are time dependent. The situation may be simpli- 
fied if one agrees to consider the energy balance not at any 
instant ¢ but during a period of oscillation of the disturbance 
2mr/ac, (in fact, half a period would suffice). This will be 
conceptually acceptable if 27/ac, is not so large that the ac- 
tual time 27h/ac,U» would still be small enough to admit the 
quasi-steady mean flow. 
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11 
— R[i(1 — sin (2ac,t) + F. — [sin? (act) — 3] 


[6] 


facial mass will clearly depend on the many variables of fluid 
mechanical nature, of which one is the mean shear flow, that 
govern the disturbance motion of the flow field. It is, there- 
fore, imperative to have an estimate of (y’/¢),-0 before 
the contribution of the tangential component gravity field 
may be discussed any further. ‘ 


Effect of the Tangential Component 
of the Gravity Field 


The linearized equation of motion for small disturbances 
in an incompressible medium but with varying viscosity co- 
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efficient is slightly different from the Orr-Sommerfeld equa- 
tion. Two terms'associated with du/dy and d*u/dy* need be 
added, but they do not appear to complicate the analytic 
nature of the disturbance equation which is obtained by di- 
rect substitution of Equation [1] into the linearized equation 
of motion after the elimination of the pressure disturbances 


2 du 1 [ 


— ary) — U"(y)¢] [9b] 


for the liquid and the gas streams, respectively. Here R.. = 
piUoh/ uw with designating the liquid viscosity coefficient 
at the interface. u(y) is the liquid viscosity coefficient in 
terms of uso that u(y = 0-) = 1. The gas Reynolds num- 
ber is defined as R,, = p,Uch/ug with py and yp taken to be 
constant throughout. The inviscid approximation to this 
equation, obtained by dropping terms of 0(R.~') compared 
with terms of 0(1) naturally remains the same as the one with 
constant viscosity 


[U(y) — e] (¢" — — U"(y)e = 0 [10] 


This is the equation that governs the bulk of the gaseous flow 
(y> 0), when R-,, is sufficiently large. For the viscous cor- 
rection near the critical layer, where U(y.) = c,, the substitu- 
tion into Equations [9] of the expansion of 


with 
n=(y—yod/e  € = (aR. [12] 


leads to the following equation for the viscous correction F(7) 
a’ a? 
— => 1 
‘a F() + age 1 [13] 


This equation is again identical with the one for constant 
viscosity. Accordingly, the argument, identical with the 
one given in (1) based upon well-known mathematical re- 
sults (3) leads to the following condition at the critical layer 


(- ( [; In a| 


i f 
ves [U(y) — + 


[14a] 


or in the limit of c;—~> 0 


-(_U*"\}! “yah 


where it is understood that the principal value of the natural 
logarithm will be taken. These relations (Eqs. [14]) will be 
valid approximations so long as the critical layer does not lie 
in the interior of the liquid layer, and the value of aR,,U’. is 
sufficiently large. 

For the “stationary disturbances,” resulting in Taylor-type 
instability, the critical viscous layer overlaps, at least in part, 
the boundary viscous layer so that 


(¢'/¢)y=0 = (¢'/P)y=ye [15] 
Substitution of Equations [14 and 15] into Equation [7] gives 
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U(y) (y) e| (9a 


Under these circumstances, whether the tangential gravit, 
field is exciting or damping the disturbance motion depend: 
only on the product (—1/F.)-(—U"o+/U’o+) 2 0. Thus 
if F, < 0, as it is on the surface of a decelerating body, the 
tangential gravity field tends to excite interfacial instability 
if (—U"o./U'o4) > 0. The simplicity of this relation is de- 
ceiving however, because the component gravity field tangent 
to the interfacial flow partly determines U"o,/U'o+ of the 
mean flow. 

It may be shown that the instability of the propagating dis- 
turbances (c, # 1) remains of the Helmholtz type, and can 
be completely suppressed by a sufficiently large increase of the 
viscosity of the medium despite the excitation of the tangen- 
tial gravity field. The Taylor-type instability of the sta- 
tionary disturbances which may be excited simultaneously 
by the normal and the tangential components of the gravity 
field cannot, on the contrary, be completely eliminated, no 
matter how large the medium viscosity may be. The in- 
stability of the stationary disturbances is of more serious con- 
cern to the engineer and will now be discussed. 


Logarithmic Increment 


In order to discuss the stability behavior of the liquid film 
on the surface of a decelerating object in any general terms 
with the hope of extracting some important qualitative in- 
formation, it will be necessary to look into the method of 
determination of the amplification rate ac; from Equations 
[8, 7 and 5], with all the steady-state properties known. By 
considering representative configurations, we hope to simplify 
the expression of the key quantity of practical interest. 
From this, qualitative information may be obtained without 
extensive numerical work. As concluded in a previous sec- 
tion, only the stationary disturbances (¢, = 1), will be con- 
sidered, for which the energy balance Equation [8] becomes 


(ac;)? “ = A+ ac(E — D) [17] 


with A given by Equation [16] and Z and D by Equations 
[5]. It is assumed that A> 0, so that the Taylor-type in- 
stability will prevail for the cases being considered. ac; is 
the amplification rate expressed in terms of the characteristic 
time h/U», which is very small, 0(10~* sec). Therefore, 
ac; will remain numerically small even for vigorously un- 
stable disturbances with the logarithmic increment of the 
order of 10 to 100 per see. As indicated by the sample calcu- 
lations reported in (1) and those for other representative con- 
ditions, the inertia terms (ac;)? M/2 and the Reynolds stress 
terms E may be ignored as small quantities compared with 
D, at least at small R.., primarily because of the small h and 
large v, of interest. Thus the simplified formula 


A ARa 
Di +Dr+ Da/a? 


[18] 


ac; => 


may serve as a convenient basis for qualitative discussions 
concerning the effect of individual parameters if R,., remains 
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0(1) orsmaller. Here 


D; = A1 — (- In (AR 
Dif f » (¢ 
Qa YoYo 
Oo du | ¢’¢* 


are the viscous dissipation contributions: 


); due to interfacial discontinuity 

qa due to distributed vorticity 

, due to variation of the liquid viscosity coefficient under 
nonuniform thermal field 


Ti.ec contributions of the gaseous layer to D, and D7 have been 
ncylected on account of its relatively small viscosity. 

{f the viscosity coefficient wu of the liquid is taken to be a 
function of temperature only, and if the temperature field 
across the thin liquid layer may be evaluated by considering 
only the conduction process, D7 may be written as 


Dy; =2 = 2B.Nu 20 
r = (4 BoN u [20] 
where 
8) = representative average value of the ratio of du/dT 
to the heat conductivity coefficient 
Nuo = Nusselt number corresponding to the heat transfer 


rate @Q» across the interface under the fictitious 
temperature gradient 


li both the viscosity coefficient 4 and the heat conductivity 
coefficient \ are assumed to vary ~exp [—o07'/To] where 
7T/T> represents the temperature, then 8 in Equation [20] 
would precisely be the same @ in this exponent. For prac- 
tical molten materials, this is unlikely. Nevertheless, the 
upper bound of (du/dT)/AX will be of the order of unity. It is 
thus apparent that, for most of the molten materials under 
consideration, D7; (being of the order of Nu) will be con- 
siderably less than unity, and may be neglected when com- 
pared with D; + Da/a?. 

The relative magnitudes of D; and D,/a? depend to some 
extent on the value of the wave number a. Since the mag- 
nitudes of R,,U'’9, for practical systems are very large, the 
logarithm of aR,,U'’o4 remains large compared with unity 
even with a as small as 107!, 10-2. On the other hand, Dz is 
of the order of 107! and is rather insensitive to the variation of 
a, as shown by the results of sample computations. The vari- 
ation of viscosity in the liquid layer results in some modifi- 
cation of the disturbance Equation [9A] of g(y). This dis- 
turbance amplitude function ¢g(y) may, nevertheless, be de- 
termined in essentially the same manner as in (1). The basic 
behavior of g(y) and Dz will be expected to remain un- 
changed. Therefore, for a of 0(1) and when (—U"o4/U‘o+) 
remains positive and of O0(1), D; will certainly predomi- 
nate. But for small a, say 0(107') or smaller, Dz/a? 
will predominate, and the amplification rate ac; will vanish 
as a’. Thus, for the practical problem of obtaining an esti- 
mate of the maximum amplification rate of the most critical 
disturbance, the region of small @ is relatively unimportant. 
The amplification rates of disturbances of different wave 
numbers @ may be calculated numerically in the same man- 
nerasin (1). It appears, nevertheless, that useful qualitative 
information about the effect of the tangential gravity field and 
the variation of the stability behavior along the surface of 
meteorlike bodies may be obtained by ignoring D./a? com- 
pared with D;. This approximation will, in fact, give quite 
satisfactory quantitative results if W/F is fairly large and 
(—U"u /U is positive and not too small. The situation 
with U"o,/Uo+ > 0 requires much more careful investigation. 
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The amplification rate of the most critical disturbances will 
then be estimated after equating the gas density to zero by 
1 cos 86 — /W sin 
10 Ra (- 2 

3 eg (+ 
where @ is the inclination of the gravity field G (and there- 
fore F) to the outward normal of the local surface. In this 
form, the effect of the tangential gravitational field as repre- 
sented by the z-component Froude number F, = —F sin 0 
becomes deceptively simple. It merely serves to add a con- 
stant shift of (1/10)T = F(ac;)max/Re, of the amount (1/2) 
sin 6 to the contribution of the normal component. In terms 
of the dimensionless group T = (10F/R.:) (a@ci) max as used in 
(1), the logarithmic increment per second of initial disturbances 
may be written as 

L.I./see = (Gh/10v) (1, cos 8 + 5 sin 8) [22] 


where [', may be read from curves similar to those given in 
(1). The simplicity of the results of Equations [21 and 22] is 
deceiving, because the presence of a tangential gravity field will 
not only introduce the second term 5 sin 4 in Equation [22], but 
will also modify the values of the interfacial mass velocity 
U> of the mean flow and of U"y,./U 9+ of the mean flow veloc- 
ity profile. The second quantity U"o,/U’o; alone is im- 
portant in determining the logarithmic increment per second 
of the initial disturbances. For very thin liquid layers of large 
viscosity, a fair estimate of —U"9,/U’o; may be made from 
the local streamwise pressure gradient Op/dx and the body 
force field from 

on which I, critically depends. Thus, the effect of the pres- 
ence of a tangential gravity field is far from being merely 
additive to the result when the tangential field is absent. 
The simplicity exists only in the dimensionless parameters 
chosen in the present investigation. 

For the practical application of a decelerating meteor- 
like body, the deceleration G is related to the pressure on the 
body nose. The pressure gradient of the nose is overwhelm- 
ingly large, so that the associated deceleration has little conse- 
quence in altering the steady-state value of —U"o,/U’o+ on 
the gas side. It is because of this fortunate circumstance 
that many important qualitative conclusions may be obtained. 


Effect of Scale and Position on the Surface 
of a Meteoric Body 


The logarithmic increment per second of the most critical 
disturbance corresponding to the local steady-state conditions 
will be considered as a measure of the seriousness of the in- 
stability of the molten layer over that part of the surface of a 
meteoric body. Accordingly, Equation [22] will serve as the 
foundation of the ensuing discussion. G is the deceleration 
of the body, a constant for all parts of the surface. The vari- 
ation of the inclination 6 of the equivalent gravitational field 
to the local normal is accounted for by the expression I. = 
cos 5sin will be taken as inversely proportional 
to (—U"o./U 0+) from the mean velocity profile on the gas 
side of the interface. The quantities of primary importance 
that vary along the surface of the body are (—U"o./U‘o+), 
hand y,. It will be assumed that the body under consider- 
ation is blunt with a sphere-like nose, decelerating along the 
axis of symmetry. The front stagnation point corresponds to 
@6=0. Anaccurate determination of any of these three quan- 
tities is an extremely complicated matter. It is hoped that 
the following simple estimate will yield the correct qualitative 
results suitable for our purposes. 

The variation of (U"o+/U c+) as a function of @ will be con- 
sidered first. Its effect is pronounced only in the evaluation 
of Ty ~ /U 0+) 73 ie., it modifies the contribution of 
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the normal component of the gravity field only. It will be 
assumed that the liquid film is very thin compared with the 
characteristic geometrical scale of the body (as is the nose 
radius R with h/R < 1) and that the liquid is very viscous. 
The inertia force is then negligible. Thus, the ratio of 
(d?U/dy*), 0+ to (dU/dy),-0+ at the interface is roughly the 
same as that at the solid wall in the absence of the molten 
layer, and may be given by 


_ (4p/dz) — [23] 
Here, the inviscid pressure gradient along the surface may be 


expressed roughly in terms of the modified Newtonian ap- 
proximation as 


M 
_ 


OR 'p max Sin 20 


This approximation is essentially correct upstream of the sonic 
point, (sin 26 < 1). Here 


yM .*p../2 = dynamic head of the incoming stream 
Mops tink = pressure coefficient at the stagnation point 
R = characteristic radius of curvature of the nose 


The tangential component of the deceleration G, = —G sin 
may be expressed as 


Ady YM max'A 


G, = i fi = 1 
sin f(y) f(y) = [24] 
where 
A = the frontal area - 
pP, = average density of the body 
L-~ = effective length of the body defined by 


prA-L= the total mass of the decelerating body 
f(y) = 0(1) isa shape factor of body nose 


The shear stress S at the liquid-gas interface will be replaced 
by the skin friction coefficient of the gas stream over the solid 
surface in the absence of the liquid layer as 


S =(yM 


Thus, when all the quantities are reduced to the present 
dimensionless forms 


_ Ute WET op sing 
Ue C, [ sin 26 — sin 0 (2 i) ] [25] 


Several aspects may be noted from this expression. Firstly, 
p/p, the ratio of the average density of the body to that of 
the surrounding atmosphere is very large ~0(10*). The 
ratio L/R may be expected to be no less than 0(1). There- 
fore, throughout the region of validity of the Newtonian ap- 
proximation around the nose (sin 26 < 1), which may be 
roughly interpreted as the entire subsonic region, the ratio 
—U"o./U'o+ is essentially determined by the pressure gra- 
dient 


IR 


h/R 

which is always positive. Both h/R and C» are of the order 
of 10~* sufficiently far away from the stagnation region. 
Thus, the value of —U"o,/U"o, may be expected to be of 
order unity. When the stagnation point is approached, Cp 
also tends to vanish as 0, with —U"o,/U'o: approaching 
a finite positive limit. The limiting value at the stagnation 
point may be estimated independently as 


” 2 1/3 
Jom Vo _Jstag 


4(y — 1) WU. 


in the present dimensionless form, where 
WU. 


Vg Rv, 


isthe Reynolds number characterizing the stagnation point gis 
flow, and U,, is the tangential velocity gradient of the invis- 
cid axisymmetric stagnation flow field (compressibility and 
dissociation ignored in this stagnation point analysis). 

The value of Uz stag has been estimated from the Newto- 
nian pressure distribution over the blunt nose, and the normal 
shock relation over the stagnation region of the nose as 


om +1 pe R L(y + 0? 


where y is the average specific heat ratio of the gas over the 
nose region. If we take the following representative value;: 


h =10-cm 
U. ~ 10 km/sec 
Rw~1m 


v, of the hot air over the nose ~1 cm?2/sec, the stagnation 
point value of —U"y,/U’o, may be as big as 5. Thus, it 
appears likely that —U"o;/U’o4 will decrease from the stag- 
nation point downstream but remain 0(1) throughout the 
subsonic nose region. This is then the region where the 
simplified results (Eq. [22]) of the previous analysis will 
apply. 

By substituting Equations [26 and 27] into [22], one ob- 
tains: 

1 For small @, i.e., in the immediate vicinity of the stag- 
nation point region 


Ea +7 [28 

sec (4.2 LU. 2 
where I, is the value of T = [10(ac;)F/Ret}]max as Shown in 
Fig. 8 of (1) when —U"o,/U"o4 is taken to be unity. y, is 
the kinematic viscosity coefficient of the gas at the interface. 


2 For @ not too small, i.e., away from the stagnation 
region 


LT. RCp.. 
— = - -— + - sin 0 29 
sec V1 nes 
where 
R= representative radius of curvature of the body nose 
Cp = the gaseous skin friction coefficient referred to the un- 


disturbed dynamic head, yM ,.*p,,/2, of the incom- 
ing stream 


The result of Equation [28] at the stagnation point is in 
sharp contrast with the result of (1), in that the stabilizing 
effect of decreasing h (the liquid film thickness) is far less 
significant in the present result [28] where it is only through 
the slight dependence of on W/F = pGh?/r, Ret = piVoh/ mi, 
etc., that A exerts its influence. This unexpected result 
originates physically from the fact that viscous dissipation 
depends very much on the shear flow condition at the inter- 
face, being nearly proportional to —U"o,/U‘o+ and, there- 
fore, to h as indicated by Equation [27].. The linear de- 
pendence of U,, on the linear scale R at the stagnation point 
over a blunt body may be recognized alternatively on a purely 
intuitive basis in that the sonic velocity will be reached 
roughly over a distance R-2/4 from the axis. By holding 
—U"o./U'o+ constant in the flat plate calculation as in 
(1), the estimate of the viscous dissipation contribution fails 
to simulate the stagnation point condition. The present 
result in Equation [28] indicates the importance of the overall 
linear scale of the body nose instead of the local film thickness 
h. Equation [28] indicates that, for geometrically similar 
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noses with corresponding physical properties 
L.I./see ~ [30] 


i.c., a large linear scale of the body implies a more serious in- 
stability condition of the interface at the stagnation point. 

Equation [29] reveals additional aspects in line with this 
scale effect away from the stagnation point region. Firstly, 
the contribution of the normal component gravity field, 
(T;/20) (RC>/sin 8), similar to the situation in the 
stagnation region, is relatively insensitive to changes in film 
thickness h (so long as it remains small), but the contribution 
increases with the body linear scale R roughly as R'/?. This is 
|, ecause C’p, which is proportional to the square root of the local 
Reynolds number, introduces a R~/? dependence. Secondly, 
the contribution of the tangential gravity field (G/v:) (h/2) x 
sin @ is linearly proportional to the local film thickness h, but 
is independent of the linear scale R of the body. 

It may, therefore, be concluded that geometrically similar 
-valing to larger sizes will, in general, aggravate the insta- 
hility conditions over the subsonic nose region. Doubling 
the diameter of the nose will possibly increase the logarithmic 
erement per second by 20 to 40 per cent. An increase from 
‘0 to 13 would mean an increase of linear increment from 104% 
to 10°, a matter of some concern. A laboratory test of a 
model of 1 cm diameter may be apparently stable with an L.I. 
per sec, say 3 or 4. But on a full-scale object of 10? cm 
diameter, the L.I./sec would be 10 to 30, if scaled according to 
R4to R'", The instability on the full-scale model is intoler- 
able, despite the apparent stability of the small-scale test in 
the laboratory. The inference from the small-scale test re- 
sults concerning the interfacial instability should, therefore, 
be very cautiously made. 

From Equation [29], the relative importance of the tan- 
gential vs. the normal field in exciting instability may also 
be visualized. The tangential field will contribute more to 
the logarithmic increment per second than the normal field if 


T, RCp 
2 
sin? > [31] 


For order-of-magnitude purposes, I, will be taken as 1 and h 
as 107! cm. Let the local flow Reynolds number be of the 
order of 108 per ft (per 30cm); let Cp be written as f-R~/? X 
10-* where f is expected to be 0(1). Then with R expressed 
in centimeters 


fRi2 2/5 1/5 
6> 0, =| x 10-*] [32] 


defines the downstream region in which the tangential field is 
more important. For Rf? = 10?, 10 and 1 cm, the values of 
6, are 0.16, 0.1, 0.06 rad (i.e., 9, 6 and 4 deg), respectively. 
It can now be visualized that over almost the entire subsonic 
nose region, except for a very small area around the stagna- 
tion point, the destabilizing effect of the tangential field is 
even more important than that of the normal field. 

The logarithmic increment per second at the point 0; will be 
roughly equal to twice the value for either the tangential or 
the normal component of the gravity field acting alone. The 
closeness of the position 6; to the stagnation point 6 = 0 and 
the symmetry of the mean properties about the @ = 0 axis 
tend to indicate that the excitation contribution due to the 
normal component of the gravity field cannot vary substan- 
tially from the stagnation region to 6;. The decrease of the 
normal field excitation, if any, would be of the order of 6,2, 
which is negligible compared with the 100 per cent augmenta- 
tion from the excitation of the tangential field. Therefore, 
unless the kinematic viscosity coefficient v, of the liquid is 
more than doubled in the region 5 to 10 deg downstream, the 
L.I./sec is likely to increase, or at least remain substantially 
the same, from the stagnation point downstream toward 
6, The local investigation in the stagnation region fails 
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to yield any decisive conclusion whether the L.I./sec will 
increase monotonically from the stagnation point down- 
stream. It is, at least, a local minimum. Accordingly, 
under ordinary circumstances, the L.I./sec and, therefore, 
the seriousness of the instability will, in contradiction to 
one’s intuition, increase from the stagnation point down- 
stream. 

The variation of the seriousness of the instability conditions 
over the blunt nose may be further clarified by comparing the 
local conditions at the sonic point and the stagnation point. 
A lower bound of the logarithmic increment per second at the 
sonic region may be obtained from Equation [29] by simply 
ignoring the contribution from the normal component of the 
gravity field which is always positive in the subsonic nose 
region. Assuming that the sonic point lies in the immediate 
vicinity of 0S 2/4 then 


(L.1./sec) sonic 0.35 Ghsonic/Vsonic [33 ] 


By comparing the estimate of Equation [33] with that at the 
stagnation point [28], it is found that the logarithmic incre- 
ment per second in the sonic region will be larger than that at 
the stagnation point go long as 


soni ox 8 Moosic Us “ls 
< 23 [ | ~15 | [34] 


Vistag Vo stag R V; stag 


Here, the value of U.z stag has again been estimated from the 
Newtonian pressure distribution over the blunt nose and the 
normal shock relation with y taken as 1.25 ~ 1.15. If the 
kinematic viscosity coefficient of the gas at the high tem- 
perature under consideration is taken as 1 cm?/sec and U, 
as 10 km/sec ~10° cm/sec and the nose radius as 1 m or 
10? cm, then this critical ratio as calculated from Equation 
[34] will be 15 or 1.5 if hsonic is 10-1 em or 3 X 107? cm. The 
corresponding figures when the nose radius is taken as 10 cm 
will be 45 and 4.5, respectively. Thus, if the kinematic vis- 
cosity coefficient of the liquid at the gas-liquid interface at 
the sonic point is no more than four to five times that at the 
stagnation point, the interfacial instability at the sonic re- 
gion, as expressed by the logarithmic increment of the initial 
disturbances per second, will be more serious than that in 
the stagnation region under most of the circumstances. Sup- 
pose some silicone-based amorphous material, like pyrex glass 
as reported in (4), is considered with an interface tempera- 
ture at the stagnation point roughly 2800 K; the kinematic 
viscosity coefficient will not then increase more than four 
times when the temperature drops to 2200 K. Such a large 
temperature drop of the interfacial liquid mass to the sonic 
point cannot be considered impossible. Even so, the sonic 
region may still be more unstable than the stagnation region 
under practical conditions. The situation clearly depends 
greatly on the fluid properties, the scale and the geometry of 
the body nose of the specific problem under consideration, 
as indicated by Equation [34]. 

By combining the information deduced from the local in- 
vestigation in the stagnation point region concerning the 
relative seriousness of the instability at the stagnation point 
and its immediate surroundings, one cannot but draw the 
following conclusion: From the point of view of estimating 
mass ablation through this instability mechanism, it appears 
hardly justifiable to derive satisfaction from the ‘tolerable 
instability” in the stagnation region. Nor does it appear 
plausible to lay confidence on the “apparent stability” as 
may be observed from the laboratory tests without any con- 
sideration of the effect of scale and of the effect of the tangen- 
tial field. It should, furthermore, be borne in mind that the 
relative area of the unstable region involved in the sonic re- 
gion is much larger than that associated with the stagnation 
point and, therefore, should receive more careful attention. 

So far, little has been said concerning the situation down- 
stream of the sonic region. This is not because it is relatively 
unimportant, but rather because little could be done to an- 
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alyze it. The analysis of the nose region, as given previ- 
ously, clearly demonstrates the importance of the aerody- 
namic situation in determining viscous dissipation and the 
logarithmic increment per second. The Newtonian pressure 
distribution is of dubious value too far downstream of the 
sonic region. The change of pressure along the surface in the 
downstream region is of the order of the undisturbed free 
stream pressure p.., and its detailed variation depends greatly 
on the geometry of the body under consideration. It is not 
possible to make any general statements for (—U"o+/U'o+) 
as was done with Equation [25] for the subsonic nose region. 
It is not possible to claim that (—U"o+/U’o+) will be 0(1) or 
even remain positive. Thus, the simplification from Equa- 
tions [18 to 21] cannot be made. Furthermore, compared 
with the nose region, very little is known about the interface 
temperature and, therefore, the kinematic viscosity coefficient 
of the liquid »;. 

Any discussion of the situation in the downstream region 
must necessarily be made with detailed calculations for spe- 
cific instances. The crucial criterion restricting the range of 
validity of the previous qualitative discussions is the validity 
of the modified Newtonian approximation for the static pres- 
sure distribution over the nose, which may, under certain 
circumstances, be a good approximation well downstream 
into the lower supersonic region. Whether the inviscid free 
stream is supersonic or subsonic makes little difference in the 
present investigation. The “sonic region” is referred to as 
the downstream limit of the present investigation primarily 
for convenience and, secondarily, because the Newtonian 
pressure distribution often fails to be good farther down- 
stream. 


Observed Phenomena 


Fig. 2 shows four photographs of meteorites from the col- 
lection of the Chicago Museum of Natural History, repre- 
senting different categories of meteorites insofar as the extent 
of ablation is concerned. The Shelburne is stony (10 X 53 X 
3 in.) with a smooth surface (accompanied by a similar but 
slightly larger piece). The Long Island is also stony, but 
much larger (34 X 20 in. cone) smooth over the nose, but 
ablated downstream. The South Bend is stony-iron (53 X 
5 X 33 in.) with a heavily ablated surface. The Williamette 
is a much larger iron meteorite, heavily eroded. The relative 
ablation intensities of these examples may be understood 
through interfacial instability in the light of results from the 
present analysis. 

The kinematic viscosity of the molten iron is much smaller 
than that of the silicone-based material of a stony meteorite. 
The smaller kinematic viscosity as well as the larger size 
tend to indicate the relative seriousness of interfacial insta- 
bility in mass ablation in the case of the Williamette meteorite, 
even if it is assumed that the history of deceleration G(t) is 
essentially the same for all meteorites under comparison. It 
is here emphasized, however, that the present analysis con- 
cerns only the initial phase of instability, and not the fully 
developed phase of shedding droplets from spikes and the 
last phase of etching into the depth of the surfaces. Ac- 
cordingly, the inference in this case is not much more than a 
plausible conjecture. Actually, the deceleration encoun- 
tered by different meteorites will vary considerably as the 
square of their entering velocities. 

The different surface conditions exhibited by the rela- 
tively smooth Shelburne and the ablated Long Island meteor- 
ites, both of stony nature, may be explained by their differ- 
ence in size, if the entering velocity is again assumed to be 
the same. 

The most crucial, and by far the most interesting, aspect is 
illustrated by the Long Island meteorite (Fig. 2B) and is 
probably the most convincing. The Long Island meteorite 
exhibits a smooth nose of relatively small extent, followed by 
an extended but not seriously ablated region. Firstly, the 
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Fig. 2 Photograghs of typical recovered meteorites. A: Shel- 
burne, stony. B: Long Island, stony. C: South Bend, stony- 
iron, 80-20. D: Williamette, iron 


appearance of the ablated surface, resembling wavy ripples, 
tends to suggest that the ablation be related to some instability 
phenomenon of surface waves, at least in the early stages 
of development. Because the liquid layer is very thin and 
the temperature difference across it is large, the wavy shape 
of the liquid-gas interface tends to produce a similar wavy 
appearance of the solid-liquid interface. Secondly, the rela- 
tively smooth nose region and the wavy appearance in the 
downstream region tend to indicate that the interfacial in- 
stability becomes increasingly more serious as one moves 
farther downstream. This last observation is in agreement 
with the prediction made earlier when the effect of the tangen- 
tial component of the gravity field was taken into account. 
According to the conclusions of the present analysis, the stag- 
nation point over the nose is, in fact, the point of minimum 
instability locally, and the seriousness of the instability is, 
most likely, increasing toward the downstream region; the 
sonic shoulder is likely to be much more seriously affected. 
As a result, the mass ablation rate over the sonic shoulder of 
a body will be larger than that over the nose. This tends to 
flatten out the shoulder and to reduce the nose radius of the 
body. This latter effect, in turn, reduces the seriousness of 
instability over the nose. The shape and the ablation pattern 
exhibited by the Long Island meteorite are shared by many 
others. These characteristics of the recovered large stony 
meteorite need not necessarily be attributed to the incidental 
geometry and composition of the meteoric body. It may also 
be noted that the waves in Fig. 2B exhibit a three-dimensional 
character with important azimuthal component. It is indi- 
cated in (1) that three-dimensionality of the waves has no 
fundamental importance. 

From the comparison of the analytical predictions and the 
“post-mortem” observation of the meteorites, it is fair to as- 
sert the plausibility of the proposed mechanism of mass 
ablation through interfacial instability, leading to loss of 
liquid mass directly from the liquid-gas interface. The loga- 
rithmic increment per second calculated from the simple for- 
mulas in this report, may well serve as a relative measure of 
the seriousness of the mass ablation rate through this inter- 
facial instability mechanism. 

It may also appear pertinent to remark on the spectacular 
phenomena of the last phase of Sputnik II plunging into the 
dense atmosphere over the Caribbean Sea as reported by 
many incidental observers (5). The process of disintegration 
was accompanied by the breaking off of “fragments” or “large 
particles’ from the main object, and these remained distinct 
after the passing of the tail. This implies that fragments or 
larger particles must have been “projected away’ to some 
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distance from the body, under some mechanism like that of 
the interfacial instability as described in (2). The possibility 
of laminar separation of the liquid layer from the body can- 
not alone be accepted as the primary mass ablation mecha- 
nism. Any particles or object breaking away from the body 
with a small velocity component normal to the body surface 
would simply be swept downstream along the body surface 
and could not possibly appear as distinctly brilliant sparks 
away from the main body brilliance, especially when viewed 
from great distances. 

{t should be pointed out, in addition, that the excessive 
miss ablation over the sonic shoulder of the blunt nose may 
lovally weaken the structural strength to such an extent as 
to cause the actual mechanical failure or breaking up of the 
be ly into large solid pieces. The phenomena are, indeed, 
very much involved. It is not suggested that interfacial 
in-tability alone will account for every detail of the observed 
phenomena. The primary purpose of the present investiga- 
tivn is to indicate the plausibility and the likelihood of the 
proposed mechanism of mass ablation, which has, so far, led 
tc predictions in no way contradicting observations. 


Conclusions 


The presence of the tangential component of the equivalent 
gravitational force field and of the variation of the fluid dy- 
namical situation of the shear flow on the surface of a blunt 
nosed body flying at hypersonic speed, considerably modify 
the relative importance of the physical variables affecting the 
seriousness of the interfacial Taylor-type instability as a 
mechanism of mass ablation. The main results of significant 
practical importance are qualitatively: 

1 The seriousness of the interfacial instability, as ex- 
pressed by the logarithmic increment of the most critical 
component of the initial disturbances, increases with the linear 
scale R of the geometrically similar body nose, everything 
else being the same. ° 

2 The stagnation region over the nose need not be the 
region most seriously affected by interfacial instability. It is 
more likely to be much less seriously affected than the 
sonic region downstream. The situation depends greatly on 
the temperature sensitivity of the kinematic viscosity coeffi- 
cient of the liquid at the gas-liquid interface. It appears that 
the temperature sensitivity of the viscosity coefficient of the’ 
liquid should be considerably larger than that of the silicone- 
based ordinary amorphous solids, in order that the stagnation 
region may be more critical on bodies of practical scales. 


Quantitatively, the logarithmic increment may be split into 
two parts, one due to the normal component and the other 
due to the tangential component of the gravity field (Eqs. 
[28 and 29]): 

1 That due to the normal component is proportional to 
R'’: and essentially independent of the film thickness h. 

2 That due to the tangential component is independent 
of the linear scale R but linearly proportional to the film 
thickness h. 

3 The tangential field contributes more than the normal 
field in exciting instability in the region with 06> @,, which 
is of the order of 0.1 radian from the stagnation point, as 
given by Equation [32]. 

4 The sonic region is more seriously affected by interfacial 
instability than the stagnation region if the ratio of the kine- 
matic viscosity coefficients of the liquid is less than the 
quantity proportional to a Reynolds number as specified in 
Equation [34]. 

Although the normal component of the gravity field does 
not modify the stability boundary of the ‘‘propagating dis- 
turbances” c, ¥ 1 for the instability of the Helmholtz type, 
the tangential component does. 

On the surface of the far downstream region of the body, 
where the Newtonian pressure distribution fails to be a fair 
approximation, the situation of the interfacial instability will 
have to be determined by detailed calculations under specific 
geometry and environment conditions. 

Examination of the surfaces of the recovered meteorites 
and the reports from the incidental observers who witnessed 
the last phase of Sputnik IT plunging into the lower atmos- 
phere, tend to substantiate the plausibility and the likeli- 
hood of the proposed mass ablation mechanism through 
interfacial instability. 
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Technical Notes 


Graphical Computation of Shock and 
Detonation Waves in Real Gases’ 


GEORGE H. MARKSTEIN? 
Cornell Aeronautical Laboratory, Inc., Buffalo, N. Y. 


The conservation relations across a shock or detonation 
wave are represented in the form of a graph applicable to 
real gases of arbitrary composition. Computation of 
shock waves reflected from a closed end is performed by 
means of another universally applicable graph. 


ETHODS for computing the state of a real gas behind a 
shock wave or detonation wave have been discussed 
extensively in the literature, and comprehensive reviews of 
the subject have been published recently (1,2). The nu- 
merical methods usually employed require time-consuming 
iteration. The present graphical procedure, similar to that 
described by Evans (3), avoids this drawback and, moreover, 
may be applied to gases of arbitrary composition. 
Transformation of the shock relations to a form that is 
independent of the individual properties of the gas is achieved 
by introducing the following dimensionless variables, which 
apparently were first used by Skinner and Squire (4) 


r = p/p [1a] 

P = p/p [1b] 

wT; [ie] 


= — = = 1/2 ld 


x= hy) [le] 


Note that in order to eliminate 7: from the shock relations, 
the definition of dimensionless velocity differs from the usual 
formulation in terms of Mach number. (It is normally not 
necessary to consider deviations of the compressibility factor 
Z = up/pRT from unity, but this may be done readily by 
incorporating a factor Z/Z, in Eq. [lc].) 


Incident Shock 


In terms of the variables [la-le] and with the notation 
indicated in Fig. 1, the equations of conservation of mass, 
momentum and energy and the equation of state for the inci- 
dent shock assume the form 


r(U; — U2) = U; [2a] 
P, wae 1 = [2b] 
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Fig. 1 Notation for shock reflection from a wall 


X% = — (1/2)U2] 
P, = 7202 


The additional information required is the relation between 
enthalpy and temperature X = X(@). This is the only rela- 
tion that depends specifically on the composition of the gas. 
It is therefore convenient to represent the relations [2a—2d | 
graphically in the 6.-x2 plane, in the form of families of curves 
that correspond to constant values of any one of the other 
variables. The P.» = const. and rz = const. curves are 
straight lines; the U; = const. curves are parabolas, and the 
U. = const. lines are hyperbolas. A plot of the two families 
U; = const. and P; = const., that are usually of primary 
interest, is shown‘ in Fig. 2. 

Shock computation is carried out by determining the X(6) 
curve for a given gas mixture and initial conditions from 
thermodynamic tables. One may then obtain the state of 
the gas behind a shock wave of velocity U; or pressure ratio 
P», from the intersection of the X(@) curve and the respective 
curve of Fig. 2, using Equations [2] for computing the vari- 
ables not contained explicitly in the graph, and Equations 
[1] for converting the results into dimensional values. 

For perfect gases, the X(@) curves are straight lines through 
the point 6 = 1, x = 0 of slope y/(y — 1), ie., 2.5 for 
a monatomic gas, 3.5 for a diatomic gas, etc. For real gases 
undergoing relaxation reactions behind a shock of velocity 


4 Sets of graphs for incident and reflected shocks drawn to a 
scale suitable for computation may be obtained from the author. 
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Fig. 2 Graph for incident shock computation 


UL, the conservation relations apply not only to the “‘frozen”’ 
state immediately behind the shock but also to any subse- 
quent state in the relaxation zone, if the flow may be regarded 
as steady in shock coordinates. Therefore, the curves U; = 
const. determine the reaction path in the 6,-X, plane (3). 


Dissociation Behind a Shock 


As an example of the application of the method to a real 
gas, the case of a diatomic dissociating gas is shown in Fig. 3. 
Yor this ease, the relations [lc and le] become (3) 


@=(1 +a) [3a] 


1 
RT, 


Typical enthalpy curves for undissociated gas, a = 0, fully 
dissociated gas, a = 1, and equilibrium composition at con- 
stant pressure, a = a,(T’, p), are shown in Fig. 3. Since the 
fully dissociated gas is monatomic, the corresponding X(@) 
curve must be a straight line of slope 2.5. The curve for the 
undissociated gas always has a larger and gradually increas- 
ing slope, and therefore the two must intersect at some point 
C that normally lies at very high temperature. As a conse- 
quence, as discussed by various authors (1,3), dissociation 
relaxation at shock velocities below the value corresponding 
to point C is accompanied by increases of pressure and 
density, as from A to B, whereas at values above that critical 
velocity, as from D to EL, the pressure and density decrease. 
At the point C the only effect of dissociation is a temperature 
drop by a factor 1/(1 + a) = 3. 


Detonation 


The application of Fig. 2 to detonation is shown in Fig. 4. 
Since for an exothermic reaction, the X(@) curve of the prod- 
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Fig. 3 Application of the graph shown in Fig. 2 to computation 
of a shock wave in dissociating diatomic gas 
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Fig. 4 Application of the graph shown in Fig. 2 to computation 
of a detonation wave 


ucts at eqilibrium must lie below that of the reactants, 
the U; = const. curves for U; smaller than a critical value 
Up do not intersect the enthalpy curve of the products. 
(Note on Fig. 2 that all the U; = const. curves go through 
the point 6. = 1,X,=0.) For U; = Up, the enthalpy curve 
is tangential to the U; = const. curve, whereas for U; > Up 
there are two intersection points. Differentiation. of Equa- 
tions [2a and 2b] shows that along the U; = const. curves the 
relation 


(Op2/Op2)u; = (ui — [4] 
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Fig. 5 Graph for reflected shock computation 


is satisfied. Since along the X(@) curve of the products the 
gas remains in thermodynamic equilibrium, the value of 
(Op2/Op2)v; at the point of tangency must be equal to the 
square of the equilibrium speed of sound. This point there- 
fore satisfies the Chapman-Jouguet condition, and the velocity 
Up is the detonation velocity. The ‘von Neumann spike” 
corresponds to the state of the reactants behind the shock 
wave of velocity Up, as indicated in Fig. 4. The two inter- 
section points of curve U; = const. > Up with the X(@) curve 
of the products correspond to so-called weak and strong de- 
tonation. The present method requires only the plotting 
of the x(@) curve of the products on the graph, Fig. 2, and 
yields the Chapman-Jouguet condition automatically with- 
out computation of the equilibrium sound speed. It should 
therefore provide a considerable saving of labor compared 
with conventional numerical methods (2). However, for 
accurate computations, a graph with much smaller steps in 
U; than shown in Fig. 2 would be required. 


Reflected Shock Wave 


In terms of the variables [la—le], the relations for the 
reflected shock, Fig. 1, become 


r3U, = r(U, + U2) 
P; (U, + U2)r2U2 


These relations cannot be represented as simply as those for 
the incident shock, because r. and U, are uniquely determined 
only if a pair of values U;, P: is given. Thus, to each inci- 
dent shock velocity U; corresponds not a single curve but a 
two-dimensional] region in the 6.-X; plane, in which a family 
of curves P; = const. may be plotted. Fig. 5 shows a portion 
of this graph. For a gas of given X(@), the incident shock 
computation yields the values of P2 corresponding to each 
U;. The state behind the reflected shock is then found as 
the intersection between the (interpolated) P. = const. 
curve and the X(@) curve, plotted in Fig. 5. ; 
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If relaxation reactions occur behind the reflected shock, 
the choice of the X(@) curve to be used for computation de- 
pends, of course, on the rate of the reactions. If relaxation 
is so fast that the gas in region 3 may be regarded essentially 
as being at equilibrium, the equilibrium curve must be 
used. Conversely, the “frozen” composition curve must he 
used if the reactions are slow, as required, for instance, in thie 
single-pulse chemical shock tube (5). The paths of relax- 
ation reactions behind the reflected shock are given by the ccn- 
dition U, + U, = const. Families of these curves may he 
plotted in the regions of interest on the 6;-X; plane for each 
U; = const., if needed. 


Nomenclature 


speed of sound 

heat of dissociation at 7’ = zero K 
enthalpy per unit mass 

molar enthalpy 

heat of reaction at 7 = zero K 
pressure 

dimensionless pressure, p/p: 

universal gas constant 

dimensionless density, p/p: 

absolute temperature 

velocity 

dimensionless velocity, (p:/p1)!/2u 
compressibility factor 

degree of dissociation 

ratio of specific heats 

dimensionless enthalpy, (p:/p:)(h — hi) 
molecular weight 

dimensionless temperature, (4:/u)(7'/T;) 
density 


NSE Q* 


Subscripts 

initial state 

state behind incident shock 
state behind reflected shock 
incident shock 

reflected shock 

atom 

molecule 

detonation 

equilibrium 
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Magnetic Boundary Layers 


H. F. LUDLOFF' 
New York University, New York, N. Y. 
HE SUPERSONIC motion of a conductive fluid of mag- 
netic Reynolds number m = 2mra0Uc and Mach number 


M°, about a thin airfoil contour h(x), carrying a relatively 
weak magnetic field H,, Hy, is determined by the equations 


— 1)gqzz — Wy = GHz, H,) 


Received May 25, 1959. 
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where 


y = disturbance velocity vector (u, v) due to the presence of 
H as well as h(x) 
G. = polynomial of second degree in the derivatives of Hz,H, 


The magnetic flux satisfies 
oH 
AH — 2m-— =0 [2] 
or 
All distances are dimensionless, measured in units of chord 
length ce. 


After introducing characteristic coordinates —&, 4 = x + 
8- into Equation [1], and fulfilling the boundary conditions 


at » = 0,v = U-h’(x) and at & = — ~, g = 0, we obtain as 
solution of Equation [1] in the upper half plane 

v = GE, — Gn, 0) + h(n) [3] 
where 


Gé,n) = dn-Gé, 7) 


This process is the same for u, and for the lower half plane. 
G van be computed after H has been determined from Equa- 
tion [2]. 

if we require as boundary condition for Equation [2] that 
H vanish at infinity, and behave like the field of a current- 
carrying wire, near the wire, then we obtain by use of a vector 
potential A, 


H, = —m-e™*.(y/r)-K,(mr) 


H, = —m-e""-Ko(mr) + 
where K, is the modified Hankel function. 
For small mr 
si s 6 
5) 


This means the magnetic lines of force are circles around the 
origin (wire) for vanishing magnetic Reynolds number (i.e., 
without wind). 

For large mr 


H, ~ +... 
H, ~ .(1 — cos 


In this case, the magnetic lines of force are given by 
A, = e™-Komr) =c 


They are, in approximation, narrow parabolae about the z-axis: 
r — x = constant. The magnetic field, which was originally 
axisymmetric, is now “blown downstream” by the relative 
wind and forms a thin “magnetic boundary layer” along the 
x-axis. The H-lines are closed curves, also on the downstream 
end (see Fig. 1). Within such an egg-shaped eurve lies the 
essential strength of the field; the outside field is negligible. 
From Equation [6], it follows that the field strength H has 
fallen off to less than 1 per cent of its value in the vicinity of 
the chord, within approximately 25 chord lengths along the 
r-axis on the downstream side; on the upstream side, its ex- 
tension is practically zero. The width of the field in the 
cross-direction is indicated by 


H? ~ 


which means that the magnetic “boundary layer thickness” 


is 
mr 


a relation which is similar to that in viscous layers. The field 
profile inside the layer is shown in Fig. 2; there are other 
varieties. 
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Fig. 3 Velocity disturbances 


The influence of the magnetic field on the velocity profile 
can be deduced from Equation [1]. Clearly, due to the non- 
homogeneity of the wave Equation [1], the velocity dis- 
turbances are not only caused by the magnetic field along 
the same characteristic, as in ordinary gasdynamics, but are 
a cumulative effect of the whole field, extending upstream 
between the characteristics & = m and & = £&, which in some 
cases may be considerable. (See Fig. 3.) 

By use of the preceding analytic expressions, the trans- 
formation of an originally axisymmetric field into a thin, 
long, magnetic boundary layer, due to the action of a strong 
relative wind, can be described in all details; this insight may 
be regarded as supplemental illustration to the often-quoted 
theorem that in a conducting fluid, the magnetic flux is partly 
carried with the fluid and partly dissipated through it. Note 
that these magnetic layers exist also in a free stream, if no 
disturbing body contour is present. The resulting velocity 
profiles are qualitatively different from those in ordinary 
gasdynamiecs in that the disturbances are zero on the sym- 
metry line and increase toward the outside and in the down- 
stream direction, as seen from Fig. 3. 

A more general study of magnetic boundary layers is to be 
based on the nonlinear equations of motion. 
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Maximum Ejection Velocity for Return 
From Satellite Orbits 


HARRY ESSES' 


General Applied Science Laboratories, Inc., 
Hempstead, N. Y. 


E PROBLEM considered concerns the minimum ve- 
locity increment required to transfer an Earth satellite 
from its initial orbit into an orbit which intersects the Earth. 
This velocity increment will be referred to as the ejection 
velocity. In particular, it will be proved that the ejection 
velocity has an upper bound »,; = 4930 fps at an altitude of 
19,325 statute miles or about 4.88 Earth radii. 

The interesting feature of this result is that the maximum 
occurs at a finite distance from the surface of the Earth, so 
that the ejection velocity first increases with altitude, reaches 
a maximum, and then decreases with altitude. At first 
glance, this result may be surprising, since one might expect 
the ejection velocity to increase monotonically with altitude, 
where a greater ejection velocity would be required for a satel- 
lite at an altitude of 50,000 statute miles, for example, than 
for one at an altitude of 25,000 statute miles. A simple 
calculation shows, however, that the ejection velocity for the 
50,000-mile satellite is less than that for the 25,000-mile 
satellite. 

The strangeness of this result vanishes, however, as soon as 
one notes that the velocity required to establish a circular 
orbit decreases monotonically with altitude, from about 26,000 
fps at the surface of the Earth to zero “at infinity.” It then 
becomes apparent that there exist altitudes for which the 
circular orbital velocity itself is less than the ejection velocities 
required for certain lower altitude orbits. Insofar as it is 
not necessary to remove the entire orbital velocity to cause 
a given orbit to intersect the Earth, it is at least intuitively 
clear that there exists a turning value for ejection velocity 
which occurs at a finite distance from the Earth’s surface. 
A formal proof of this statement follows, which is based upon 
orbital motion about a spherical Earth, in the absence of the 
atmospheric envelope and all other perturbating forces. 

The apogee velocity for a satellite traversing an elliptical 
orbit is well-known to be 


ve = ‘) [1] 


a\l +e 


= semimajor axis 
e = eccentricity 

go = acceleration due to gravity at the surface of the Earth 
Ro= radius of the Earth 


The ratio of perigee radius r, to apogee radius r. is 


Tp (; —e 
re \l+e (2) 
and the semimajor axis is 
Tp + f. 


[3] 


The mathematical requirement for an elliptical orbit to 
intersect the Earth is that the perigee radius be equal to the 
radius of the Earth. Inserting Equations [2 and 3] in Equa- 
tion [1] withr, = Ro gives 


Ye = [4] 
Ta(Ta + Ro) 


Received May 20, 1959. 
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The circular velocity for a satellite at an altitude h = r. — Ry 
is 


5) 


Ta 


= 


Clearly, the difference (v. — va) is the ejection velocity, or 
the velocity decrement required to reduce a given circular 
orbit into an elliptical orbit which intersects the Earth. 
Moreover, insofar as any elliptical orbit can be inscribed in a 
circular orbit with a radius equal to the apogee radius of the 
ellipse, the ejection velocity for any orbit will never exceed 


Vey = (Ue — Va) = (: 6} 


The altitude at which the ejection velocity assumes its maxi- 
mum value is obtained in the usual way by setting the ce- 
rivative of Equation [6] with respect to ra, equal to zero. 
This results in a cubic equation for r., given by 


ra? — 5Ror,? 5Ro*rg = 0 


or, withr, = KR, 
— 5K? — 5K —1=0 7] 


Equation [7] has one real root, K = 5.88, so that the maxi- 
mum value for v,j occurs at 23,285 statute miles from thie 
center of the Earth or at 19,325 statute miles from the surface 
of the Earth. It then follows from Equation [6] that tlie 
minimum ejection velocity required for any satellite or)it 
never exceeds 4930 fps. 

These numerical values agree with similar results reported 
without analysis in a Rand Corp. report? which, incidently, 
gives a very fine treatment, in lay terms, of almost the entire 
field of space technology. 


2 “Space Handbook: Astronautics and Its Applications,” 
Staff report of the Select Committee on Astronautics and Space 
Exploration, prepared by the Rand Corp., Santa Monica, Calif. 


Achieving Satellite Rendezvous 


RICHARD J. WEBER! and WERNER M. PAUSON’ 
Lewis Research Center, NASA, Cleveland, Ohio 


ANY proposed space missions will require achieving 

rendezvous of two bodies in an orbit about a planet. 
For example, an interplanetary transport vehicle may wish 
to transfer passengers to an orbiting space station. It is 
generally most efficient to control the spaceship during its 
approach to the final orbit altitude in such fashion that 
rendezvous is accomplished immediately upon entering orbit. 
(Methods of transferring between dissimilar orbits with simul- 
taneous rendezvous are discussed in (1 and 2).*) However, as 
a result of guidance errors, unaccounted-for perturbative 
forces, ete., it is unlikely that perfect rendezvous will be 
achieved in this fashion. The situation may then be envi- 
sioned where the spaceship has established itself in the de- 
sired orbit (assumed to be circular), but is separated from the 
target body by some angle. 

Several procedures are conceivable for correcting this 
angular error. With high thrust engines, it is possible for the 
spaceship to be transferred by means of an impulsive thrust to 
an elliptical trajectory that returns to the circular orbit at 
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Fig. 1 Forces acting on spaceship 


the proper instant for rendezvous (2). With low thrust 
engines (such as ion or plasma rockets), however, thrust must 
be applied continuously during an appreciable time interval 
in order to change the velocity or path of the spaceship. One 
expedient is to continuously direct the thrust in a more or 
less circumferential direction in order to spiral some distance 
away from the circle. The distance is chosen such that, upon 
spiraling back to the circle, rendezvous is achieved. 

In this note an alternative method is considered in which 
the spaceship remains in the circular orbit while it accom- 
plishes rendezvous with the target. The velocity increments 
and times required for this method are compared with those 
for the techniques mentioned above. 


Analysis 

The initial thrust applied to the spaceship is directed cir- 
cumferentially. As the velocity of the spaceship changes, 
the thrust vector must be rotated so that its radial component 
compensates for the resulting imbalance between the centrif- 
ugal and gravitational forces in order to maintain the space- 
ship in the circular orbit (Fig. 1). 

If it is assumed that the thrust-to-weight ratio of the 
spaceship is constant during the rendezvous maneuver, the 
differential equations of motion may be solved explicitly for 
this method of flight, vielding the relations 

Qa = a — 20 [1] 
d9 
ail, V1 + (F/W) sin 20 

Initially, @ is zero deg if the target is ahead of the spaceship 
(180 deg if the target is behind). Equation [1] shows that 
the thrust angle varies linearly with the distance traveled. 
Note that upon completion of 45 deg of spaceship travel, the 
full thrust vector is oriented radially in order to hold the 
spaceship in the circular orbit, and the vehicle can no longer 
be accelerated. 

Equation [2] can be transformed into an elliptic integral of 
the first kind and evaluated with the use of tables. In the 
period of time At, the spaceship and the target travel through 


At = 
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Fig. 2 Maximum possible decrease in angular separation 
during acceleration and deceleration phases as a function of 
thrust-to-weight ratio 


-3 ORBIT 
28x10 ALTITUDE, 
MI 
500 
Qa 
a 1000 
z 
wn 
S 
a 
aan 
oO 
oO 
<a 
fe) 2 4 6 8 


THRUST-TO-WEIGHT RATIO 


Fig. 3. Rate of closure during coasting as a function of thrust- 
to-weight ratio for two orbit altitudes 


the angular distances @ and v.At/r, respectively. The differ- 
ence between these displacements is the amount by which 
the initial angular separation 6) is reduced. The thrust vector 
at time At is at an angle a, with the circumferential direction. 
To restore the spaceship to the circular orbit velocity, the 
thrust vector is reversed (mirrorwise) to an angle of 180—a; 
deg and then rotated back to a circumferential position, 
ending, however, opposite to the initial thrust direction. 
The further reduction in angular separation during this 
deceleration phase is equal to the reduction accomplished 
during the initial acceleration period; also, each phase of the 
maneuver takes time At. The maximum amount 6, by 
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Fig. 4 Rendezvous time as a function of separation angle for 
several values of F/W. Altitude = 500 miles 


which the initial angular separation 6) may be reduced during 
the acceleration period plus a similar deceleration period is 
shown in Fig. 2. 

If the initial angular separation is greater than the value 
shown in Fig. 2, the thrust may be held in the intermediate 
radial position and the spaceship allowed to “coast” at con- 
stant velocity for a long enough period of time to make up the 
difference between 6) and 6m. The rate of closure during the 
coasting phase is given by Equation [3] and plotted in Fig. 3 


_ F 


As a measure of the fuel expenditure necessary to carry out 
a rocket mission, it is convenient to consider the character- 
istic velocity increment that would be achieved by the 
vehicle in linear gravitation-free flight. This velocity incre- 
ment may be calculated for the constant radius rendezvous 
maneuver from Equation [4] 


AV = f [4] 
0 m 


where trotar = 2At + teoast- 

The times and velocity increments required to correct for 
various angular separations are shown in Figs. 4 and 5. Note 
that both of these quantities are quite large for all but the 
smallest values of 6) in the range of thrust-to-weight ratios 
anticipated for continuous thrust propulsion systems. 
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Fig.6 Velocity increments required for rendezvous as a function 
of rendezvous time. Altitude = 500 miles 
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Discussion 

A comparison among the various methods described for 
attaining rendezvous is made in Fig. 6 for several separation 
angles. (The velocity increments and times shown are for 
separation angles in which the target leads the spaceship.) 
The constant radius technique is definitely inferior to the other 
methods in terms of both the time and the velocity increment 
required for rendezvous. For a separation of 0.25 deg (which 
is only about 20 miles at an altitude of 500 miles) a AV of nearly 
300 fps is necessary for a thrust-to-weight ratio of 10-4. With 
sufficiently high values of specific impulse, the value of AV 
needed for rendezvous may not correspond to a very large fuel 
consumption; however, the spaceship crew would not care to 
take the indicated time of 24 hr to complete the last 20 miles 
of their journey. 

‘The constant radius method of achieving final rendezvous 
does not, therefore, appear desirable unless available thrust- 
to-weight ratios are much higher than 10 ~‘, and then only if 
further study indicates some advantage (e.g., guidance) over 
the spiral trajectory method. 


Nornenclature 


F = engine thrust 

m = spaceship mass 

r = distance from center of Earth 

t = elapsed time 

v. = circular orbit velocity 

W = weight of spaceship 

a = angle between thrust and velocity vectors 
6 = angular separation between the target and spaceship 
6 = angular distance traveled 

nh = gravitational force constant for Earth 
References 


1 Paiewonsky, B. H., ‘‘Transfer Between Vehicles in Circular Orbits,” 


WADC Rep. TN 57-267, Aug. 1957. 
2 Brunk, W. E. and Flaherty, R. J., ‘‘Methods and Velocity Require- 
ments for the Rendezvous of Satellites in Cireumplanetary Orbits,’’ unpub- 


lished NASA data. 


Stability of a Globular Cluster in the 
Oblate Earth’s Equatorial Plane 


MARVIN H. HEWITT! 
The Martin Company, Denver, Colo. 


An expression is obtained for the stability of a globular 
cluster whose center of gravity describes a circle in the 
oblate Earth’s equatorial plane. It is found that the sta- 
bility criterion is not significantly changed if a spherical 
Earth is assumed. 


HE CENTER of gravity of an aggregation of particles 
describes a circle in the equatorial plane of the oblate 
Earth. Let C denote this center, and let E be any element of 
the cluster. Also, r and @ will denote the equatorial polar co- 
ordinates of C referred to the Earth center 0, and &, 7 are the 
coordinates of E referred to C as origin, the axis of £ being the 
prolongation of OC. Let M be the mass of the Earth. Neglect 
the Earth’s atmosphere and assume that the cluster, as a 
first approximation, is homogeneous and spherical. 
The magnitude of the globe’s attraction at an internal point 
is given as 


pz = —Cmz/b {1} 
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where 
m, b = mass and radius of the globe, respectively 
z = CE 
We obtain? 
—Gm/b' [2] 


The approximate potential function for the gravitational 
attraction of a solid homogeneous oblate Earth upon a dis- 
tant particle in the equatorial plane may be given as? 

VU M(a? — c?) 


Us -— 3 
10x? [3] 


a = semimajor axis of the Earth 
¢ = semiminor axis of the Earth 
zx = distance from the Earth’s center to distant particle 


The corresponding acceleration g is then given by the expres- 


sion 
1 3(a? — c?) 
= -GM|-— + ——_— 4! 
g ] |4! 


It is seen that this acceleration represents a central non- 
Newtonian force. 
The equations of motion are given now as 


l ) 1d d6 


dt? a) nat 
[5] 
€") 
—-GM 
(r + &)! | 
[6] 


3n(a? — c?) 


10(r + 


These equations also hold for the motion of a cluster element 
at C, where € = 0, 7 = 0. Consequently, when Equations 
[5 and 6] are expanded in powers of & and 7, all the terms in- 
dependent of — and yn must cancel. We obtain 


me an dd (2 
dt dt de 


dt" dt de ap 
1  3(a? — c?) 
1075 


Since the center of gravity of the cluster describes a circle 
about the Earth, we have r = d (constant), dé/dt = n (con- 
stant). Centripetal acceleration in the oblate Earth’s equa- 
torial plane is of the form 


y? 1 3(a? — c*) 
dn? = GM E | 


The centripetal acceleration in any orbital plane of the spheri- 
cal Earth is given by the expression 


= dno? = —— [10] 


2 Webster, A. G., “Dynamics of Particles and of Rigid, Elas- 
tic, and Fluid Bodies,” G. E. Stechert and Co., N. Y., 1942, 3rd ed. 
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where 
a 
e = 
(7] 
2 
-4 
000 
tion 


where np = d6/dt. Consequently 
3MG(a? — c?) 


(11) 


n? = No + 


Introducing Equations [10 and 11] into Equations [7 and 8] 


yields the equations 
+ — 5n* — p) = 0 [12] 
dt? 

dé 

ae Wan? + = 0 [13] 
Let 

— = A cos (pt + a) [14] 
= B sin (pt + a) [15] 


where A, B and a are real constants. We obtain from Equa- 
tions [12 and 13], after eliminating A and B, the determinan- 
tal equation 
(p? + 2n? + w)(p? + bn? — 2m? + — 4n%p? = [16] 
If p? is real and positive, the globe will be stable. This con- 
dition will be satisfied if 


|p| > 5n® — Qo? {17] 
or 
3MG(a? — c? 
|u| > Bro? + {18] 
We may write Equation [18] as 
m 3m? 3M(a* — c*) 


Let D denote the density of the globe. The condition that the 
globe is stable is therefore 

9M 9M(a? — 
4rd’ 


D> [20] 


It is of interest to compare the first and second terms of 
Equation [20]. We have 


9M(a? — c?)/8md* 


(21), 
i A stable globe in the equatorial plane at a central distance 
3 equal] to four times the equatorial radius would be required to 
1 have a density greater than 0.05 gm cm~, whereas a stable 


equatorial globe at a central distance equal to 11/10 the equa- 
torial radius would have a density greater than 11 gm cm=%, 
which is about the density of lead. 


Optimum Radiator Temperature for 
Space Power Systems 

7 EDWARD T. PITKIN! 

” Marquardt Aircraft Co., Van Nuys, Calif. 


NUMBER of articles have recently appeared (1 to 3)? 
in which the optimum radiator temperature for a power 
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In (1), the de- 
sired radiator temperature to obtain maximum overall effi- 
ciency of a solar powered converter is discussed, while (2 and 
3) are concerned with minimum radiator weight for a given 


converter to be used in space is discussed. 


power output. The purpose of this note is to discuss the 
radiator temperature for minimum system weight when 
power source weight varies as a function of power output, 
and the radiator weight is a function of temperature. 

Consider an energy conversion system which takes therm: 
energy from a heat source, converts a portion to available 
mechanical or electrical energy, and discharges the unavail- 
able energy through a radiative heat sink. This analysis 
will apply to any system where the useful power output /’, 
is equal to the product of a conversion efficiency factor 7, tlie 
Carnot efficiency ». and the thermal power output of the 
source P, 


P, = = [1] 
where 


This relation will hold for most systems which employ re- 
circulating working fluids with conventional heat engine 
cycles. It is further assumed that the energy source weig)it 
is proportional to the thermal power output raised to an ex- 
ponent, and that the radiator weight is proportional to its 
surface area. Other weight-producing components of tlie 
system will be motors, generators and connecting structure. 
The total weight of these components may be assumed di- 
rectly proportional to the useful power output as a first order 
approximation. 
The energy source weight may then be expressed 


P./n. |" 
Ws = W, W, [3 
) (P./m 


where 7:, the thermal efficiency, is the product of the Carnot 
efficiency and the conversion efficiency. Let the reference 
condition denoted by subscript 1 be the case of 100 per cent 
thermal efficiency, and divide by an arbitrarily fixed power 
P., while noting that this power is identical to the thermal 


Sep of the source at the reference condition (i.c., P, = 


Ps Py) 


An energy balance on the system gives the heat rejected by 
the radiator 
1 — 
P, = ——' p, [5] 
Ne 
This is equated to the heat removed by radiation by use of 
the Stefan-Boltzman relation while assuming the radiator 
surface and fluid temperatures are essentially the same. Also, 
the heat sink is assumed to be at absolute zero temperature 
as a reasonable approximation for the space environment 


P, = A,oeT,* [6] 


If a temperature gradient exists over the radiator surface, 
Equation [6] may be interpreted as defining an appropriate 
mean temperature. This might apply to systems utilizing 
noncondensing fluids. The radiator area is easily determined 
as 

(1 


A, = 

Let K, be the proportionality constant between radiator 

weight and area, and K, be the constant relating weight of 

motor, generator and connecting structure to power. The 

weight to power ratio of the conversion system then becomes 
W Wa 1 K,(1 — 


P, Pa cen: T K, 8) 
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Fig. 1 Effect of converter efficiency on optimum radiator 
temperature 


Differentiating with respect to 7, at constant power to 
obtain a minimum gives 


Pa oeT Ol, UG 


_aW/P.) 


oT, 
Also remember that 


T, 
- 7) [10] 


The conversion efficiency 7 is not likely to be a function of 
T, so 


[9] 


= 0 


[11] 


Substitution of Equations [10 and 11] into Equation [9] and 
algebraic rearrangement yields an implicit expression for the 
radiator to maximum temperature ratio 7,/T,, which gives 
minimum system weight 


Wa _ T, \¥ 


Equation [3] may now be used to show 


Wa W, 
Py P, [n(: [13] 
This is then substituted into Equation [12] to yield 


Ws o€TmtN _ ( Tm 


This equation is best solved for 7,/7', by plotting the right- 
hand side for given values of 7',/T,, and ». Typical values of 
yptimum temperature ratio are plotted in Fig. 1 vs. the 


[14] 
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Fig. 2 Weight characteristics of typical gas-cooled reactors 


parameter given in the left-hand side of Equation [14] for 
different values of efficiency factor 7. The independent vari- 
able is composed of factors which can initially be estimated 
quite closely in a preliminary design study. As an example 
take N = 0.2, and assume the following values of the other 
parameters: 7,, = 2000 R, « = 1.0, K, = 1 lb/ft? and W,, P, = 
10~? lb-hr/Btu. These values combine such that Fig. 1 
is entered at an abscissa of 55. Then 0.40 < T,/T,, < 0.48, 
depending on the conversion efficiency. The Carnot effi- 
ciency can then easily be determined from Equation [2], and 
the power output is established by Equation [1]. 

It is interesting to look at the case where N = 0, i.e., power 
source weight is independent of output. In this case, mini- 
mum system weight corresponds to minimum radiator area, 
the case treated in (2 and 3). Letting N go to zero in Equa- 
tion [14] gives 


4 - - |-1=0 [15] 
This can be rearranged as a quadratic in 7/7’, 
The solution is 
Tr _ —5 + — 


8n 


where only the positive sign is retained, since negative values 
of temperature will otherwise occur. Substitution of effi- 
ciency values 0 < n < 1 into Equation [17] will show that 
the radiator to maximum temperature ratio for minimum 
radiator surface area varies between the relatively narrow 
limits 0.80 < T,/T,, < 0.75. Note that Equation [17] be- 
comes indeterminate for 7 = 0, requiring use of L’Hospital’s 
rule to obtain the limit. This numerical result has been pre- 
viously stated in (2) without a supporting derivation. 
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The value of N for a given design study must be deter- 
mined from the power source weight characteristics. As an 
example, the thermal power to weight characteristics of some 
typical nuclear reactor designs taken from (4) are plotted in 
Fig. 2. The proper value of N for use in any particular de- 
sign study may be obtained from the slope of the thermal 
power to weight characteristic of the power source in the 
power range of interest. Slopes for different values of N are 
shown on the figure. It may be expected that in the case 
of nuclear power sources, N will approach zero with low 
power reactors where criticality requirements govern size. 
In very large reactors, one may expect heat transfer surface 
area, which is proportional to thermal power output, to 
govern weight. In this case N could approach 1. 


Nomenclature 

A = area, ft? 

K = weight constant, lb/ft? and Ib/Btu/hr 

P = power, Btu/hr 

N = weight exponent 

T = temperature, R 

W = weight, lb 

o = Stefan-Boltzman constant, 0.173 K 10~° Btu/hr-ft?-R‘ 
€ = emissivity 

7 = efficiency 

Subscripts 

1 = reference condition of 100 per cent thermal efficiency 
c = Carnot 

g = generator + motor 

m = maximum 

o = output 

r = radiator 

s = source 

t = thermal 
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Optimum Proportioning of Two 
Propellants to Obtain Maximum 
Burnt Velocity 
R. H. OLDS! 

Lockheed Aircraft Corp., Palo Alto, Galif. 


Theory indicates that certain combinations of two pro- 
pellants have greater propulsive capability than either pro- 
pellant alone. Optimum performance results from a 
suitably programmed variation of propellant proportions 
throughout burning. Results applicable to the first stage 
of a rocket are presented. 


HE USE of variable propellant characteristics during 
burning has potentialities that have not yet been exploited 

to any great extent. Lorell and Hibbs (1)? by variational 
methods demonstrated benefits resulting from the use of 
varying proportions of two propellants during burning. The 
Received May 13, 1958. 
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variational problem requires consideration of four cases. In 
(1), two of the cases were treated. This paper reports results 
from an extension of the analysis to include all four cases. 
Continuous variation of the relative proportions of two 
propellants during burning is schematically represented by 
Fig. 1 in which z corresponds to the volume fraction of total 
propellant consumed, and y designates the volume fraction 
of the more dense propellant in the burning mixture. I) 
actual practice a somewhat more sophisticated means of com- 
bining the two propellants would probably be required. 


Basic Equations 

The evaluation of performance is based upon the motion c/ 
a rocket in the absence of gravity and aerodynamic drag. 
The effective exhaust velocity varies during burning, so th» 
familiar expression for the burnt velocity in terms of th: 
logarithm of the mass ratio is not valid here. It is necessar: 
to go back to first principles to express the burnt velocity 


dm 
vu = -f {1 
m 


The mass of the rocket at any time during burning is expresse« 
by 


1 
m= m+ Vf ode [2 


Specific impulse is an experimentally measurable propellant 
characteristic which we choose to relate to the effective ex- 
haust velocity and the effective specific enthalpy in the 
following way 


ur = V2JH = [3] 


In dealing with mixtures of two propellants, let subscript 
A identify the more dense propellant while subscript B desig- 
nates the less dense propellant. The density and effective 
specific enthalpy of a mixture are assumed to be volumetri- 
cally related to the densities and effective specific enthalpies 
of the separate propellants by 


p = pay + pa(l — y) [4] 
pH = paHay + psHa(1 — y) [5] 


in which y represents the volume fraction of propellant A. 
It is convenient to employ the following dimensionless 
parameters 


(6) 
pe 
paHa 
b = a | 7 
[7] 
my 
= — 8 
c [8] 


Propeliant B 


Propellant A 


Fig. 1 Schematic representation of a rocket with a variable 
mixture of two propellants burning simultaneously in the same 
chamber 
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Fig. 2 Qualitative behavior of the function ,(x) 


By means of a combination of the foregoing relationships, 
Equation [1] can be expressed in the form 


1 + + by) 
dx 


1 
6+ + ay)dx 


The problem is to find the function y(x) that maximizes the 
integral in Equation [9] and to find the maximum value of 
vy. 

The four cases to be considered are represented in Fig. 2. 
They are distinguished by various combinations of boundary 
conditions involving one or both of the propellants at the be- 
ginning and the end of burning. Cases I and III were treated 


in (1). 


Results 
Case I 
a+b 
Vi = [10] 
(a+ble-a+b 
“ala — b — 2be) 
yr) = be) 
(a — b)*x? — 4be(a — b — be) [12] 
Case II 
+b + 
a—b a(l + a) (1 + 6) 
(a+ble-a+t+b 
a(a — b — 2bc) 
When <x <1 
(a — 


2ab 2ab*/(a — b)*x? — Abe(a — b — be) 
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[9] 


X 
[17] 
a—b—be b 
ac re 2a(1 + 6) 
Case III 
a+b 
Yi > ab [18] 
(a + b) (1 +) 
[19] 
When 0 < x < x 
(2) a+b (a — b)x (20) 
2ab 2bV/ atx? — +c)? 
(a — b) (1 +e) 
y= + In [21] 
Case IV 
(a + b + 2ab\(1l + ¢) [29] 
2 
2aV(1 + + 6) 
23] 
When x, < < x2 
2ab 2bV — ab(l +c)? 
vy = + In 
2ac [25] 


(: + In + 
Discussion 


Which of the four cases one has to deal with in any particu- 
lar situation is most easily determined by evaluating the 
quantities y;, ys, % and x2. Only those with values between 
zero and unity have physical significance. This requirement 
generally eliminates all but one case. There are extreme 
situations in which none of the y’s or 2’s have physically per- 
missible values. In these situations, only the propellant with 
the greater impulse per unit volume should be used. 

It is of interest to note that the expression for y; represents 
the mixture of two propellants with the greatest impulse per 
unit volume. In all cases, the propellant or mixture of pro- 
pellants with maximum impulse per unit volume is used at 
the start of burning. 

Tt can be shown that the effective exhaust velocity increases 
at exactly the same rate as the velocity of the rocket during 
those intervals when both propellants are being burned. 
This is a curious result that is no doubt closely related to a 
similar result obtained by an energy optimization process 
treated in an earlier paper (2). 

The possible benefits from practical applications of results 
reported here become significant when one admits the use of 
propellants with greatly different densities. The feasibility 
of doubling the density of a solid propellant by addition of 
large amounts of heavy material, e.g., lead oxide, while re- 
taining good physical properties and combustibility has been 
demonstrated (3). Even though the resulting propellant has 
a rather poor specific impulse, it can be shown that combina- 
tions of the dense propellant with the initial undoctored pro- 
pellant have propulsive capabilities greater than either pro- 
pellant alone. 

The analysis that yielded the results reported here has been 
extended to include the second stage of a multistage rocket 
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(4). It appears that the range of useful applicability of vari- 
able propellant properties is much more restricted in the 
second stage than it is in the first stage. 


Nomenclature 


a, b, c, = dimensionless parameters defined by Equations {6, 
7 and 8], respectively 

Jo = factor for converting specific impulse to velocity, 
32.174 (fps)? 

H = effective specific enthalpy, Btu/Ib 

Tey = specific impulse, sec 

J = factor for converting specific enthalpy to velocity, 

2.503 10‘ (fps)? lb/Btu 

rocket mass 

exhaust velocity, fps 

rocket velocity, fps 

initial propellant volume, ft* 

volume fraction of propellant consumed 

volume fraction of propellant A in burning mixture 

= propellant density, lb/ft* 


Subscripts 


the more dense propellant 

the less dense propellant 

final value 

initial value 

beginning of variable propellant mixture for cases II and 


end of variable propellant mixture for cases III and IV 


to 
ll 
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Augmentation of Nuclear Rocket 
Specific Impulse Through 
Mechanical-Electrical Means 


MARTIN GOLDSMITH! 
Rand Corp., Santa Monica, Calif. 


N A RECENT paper by Ackeret (1)? it is pointed out that 

a conventional chemical or nuclear rocket operating in 
space has available to it an essentially infinite expansion 
ratio. Theoretically, therefore, one could obtain complete 
conversion of the thermal energy of the hot propellent gas 
into directed kinetic energy. This would, of course, require 
a nozzle of infinite area ratio, hence this ideal situation is pre- 
cluded by practical considerations. 

Ackeret further points out that isothermal expansion of the 
rocket exhaust leads to improved performance as compared 
to adiabatic expansion for a given expansion ratio. This, of 
course, requires the addition of heat to the gas during the 
expansion. This result is difficult to achieve in a practical 
system, especially if one is limited by permissible material 
temperatures. 

However, Ackeret suggests a stratagem by which the avail- 
ability of the large permissible expansion ratio is exploited 
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to obtain a quasi-isothermal expansion of the exhaust gas 
in a series of turbines. The gases are assumed to be reheated 
by a temperature-limited nuclear reactor. The turbine work 
is then converted into electricity, which is used directly to 
heat the gases further after their last heating in the reactor 
and before they are expanded in the exhaust nozzle. In this 
way, the temperature limitations imposed by the reactor 
materials are partially circumvented. 

To clarify this proposed cycle, the sequence of events is 
restated. The propellant, perhaps hydrogen or helium, is 
passed through a reactor or a part of a reactor at some initia! 
pressure. The heated gas is then expanded through a turbine 
and is then returned to a reactor at the lower exhaust pressur:: 
for reheating to the original limiting reactor temperature 
This cycle may be repeated many times, mechanical com. 
plexity and size being the limiting criteria as the pressure 0’ 
the gas is reduced through continued expansion. 

The turbine work is converted into electricity by conven- 
tional means. After the gas has made its last pass throug] 
the reactor, it may then be heated to a temperature abov: 
that possible in the reactor by discharge of the electricit, 
through it. The gas is then expanded through a thrust nozzl 
in the conventional fashion. 

The purpose of this note is to calculate the augmentatio: 
of specific impulse possible through use of Ackeret’s cyck 
for some realistic designs. Values will be assumed for th« 
various parameters that enter into the problem. These in- 
clude: Number of reheat cycles, turbine inlet temperature 
initial pressure, turbine pressure ratio, turbine adiabatic 
efficiency, maximum reactor gas temperature, nozzle pressur: 
ratio. Obviously, covering a wide range for each variable 
in the calculations would lead to a prohibitive number of 
cases; hence typical single values will be chosen for most of 
them. The use of hydrogen as a working fluid will be as- 
sumed. The appropriate equilibrium value of dissociation 
is assumed to occur through the cycle. 

For example, the thrust nozzle expansion ratio is assumed 
to be 1000:1 in all cases. This corresponds to a nozzle area 
ratio of roughly 40:1. Each turbine is taken to operate at 
70 per cent of isentropic efficiency with a pressure ratio of 
10:1. Maximum permissible turbine inlet temperature is 
taken to be 1500 K. The turbine output is considered to be 
entirely converted into electricity. Pumping and flow losses 
are neglected. 

Specific impulse is calculated for a range of reactor exit 
temperatures from 2500 to 3500 K. Two initial pressures to 
the nozzle are considered, 1 atm and ;'5 atm. The perform- 
ance for one and two reheat cycles, as well as the conventional 
case of no reheat cycling is shown in Fig. 1 for the listed condi- 
tions. Sample calculations are shown below for two typical 
cases, using the thermodynamic properties of hydrogen given 
in (2). 


Case 1 


Initial nozzle pressure = 1 atm 

Reactor exit temperature = 3500 K 

No reheat cycle: 
h (enthalpy) at 1 atm and 3500 K = 29.4 keal/gm 
h after adiabatic expansion to 10-* atm = 8.6 
enthalpy difference (Ah) = 20.8 
I.p = (1/g) V2 X 4.186 X Ah = 1350 sec 


Case 2 


Initial nozzle pressure = 0.1 atm 
Reactor exit temperature = 3000 K 
Two reheat cycles: 
h at 10 atm, 1500 K = 5.2 keal/gm 
h after isentropic turbine expansion to 1 atm = 2.6 
isentropic turbine work = 2.6 
heating energy available = 0.70 X isentropic turbine 
work = 1.8 
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Fig. 1 Specific impulse performance of the Ackeret cycle 


h at 1 atm, 1500 K, after reheat = 5.2 

A after isentropic turbine expansion to 0.1 atm = 2.8 

isentropic turbine work = 2.4 

heating energy available = 0.70 X isentropic turbine 
work = 1.7 

h at 0.1 atm, 3000 K = 24.9 

total electrical heating energy = 3.5 

total h at 0.1 atm before nozzle = 28.4 

h after nozzle expansion to 10~4 atm = 9.0 

enthalpy difference (Ah) = 19.4 

I.p = (1/9) V2 X 4.186 X Ah = 1300 sec 


Certain trends may be observed in the data shown in Fig. 1. 
The reheat cycle is of most benefit at lower reactor temper- 
atures. The large energies absorbed by dissociation at higher 
temperatures reduce the effect of the electric afterheating. 
Likewise, at the lower nozzle inlet pressure, the consequent 
increased level of dissociation further accentuates this trend. 

Therefore, tentative conclusions may be drawn as to the 
value of the reheat cycle in increasing engine performance. If 
the permissible reactor temperature is quite high, only a 
marginal increase in performance is available, unless many 
reheat stages are used. Mechanically this is probably not 
feasible. 

At the lower end of the temperature spectrum, it is seen 
that even two stages of reheat offer a worthwhile increase in 
specific impulse. Only a complete analysis of the perform- 
ance of the engine and vehicle will indicate whether overall 
performance is increased and, if so, whether to a degree com- 
mensurate with the necessary increase in engine complexity. 

It is interesting to compare these results to the performance 
obtained through isothermal expansion of hydrogen gas in a 
nozzle, where heat is continually added to the flowing gas. 
For an infinite pressure ratio in isothermal expansion, the 
specific impulse, according to simplified analysis, becomes 
infinite. However, in Table 1, the realistic case of a 1000:1 
pressure ratio, starting at a nozzle inlet pressure of 1 atm, is 
shown for several temperature levels. 
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Table 1 Specific impulse for isothermal expansion of 


hydrogen 
Temp, K Tsp, sec 
2500 1270 
3000 1610 
3500 1910 


A substantial amount of the heat added during isothermal 
expansion is consumed in the dissociation energy of the gas. 
Unlike the case of adiabatic expansion, this energy is not 
useful in the expansion process and represents dead loss. 
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Beamed Electromagnetic Power as a 
Propulsion Energy Source 


MARTIN I. WILLINSKI' 


Hughes Aircraft Co., Culver City, Calif. 


A novel rocket propulsion system, utilizing an external 
energy source in the form of beamed electromagnetic 
power, is proposed. A comparison is made with solar 
energy. Its application to propulsion fluid heating and 
electric rockets, as well as to the conversion of a portion of 
the beamed energy into propellant mass, is discussed. 
The problems of electromagnetic power beaming are re- 
viewed, and some potential solutions discussed. The 
possibility of utilizing iconospheres as a means of beaming 
power over astronomical distances is suggested. 


IHE PURPOSE of this paper is to suggest and discuss a 

novel rocket propulsion system which might permit the 
achievement of currently difficult objectives of space travel. 
Although this discussion is speculative in nature, it is believed 
appropriate in view of anticipated technological progress and 
the long lead time required between the conception and real- 
ization of advanced systems. 

Two basic limitations to the achievement of many of the 
future objectives of space travel are: 

1 The limit of the total quantity of energy which can be 
stored in a space vehicle, whether it be nuclear or chemical, for 
eventual conversion into power in the vehicle’s reaction jet. 

2 The limit of the vehicle specific power which may be ex- 
pressed as power in reaction jet divided by vehicle weight 
(1,2).? 

For example, the magnitude of both of these propulsion 
system characteristics must be extremely high for vehicles 
capable of short duration interplanetary (3) and interstellar 
flights. Even the most energetic power sources, such as the 
fission and fusion fuels, have been shown to be limited in these 
respects (4). 

A method for circumventing these limitations is the utiliza- 
tion of a rocket propulsion energy source external to the ve- 
hicle. This bas been suggested previously only in the form of 
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natural energy sources, such as solar energy (5, 6). Although 
this utilization of solar energy would overcome the limitation 
of total required energy, it does not appear capable of over- 
coming the specific power limitation. This is a result of the 
constant value of intercepted power per unit solar energy 
collector area at any given distance from the sun. 

A different type of external propulsion energy source is 
proposed here to overcome the specific power as well as total 
energy limitations previously discussed. This is an artificial 
energy source in the form of beamed electromagnetic power. 
For a beamed power energy source, the intercepted power 
may be increased as long as the transmitted power is in- 
creased. A comparison of beamed and solar energy on an 
intercepted power basis is shown in Fig. 1. It should be noted 
that the solar energy regime would vary in area with varia- 
tions of distance from the sun. 


Potential Utility 


Beamed power is seen as being utilized as a rocket propul- 
sion energy source in two ways: 

1 The beamed power would be utilized to heat a propulsion 
fluid, such as hydrogen, ammonia, gasoline or even water, 
which would then be expanded through a nozzle to produce 
thrust. 

2 The beamed energy would be converted into electrical 
power suitable for powering an ion or plasma rocket. It is 
noteworthy that the use of beamed electromagnetic power as 
a propulsion energy source for an airborne vebicle powered by 
an aeropulse-type engine has been proposed (7). 

The first method appears to offer the performance ad- 
vantages of a nuclear propulsion fluid heating type rocket with 
a considerable reduction in onboard shielding problems. In 
addition, any method of utilizing a beamed power propulsion 
energy source has the advantages of a fixed primary power 
installation. 

The receiving antenna could take the form of a large in- 
flated balloon, in which case drag considerations would limit 
such a system to essentially vacuum conditions and hence 
vehicle upper-stage operation. Also, the skin of the vehicle 
itself could possibly be utilized as a surface antenna, thereby 
allowing operation in the atmosphere. Beam transmission 
ranges up to several hundred miles, and 500 to 1000 mega- 
watts of power received and utilized for propulsion would be 
suitable for an upper stage of a satellite booster using such a 
system. 

The application of beamed power to ion rocket propelled 
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Fig. 1 Solar and beamed power regimes 


vehicles may offer a means of attaining vehicle accelerations 
orders of magnitude greater than those which have been 
proposed in past design studies of such vehicles (2, 8, 9). One 
such vehicle might be a high velocity “artificial comet’’ for 
reconnaissance in the vicinity of the outer planets of the solar 
system. For this application, beam transmission ranges of up 
to 100,000 miles with approximately 10 megawatts of power 
received and utilized for propulsion would be required. In 
this case, the power would be beamed from multiple stations 
during the vehicle’s single powered maneuver in the vicinity 
of the Earth. 

The beam ranges and power levels required for minimum 
duration interplanetary trajectories and interstellar flight 
trajectories could possibly be orders of magnitude greater than 
that required for an “artificial comet.’”’ A use of beamed 
power for the more distant future might involve the conver- 
sion of some fraction of the intercepted beamed energy into 
propellant mass which would be accelerated from ion motors 
powered by the remaining energy intercepted from the beam. 


Problem Areas 


The current engineering limitations on the problem of 
beaming large amounts of electromagnetic energy appear most 
formidable. This becomes apparent from the basic beamed 
power equation® 


[1] 


Ideal far-zone antenna gain, as used in Equation [1], is given 
by 


G = [2] 


It has been shown that far-zone gain may be retained under 
certain circumstances in near-zone operation of antennas by 
application of a focusing principle (10). It has also been 
shown that there is an optimum size for an antenna operated 
under otherwise fixed conditions (11). 

As can be seen from Equations [1 and 2], the quantity of 
power received varies inversely as the square of the antenna 
separation distance. Some quantitative indication of possible 
resulting power efficiencies is shown in Fig. 2 where the ratio 
KP,/Pr is plotted as a function of the antenna separation 
distance for a range of values of transmitter antenna gain and 
receiving antenna diameter. 

Although the antenna gains shown in Fig. 2 are not repre- 
sentative of what has been achieved, they do represent values 
which are theoretically attainable under certain conditions. 
There are many problems which will be encountered in the 
development of antennas having such high gain. For ex- 
ample, the actual gain of an antenna is dependent on the an- 
tenna fabrication tolerances. For maximum gain, a para- 
baloid dish antenna reflector should be constructed with a 
tolerance of approximately ;*, of the operating wave length. 
This tolerance might be difficult to achieve for very large size 
antennas. The maximum gain is also affected by the en- 
vironment in which the antenna is operating; higher gains 
should be attainable in space than within the Earth’s at- 
mosphere. 

The desirability of lower wave length is apparent from 
Equation |2]. However, there are factors which tend to limit 
the maximum usable transmitting frequency. For example, 
the maximum power obtainable from any particlular type of 
known power transmitting tube decreases with increasing fre- 
quency (12). Also, power attenuation by the water vapor and 
oxygen in the atmosphere generally increases with decreasing 
wave length starting at a wave length of approximately 10 
em. 


’ This equation is derived from geometrical optics, and there- 
fore the ratio Pr/P7r will, in practice, be reduced by the antenna 
aperture efficiency factor K which depends upon several of the 
design parameters of the antenna system. 
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Fig. 2 Beamed power reduction with distance 


An additional problem area is the voltage breakdown of the 
mospheric constituents. This phenomenon would limit the 
aximum power density of the beam. Completely new ap- 
coaches may be required for the solution to some of these 
id other problems associated with beaming electromagnetic 
wer. 


(Cutlook 


There are several prospects for obtaining solutions to these 
cifficult problems. Current and suggested future studies con- 
cerning space communications may result in techniques for 
constructing larger high gain antennas in space than may be 
possible on Earth. This, in turn, suggests the possibility of 
‘frequency switching,’ in which one frequency, selected with 
consideration of atmospheric attenuation, would be used for a 
“shorter” surface-to-space antenna link and another higher 
frequency for the “longer” space antenna-to-space vehicle 
power transmission link. 

If enormous auantities of cheap electrical power become 
available for propulsion purposes, extremely low beam power 
efficiencies would become usable. Thus the ratio of received- 
to-transmitted power may be reduced to an extremely small 
quantity, thereby permitting operation at extreme ranges. 
The eventual harnessing of fusion energy provides a hopeful 
solution to the problem of providing enormous quantities of 
electrical power for propulsion purposes. It is believed that 
when controlled fusion power becomes a reality, fusion power- 
plants will produce considerably more electrica] power than is 
produced by conventional or nuclear fission-type powerplants 
(13). 

One approach to the problem of antennas for beaming 
power over astronomical distances might be the utilization of 
a natural electromagnetic radiation focusing device as sug- 
gested by Berger (14). This natural focusing device is termed 
an iconosphere by Dr. Berger and may be defined as a series of 
spherical, optically-real focal surfaces concentrically surround- 
ing a spherical body. Thus, the spherical body acts as a posi- 
tive lens on incident electromagnetic radiation by refraction, 
diffraction and/or gravitation. He suggests the use of the 
iconosphere of the Earth, or even the sun, as a means for 
communicating over interstellar distances. The prospect of 
applying iconospherics to the problem of beaming power for 
propulsive energy purposes appears intriguing; indeed, icono- 
spheric lenses, combined with beamed power, might prove to 
be the key to the accomplishment of interstellar travel. 


Conclusion ‘ 

The realization of propulsion systems utilizing beamed power 
energy sources would represent a significant advance in the 
areas of propulsion and space travel. It is hoped that the con- 
cepts discussed here will stimulate research toward the solu- 
tion of the many problems associated with the successful de- 
velopment of beamed electromagnetic power as a propulsion 
energy source. 
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Nomenclature 

Ar = effective area of receiving antenna 
Ar = effective area of transmitting antenna 
d = distance between the transmitting and receiving antennas 
Dr = physical diameter of receiving antenna 
G = transmitting antenna gain 

K = antenna aperture efficiency factor 

Pr = power received 

Py = power transmitted 

\ = wave length of transmitted power 
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Utilization of Stored Energy in Solid 
Material as an Added Energy Source 


JOHN J. NEWGARD' 


Thiokol Chemical Corp., Reaction Motors Division, 
Denville, N. J. 


A’ ATTEMPT has been made to survey two types of 
solid-state phenomena in order to determine the ca- 
pacity inherent in these phenomena for energy storage (1).? 
The first process to be discussed is analogous to energy stor- 
age via an order-disorder transformation known to exist in 
some metallic allovs: Disordered material is ordered; energy 
is added and stored in the ordered structure; upon combus- 
tion the ordering energy as well as the combustion energy 
is released. The second mechanism is that of interstitial- 
vacancy formation. Interstitials and vacancies are created 
by heat and irradiation. Somehow these are retained until 
combustion. Upon combustion, the energy associated with 
the interstitials’ recombination with vacancies is released 
along with the normal combustion energy. 

Specific examples and orders of magnitude associated with 
both phenomena are given. The conclusions of this study 
can be given now. Even if an order-disorder or interstitial- 
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vacancy phenomenon existed in readily combustible mate- 
rials, the energy added to combustion is small (approximately 
30 per cent) under the most reasonable conditions assumed 
in the addition of stored energy to these materials. 

If one assumes the existence of a perfectly ordered com- 
pound AB, it can be shown that the ratio under ideal circum- 
stances (that the interaction energies of any postulated lattice 
configuration for nearest neighbor atoms in ordered alloy 
Vaz is much greater than the interaction energy between 
nearest neighbor atoms V 44 and V gz in the disordered alloy), 
the ratio of ordering energy to internal energy before ordering 
is equal to unity. Internal energy before ordering may be 
expressed as 


Eup Er 
(AB )aisordered solid (A point (AB) 


where 


Exp = energy required to reach melting point of AB 
Er = energy required to melt AB (fusion energy) 
Uo = Eup + Er 


where now U, is the internal energy of the solid in question. 
It can also be shown that the ordering energy of AB is greater 
than ordering energies associated with AB, where n is an 
integer. 

If one assumes some order of magnitude parameters, such 
as average melting points of 100 to 800 C, heat capacities of 
solid material from 0.1 to 1 cal/gm, heats of fusion of 10 to 100 
cal/gm and average molecular weights of substances from 
10 to 400, ordering energies can vary from 2 X 10? to 1.2 x 105 
cal/mole. 

Combustion phenomena give heats of combustion which 
vary from 2 X 104 to 2 X 10? cal/mole with most substances 
giving energies upon combustion in excess of 105 cal/mole. 
Hence, the upper order of magnitude of ordering energy and 
the lower value for heat of combustion of chemical substances 
have the same magnitude. 

Because of the assumption of ideal circumstances in the 
order-disorder phenomenon, it has been stated that ordering 
energy is of the same magnitude as internal energy. Based 
upon more realistic assumptions that Vas = Ves = 3 Vaz, 
ordering energy now becomes one third of the internal energy, 
and hence, if an order-disorder phenomenon can exist, at most 
a 30 per cent increase over the ordinary combustion energy 
can be expected in those cases where the combustion energy 
is as low as 10° cal/mole. 

Order-disorder phenomena have appeared only in metallic 
alloys. If a useful propellant alloy, such as MgAl (magnal- 
ium) would exhibit order-disorder (to the author’s knowledge 
this alloy hasn’t shown the order-disorder phenomenon), the 
additional energy added to the combustion energy would 
only be about 15 per cent, since the heat released on ‘“burn- 
ing” MgAl is greater than 2 X 105 cal/mole. In conclusion 
then, order-disorder energy addition to combustion energy 
makes only a small addition to the useful energy attainable 
from a combustion phenomenon. 

Two other stored energy phenomena are those which store 
energy by Schottky and Frenkel defects. Both these mech- 
anisms store energy by movements of lattice atoms. If 
the lattice atoms move to the surface, creating a surface 
atom and a vacancy, a so-called Schottky defect is created. 
Movement of lattice atoms into interstitial positions in the 
lattice, creating interstitial atoms and vacancies, is called a 
Frenkel defect. In both cases this stored energy is assumed 
released upon combustion. 

It has been shown (1) that for the creation of a Schottky 
defect approximately 1 ev is required (or stored), and approxi- 
mately 2 ev are required to create a Frenkel defect. However, 
the number of defects created in either case is proportional 
to e~*'*", where E is the energy required to create the defect, 
and T is the temperature. It can be shown for both Schottky 
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and Frenkel defects that, in order to increase the energy upon 
combustion by 8000 cal/mole, approximately one third of all 
lattice atoms must appear on the surface or in interstitial 
positions (assuming no recombination, i.e., surface atoms or 
interstitial atoms returning to their lattice positions), and the 
temperature at which this occurs is 11,600 K. At ordinary 
temperatures (300 K), only one atom in 10° is either a Schottky 
or a Frenkel defect. Hence, these phenomena, even if they 
did exist in useful combustible material at reasonable tem- 
peratures, would represent a minute fraction of the energy 
available from ordinary combustion. 

A special case may be made for neutron radiation-induced 
energy storage in graphite (2). High energy neutrons in 
passing through a graphite lattice elastically transfer ap- 
proximately 1/1000 of their initial energy to carbon atoms. 
These subsequently move to higher energy interstitial posi- 
tions (creating vacancies) not unlike Frenkel defects. By 
this means energy storage as high as 540 cal/gm at 30 C has 
been obtained experimentally, using reactor neutrons at 
integrated fluxes of 2 * 10?! nvt (reactor flux times time in 
units of 1/em?), 

Graphite is the only substance showing such a large 
energy storage effect at ordinary temperatures without the 
effects annealing out. A semi-empirical model (2) predicts 
saturation in energy storage between 600 and 700 cal/gm. 
For normal graphite oxidizing in air at 25 C, the reaction is 


C + CO. 


and the heat of reaction is 2140 cal/gm. Hence, the max- 
imum increase to be expected in ‘combustion’ energy at 
25 C if graphite is to store energy due to neutron irradiation 
is about 30 per cent. 


References 


1 Kittel, C., ‘Introduction to Solid State Physies,’’ John Wiley and 
Sons, Ine., New York, 1956. 
2 Newgard, J. J., “Simple Semi-Empirical Model for Neutron Induced 


Stored Energy in Graphite,” J. Appl. Phys., in publication. 


Optimum Magnetic Field for Stagnation 
Heat Transfer Reduction at Hypersonic 
Velocities 
ROBERT GOULARD! 

Bendix Products Division, South Bend, Ind. 


The application of a magnetic field to the stagnation 
area of a blunt body has two contrary effects: It reduces 
convective energy transfer, but increases radiation energy 
transfer. Consequently, for very high speeds, when this 
latter effect is existent but not dominant (q; < 0.4q,°’), 
an optimum magnetic field exists which insures minimum 
energy transfer. At still higher speeds, radiation is 
dominant, and the application of a magnetic field always 
increases the energy transfer to the stagnation area or the 
nose. This argument closely parallels the well-known 
argument on the existence of an optimum radius. The 
existence of an optimum magnetic field is shown, in fact, 
to depend on the choice of a nose radius below its optimum 
value. 


a* SHOWN in recent papers, the application of a mag- 
netic field at the blunt nose of a hypersonic vehicle re- 
duces the convective heat transfer to the surface of the nose. 
Received April 21, 1959. 
' Consultant, Applied Sciences Laboratory; also Assistant 
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This is due essentially to the retarding effect imposed on the 
particles streaming past the wall behind the detached shock, 
resulting also in an increased shock layer thickness. 

At high temperature, however, radiative flux becomes an 
appreciable part of the total energy flux, and its value in- 
creases proportionately to the shock layer thickness (3).? 
The presence of a magnetic field will therefore have two op- 
posite effects on the two components of the total energy flux. 
As a result, a minimum energy flux will usually be achieved 
for a certain optimum value of the magnetic field. Higher 
values will, in fact, cause an increase of the energy flux at the 
wall, 

These two opposite effects have been analyzed in (1), where 
they are illustrated in two cases of hypervelocity flight 
= 8and10km/sec). An order of magnitude discussion, 
sugested by Miele, will be used here to present in simple 
form the main conclusions of (1). 


Optimum Magnetic Field 


't has been shown in (2) that the application of a magnetic 
fie'd at the stagnation point of a blunt body reduces the con- 
ve: tive heat transfer g. according to the expression 


= qe [1] 


where Z is a magnetic parameter. We note that Z, and there- 
fore ge, decreases when the applied magnetic field B increases. 
Similarly, it is shown that the velocity gradient varies in 
proportion to Z. Hence, by continuity, the shock layer 
thickness is roughly inversely proportional to Z, and so is the 
energy radiated by the shock layer 


One recognizes the opposite effects of the magnetic param- 
eter Z on the two components of the total heat transfer q. 
A minimum value of q is obtained by the condition dq/dZ = 0, 
which is satisfied by the optimum value of Z 


gp 


with 0<Z<1. 

Hence, whenever radiation heat transfer is present, we con- 
clude to the existence of an optimum value of the applied mag- 
netic field corresponding to a minimum total heat transfer. 
For instance, it is shown in (1) that for u.. = 26,200 fps, 
altitude = 120,000 ft, R = 1 ft, the optimum value of B 
will be of the order of 3 X 104 gauss ({=4), corresponding to 
a 20 per cent reduction of the total energy transfer to the 
hemisphere nose. 

We also note that, by definition, Z < 1. It follows from 
Equation [3] that for flight configurations such that 


qr /qe > 0.4 [4] 


there will be no such optimum; in this case, the application 
of a magnetic field will actually increase the total heat transfer 
to the nose. 


Optimum Nose Radius 
The values of q. and gg, on the other hand, are known 


2 Numbers in parentheses indicate References at end of paper. 
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(3) to be of the form 


1/2 
qe & a [5] 


~ co R [6] 


By inspection, it can be concluded that condition gp2/q. < 
0.4 will not be maintained for “black” surfaces at very high 
velocities, large radii or low altitudes. A magnetic field, in 
this case, would not reduce the total heat transfer. This 
situation is shown in (1 and 8) to be typical of the direct re- 
entry of circumlunar vehicles (~10 km/sec, gp ™ q,). 

It is interesting to note, in Equations [5 and 6], that for a 
given aerodynamic situation, the nose radius R has also an 
opposite effect on the two terms of the total heat transfer 
gq. This leads to the well-known argument on the possible 
existence of an optimum radius, which exists for 0g/OR = 0, 
corresponding to 


qr /qe = 0.5 [7] 
Conclusion 


Comparison of Equations [4 and 7] shows that the applica- 
tion of a magnetic field at the nose will be conducive of energy 
transfer reduction only if the nose radius is less than a value 
inferior to its optimum value. In other words, if a decrease 
of q is to correspond to an increased magnetic field B, the 
radius of the nose must be such that further thickening of the 
shock layer will bring a decrease of the energy flux. 

This conclusion illustrates, in particular, the common 
effect of radius increase and magnetic field application, which 
is the thickening of the shock layer. In terms of energy 
transfer, the presence of a magnetic field can be considered 
as equivalent to the fictitious substitution of a larger nose 
radius. This magnetic field can be adjusted to the various 
optimum requirements of a given flight. The overall effect 
on the drag and the trajectory is briefly discussed in (1). 


Nomenclature 


constants 
magnetic field 


Qc, ar 


Ger Ve convective heat flux, with and without magnetic field 
qr; Gr” net radiation heat flux, with and without magnetic 
field 
q = total energy flux, g. + qr 
R = nose radius 
u = velocity 
x = distance from the stagnation point along the nose 
Z = + (1 + 
Ay = wall absorptivity coefficient 
\ 0x /B=0 
p = density 
o = electrical conductivity 
Subscripts 
es = stagnation 
co = upstream of shock 
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Technical Comments 


Evaluation of Conical Nozzle Thrust 
Coefficient 


G. V. R. RAO! 


Rocketdyne, Canoga Park, Calif. 


HE THRUST of a nozzle based on one-dimensional flow 
of ideal gas is given by 


F; (w/g)V, + (DP, {1] 


This formula is derived under the assumption that at the exit 
section A, we have uniform flow of velocity V, parallel to 
nozzle axis. The pressure of the gas in the exit plane and the 
ambient pressure are denoted by p, and p,, respectively. In 
(1)? it is shown that the exit flow of a conical nozzle can be 
represented by a source-type flow, i.e., across the spherical 
area A,, as shown in Fig. 1, one can assume a constant 
velocity V, and flow direction normal to the spherical surface. 

Upon integration, one obtains thrust of a conical nozzle as 


Fonical = A(w/g)V, + (Pp, [2] 
where 
1+ cosa 


In this equation, V, and p, are the velocity and pressure cor- 
responding to expansion of the gases to area A,, given by 


A, = (1/A)A, 


The use of Equation [2] was clarified in (2), and the effect of 
using incorrect area A, was discussed. One obtains different 
values for exit momentum depending upon whether the ex- 
haust gases starting from the throat section are assumed to 
expand to exit area A, or A,. However, the importance of 
the critical flow area, which is generally less than the geo- 
metric throat area, in determining V, and p, was not dis- 
cussed in (2). 

The influence of the nozzle wall at the throat causes non- 
uniform flow conditions in the throat region. The conse- 
quent sonic line configuration is treated in (3). In the lin- 
earized treatment of the transonic flow, the shape of the sonic 
line turns out to be parabolic and depends upon the ratio of 
throat radius to the radius of curvature of the nozzle wall. 
One can compute the mass flow crossing the throat section 
and find that it is less than the critical flow corresponding to 
uniform flow conditions across the plane area A, Let A* 
denote the critical area corresponding to the mass flow w/g 
through the nozzle 


R 
= p*V*A* = pV cos 62rRdR = Dp*V*A, [3] 


The coefficient D in Equation [3] takes into account the loss 
in theoretical mass flow obtained for ideal gases, because 


Received June 11, 1959. 
1 Engineering Specialist. Member ARS. 
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of the influence of nozzle walls. This coefficient should be 
distinguished from the discharge coefficient of a rocket nozzie, 
caused by nonideal properties of gases in addition to the 
wall effect. 

The flow coefficient D is always less than unity. Assuming 
y = 1.23, this coefficient is computed for various values of 
nozzle wall curvature and is shown in Fig. 2. For the sake of 
comparison, the values of D for y = 1.4 are also shown in te 
figure. 

Since critical area is less than throat plane area, tlie 
theoretical estimate of mass flow can be modified as 


+ 


where p, and 7’, denote the chamber pressure and temper:- 
ture, respectively. + is the isentropic exponent, and R is tlie 
gas constant of the exhaust gases. 

Rearranging terms in Equation [2] one obtains 


or in terms of thrust coefficient 


In the above equation C’,,,, represents one-dimensional thrust 
coefficient of a nozzle discharging into a vacuum at an area 
expansion = A;/A*, 

The following formula is widely used in computing the 
thrust of a conical nozzle 
CP p,A,) = P.)A, [6] 
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Fig. 1 Conical nozzle parameters 
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Fig. 2 Variation of flow coefficient with wall curvature 
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where Cr,,, represents one-dimensional vacuum thrust for area 
expansion e = A,/A, To evaluate the difference between 
Equations [5 and 6], one observes the following relation de- 
rived from the dependence of one-dimensional thrust co- 
efficient on area expansion 
where M, is the Mach number corresponding to area ratio e. 
From the definition of vacuum thrust 


= (pe/p-Je(1 7M,’) [8] 
sul stituting Equations [7 and 8] into Equation [5] gives 


e’ and € can be expressed in terms of A, and A, in this equation, 
which reduces to 


F conical = CP vac + AD Pal Pe) [9] 


The thrust computation of a conical nozzle as given in 
Equation [9] is more accurate than using [6], under the as- 
suinptions involved, since Equation [9] includes the effect of 
nozzle walls in the throat region. Ignoring this effect means 
assuming critical flow through throat area A, Then the 
factor D becomes unity, and the above equation reduces to 
Equation [6]. In the first approximation, we will have 

F conical D-F, 


conical 


the effect being mainly due to change in the critical mass flow 
as given by Equation [4]. 
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Stabilized Gaseous Detonation Waves! 


J. A. NICHOLLS? 
The University of Michigan, Ann Arbor, Mich. 


ECENT years have witnessed the steady increase of in- 
terest in the area of detonative combustion and, of par- 
ticular interest here, the concept of stationary detonation 
waves. Recently, Gross (1,2)* reported on the attainment of 
strong stabilized detonation waves and made reference to 
some of our work. In view of his comments, I feel obligated 
to elaborate on our experimental procedure and results, to 
question some of the conclusions he has drawn, and to point 
out that a standing detonation wave (possibly of the Chap- 
man-Jouguet type) was first achieved at the University of 
Michigan Aircraft Propulsion Laboratory in May 1958. 
Personnel of this laboratory have been studying this 
problem under sponsorship of the Air Force Office of Scientific 
Research since 1954. This work has resulted in an experi- 


Received May 28, 1959. 
! This research was supported by the U.S. Air Force under 
contract AF 18(600)-1199 monitored by the Air Force Office of 
Scientific Research, Combustion Dynamics Division, Air Re- 
search and Development Command. The work is being contin- 
ued under contract AF 49(638)-562. 

2 Research Engineer, The University of Michigan Research 
Institute. Member ARS. 
* Numbers in parentheses indicate References at end of paper. 


Aueust 1959 


mental facility for generating detonation waves, and much ex- 
perimental information has been obtained. The essential 
features of the experimental arrangement are shown in Fig. 1. 
The nozzle is operated highly underexpanded and discharges 
to the atmosphere. Typical operating conditions are as fol- 
lows: 

Stagnation pressure of the air and hydrogen, 350-765 psia. 

Stagnation temperature of the air, 1950-2600 R. 

Stagnation temperature of the hydrogen, 550 R. 

Mach number into the normal shock, 3.5-6.2. 

Static pressure behind the shock, 20 psia. 

Equivalence ratios, 0.33-1.5. 

Premature combustion of the gases in the nozzle or up- 
stream of the shock is ordinarily precluded by rapid expansion 
of the gases. Many runs have been made, and as indicated in 
the figure, a combustion zone is established which is normally 
located somewhat downstream of the shock front. This con- 
figuration is very steady and reproducible. The attainment of 
this phenomenon was first reported at the 2nd AFOSR Air 
Breathing Combustion Symposium, held at Fairchild Engine 
Division, N. Y., in April 1958. A formal paper covering the 
results of these experiments was presented at the Seventh 
Combustion Symposium, London, England, August 1958 (3). 

The separation between shock and combustion mentioned 
above corresponds to the ignition time delay for the given 
pressure, temperature and mixture ratio existing downstream 
of the shock. This delay is to be expected, and recent 
theoretical work by the author has substantiated, reasonably 
well, the delay times observed. These results will be pub- 
lished in the near future. Now ignition time delay is an ex- 
ponential function of reciprocal temperature, so it is not sur- 
prising that in those cases in which partial premature com- 
bustion in the nozzle was experienced (3) the delay decreased 
to shorter times than could be resolved experimentally. In 
effect, the premature burning served to increase the stagna- 
tion temperature of the stream and, accordingly, the tempera- 
ture downstream of the shock. The delay was thus ma- 
terially decreased. Another effect of the premature burning is, 
of course, a change in the gas composition which tends to ob- 
scure the reaction under study and which may or may not 
expedite the initiation of the combustion reaction. 

The question naturally arises as to whether the configura- 
tion mentioned is to be properly classified as a detonation 
wave or merely a case of shock ignition. In many cases we 
have recorded wave patterns wherein the combustion zone has 
had little or no effect on shock wave position and hence the 
Mach number of propagation of the wave. If one were to 
attempt to classify such waves on the basis of the fuel-air 
ratio and Mach number of propagation (on a one-dimensional 
basis), it would be found that some waves corresponded to 
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strong detonations, some to Chapman-Jouguet, and some 
would yield no physical solution. It seems apparent that the 
assumption of a constant-area stream tube breaks down in 
cases of this type, so that it becomes meaningless to attempt 
any designation in the conventional manner. Furthermore, 
with no strong interaction between shock and combustion, it 
hardly seems plausible that these could rightfully be called 
detonations. In fact, we have preferred to consider such 
phenomena as shock ignition. On the other hand, we have 
observed cases where the onset of combustion has driven the 
shock upstream to a new lower Mach number of propagation. 
Such a wave was first observed in May 1958 in which the 
stagnation temperature of the air was 2600 R, and the stagna- 
tion pressure of each stream was 765 psia. The initiation of 
combustion caused the wave to move upstream from a Mach 
number of 6.1 to about 5.7 This latter Mach number is close 
to the Chapman-Jouguet Mach number of detonation under 
the given conditions This I would refer to as a detonation 
wave, but will not attempt to define it precisely according to 
strength in view of difficulties mentioned. 

Gross (2) has pointed to the fact that in our experiments 
light emission is realized for appreciable distances downstream 
of the wave. On this basis, he expresses some doubt as to 
whether this could be a detonation wave. Actually, the per- 


sistence of the light emission can be attributed to the fact 
that we inject sodium chloride into the air stream. Further- 
more, there will be gradual combustion of the outer portions 
of the flow that did not pass through the normal shock. Thus, 
the light emission is irrelevant for either the verification or 
refutation of the existence of detonation. 

As is evident, there is considerable ambiguity surrounding 
the appropriate description of stabilized shock combustion 
waves. Certainly, much of this ambiguity is introduced by 
the finite reaction length required, which subjects the combu-- 
tion process to externally imposed pressure fields. It is fel:, 
then, that only when the reaction length is small compared to 
other characteristic dimensions will a true classification of 
standing detonation waves be possible. This condition can 
be met at higher temperatures where the kinetics are su’- 
ficiently rapid to allow the wave to be treated as a discor- 
tinuity. 
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Identification of Inert Objects in Space 
LOUIS G. VARGO! 


Aeronutronic, A Diy. of Ford Motor Co., Newport Beach, 
Calif. 


ANMADE objects in space which have ceased to serve 
their original purpose are variously referred to as “junk”’ 
or “debris.” These objects are non-emanating and may be 
considered as lost in the sense that tracking is impossible or 
unwarranted, or both. There exists a very real possibility, 
however, that future observations will disclose such objects 
either by chance or in a search program designed to utilize the 
inert equipment or structure for other purposes. The period 
of time that the body remains unobserved may be long indeed. 
During this period the complex effects of perturbative forces 
may become of large magnitude, so that radical departures 
from the original Keplerian orbit will result. It appears im- 
portant then, to inquire if such objects have some type of 
“signature” by which they can be positively identified. 

Using the reflected radiation characteristics of the body, 
for example, the radar cross-sectional area, appears not too 
attractive for several reasons. We will not explore this further 
except to note that identification by this method has basic 
difficulties associated with limited resolution capability and 
the attendant problems of ambiguity. 

In 1896, Tisserand? addressed himself to a problem which 
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arose in the field of celestial mechanics similar to the one 
considered here. Cometary orbits, perturbed by a major 
planet, such as Jupiter, experience large variations in their 
heliocentric elements between successive observations near 
perihelion. Since comets are virtually indistinguishable by 
optical means, the astronomer proposed using the value of an 
invariant function of the elements, the Jacobi integral, as the 
identifying characteristic of the body. With certain assump- 
tions, the integral takes the simple form 


J =~+2-+/a(1 — e2) cosi = constant (1) 


a 


Here a, e and 7 are, respectively, the semimajor axis, eccen- 
tricity and inclination of the Keplerian motion relative to the 
parent body. The unit of distance is the mean distance of the 
disturbing body. One merely catalogues values of J for future 
reference with the knowledge that any variation in J will be 
orders of magnitude down from changes in the elements them- 
selves. 

Several types of missions are currently being considered 
which result in a situation similar to Tisserand’s problem. 
Among these are terrestrial satellites with perigee distances of 
several Earth radii and lunar probes that do not impact or 
escape. These missions give rise to Earth-centered orbits 
which are dominantly perturbed by the moon. Also, the 
heliocentric orbit of a nonimpacting planetary probe might 
well be expected to undergo changes due to the target planet 
or the Earth on subsequent passages. To date, no artificial 
objects fall into these classes. One may reasonably anticipate 
the necessity of having identification methods available. 
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Dyna-Soar (for dynamic soaring) is a joint 
project between the Air Force and the NASA, 
and is an attempt to solve the technical prob- 
lems of manned flight in the sub-orbital regions. 
Advance knowledge on the project indicates 


- how a boost-glide vehicle can operate from the 


outer fringes of the atmosphere where it can 
maneuver and be recovered undamaged. Studies 
show that by varying the original rocket boost, 


Tmk(s) ® 


and thus the velocity, and with the control 
available to the pilot, the Dyna-Soar aircraft 
can circumnavigate the earth, followed by a 
normal and controlled landing. Boeing Airplane 
Company, one of the competing companies for 
the development contract for the complete 
boost-glide system, has delegated to RCA the 
responsibility for the development of important 
electronic components of Dyna-Soar. 
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New Patents 


Jet engine burner having means for 
reading the pilot flame (2,872,785). 
. E. Barrett and W. F. Payne, Cedar 

a N. J., assignors to Curtiss-Wright 


channel members extending 
from the pilot burner tube to annular 
channels having an open side facing 
downstream. Channels have a diameter 


larger than the tube. 
Thermoelectric cooling apparatus (2,872,- 
788). N. E. Lindenblad, Princeton, N. 


J., assignor to Radio Corporation of 
America. 

Thermocouple element for a heat pump. 
The element es a pair of bodies of dis- 
similar thermoelectric material mounted 
on a plate of conductive material. Mem- 
bers transfer heat between the element and 
the surrounding air. 

Rocket base fuze (2,872,869). V. K. 
Tasmussen, San Marino, Calif., assignor 
to the U.S. Navy. 

Deceleration discriminating fuze con- 
taining a fixed explosive train and a deto- 
nator plunger initially held against longi- 
tudinal movement. The fuze is detonated 
at a predetermined reduced deceleration 
by trigger spring releasing a firing pin. 
Igniter squib (2,872,870). W. A. Gey, 
China Lake, Calif., assignor to the U. 8. 

avy. 

Flame inhibiting means placed upon a 
cylindrical booster grain of propellant 
within a housing. The central bore of 
the grain contains compressed black 
powder pellets protected from deteriora- 
tion. The closure is destroyed by burning 
the igniter propellants. 

Solar energy apparatus (2,872,915). W. 
S. Bowen, Westfield, N. J. 

Transparent semicylindrical parabolic 
mirror forming a closure over a primary 
semicylindrical parabolic mirror. Both 
mirrors have the same axis of principal 
focus. Heat radiated from a heat receiv- 
ing tube is reflected from both mirrors. 
Flight control system (2,873,074). H. 
Harris, E. B. Hammond Jr., W. T. White 
and R. M. Giulianelli, North Bellmore, 
N. Y.. assignors to Sperry Rand Corp. 

On a dirigible craft, signals are caused 

by the displacement of the longitudinal 
axis from a reference direction, and by ac- 
celerations normal to the axis. Combined 
signals are used to control orientation of 
the craft. 
Remote control fluid pressure responsive 
actuator (2,873,648). C. W. Musser and 
F. W. Dietsch, Philadelphia, Pa., assignors 
to the U. S. Army. 


2,878,648 


Epitor’s Note: Patents listed above 
were selected from the Official Gazette of 
the U.S. Patent Office. Printed copies of 
patents may be obtained from the Com- 
missioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents. 
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Apparatus for mechanically moving a 
load i remote control, such as forcibly 
ejecting an aircraft canopy in an emer- 
gency escape. A cartridge actuated 
thruster is connected by tubing to a car- 
tridge actuated initiator located adjacent 
to the operator. 

Airfoil tes rature control means (2,- 
873,933). . A. Fanti, Springfield, Mass., 
assignor to United Aircraft Corp. 


LW. 


Metal airfoil adapted to be exposed to a 
high external heat source. Means located 
in the airfoil for varying the temperature 
over a minor portion of the chord to main- 
tain a more uniform temperature through- 
out the chord. 


Continuous combustion power plant 
(2,874,539). H. M. Fox, Bartlesville, 
os assignor to Phillips Petroleum 


ig oxidizer 


Method of operating a powerplant with 
a liquid oxidizer selected from the group 
consisting of white and red fuming nitric 
acids. The oxidizer and fuel are intro- 
duced into the combustion zone in streams 
at a pressure of 15 to 40 psi greater than 
the pressure in the combustion zone, so 
that gentle mixing results. 


Substituted furans as hypergolic fuels 
(2,874,535). A. L. Ayers and C. R. 
Scott, Bartlesville, Okla., assignors to 
Phillips Petroleum Co. 

Method of developing thrust by the 
combustion of bipropellant components 
in a reaction motor. Fuel comprises a 
hydrocarbon substituted furan. The total 
number of carbon atoms present in the 
formula is not greater than 16. 


Cooling means for tailpipe (2,874,536). J. 
A. Benson and G. W. Crosby, Saugus, 
Mass., assignors to General Electric Co. 

Thrust augmentation for a jet engine. 
Means for permitting axial expansion of 
the liner between the hot flame area and 
the afterburner. 


Thrust reverser for jet engine (2,874,538). 
R. G. Laucher and F. A. Coss, Simi, 
oo assignors to Marquardt Aircraft 


Variable area exit nozzle composed of 
segments pivotally connected to the end 
of the engine. Rearward movement of a 
shaft forms an opening between a support 
and a cover through which blocked ex- 
haust gases are expelled to obtain reverse 
thrust. 

Aircraft (2,874,920). 
St. Louis, Mo. 


G. E. Mallinckrodt, 


George F. McLaughlin, Contributor 


2,674,920 


First and second assemblies connected 
by bearing means adapting them for 
opposite rotation about an axis of spin. 
An engine with first and second reactive 
thrust components, the first attached to a 
a and the second comprising a 
yarrel-shaped wing surrounding the spin 
axis. 


Fuselage shaping to reduce shock waves 

at transonic and supersonic speeds 

——— R. T. Whitcomb, Hampton, 
a. 


A first cross-sectional area is taken in a 
plane perpendicular to the longitudinal axis. 
A second area is taken in planes tangent 
to the Mach cone whose axis is on the longi- 
tudinal axis, and projected to the air- 
stream, gradually increasing to a maximum 
and then decreasing from nose to tail. 


Variable gear (2,875,628). H. J. Walter, 
Huntsville, Ala. 

Two conically-shaped sets of gears with 
gradually increased diameters and a single 
interposed spur gear. One set is mounted 
on one shaft driven by the driving shaft; 
the other is mounted on a parallel shaft. 
The spur gear is shiftable stepwise. 


Jet propulsion power plant (2,875,576). 
J. Endres, Wackersberg, Germany. 


A chamber defined between an inner 
and an outer cylinder and open at one end 
to form a discharge nozzle. A pair of 
—— pistons operable in the inner 
cylinder form a_ central combustion 
chamber between the pistons and end 
combustion chambers. 


et spoiler (2,876,619). A. P. Johnstone, 
rimley, England, assignor to Power 
Jets (Research and Development) Ltd. 
Pair of semicircular members at the 
discharge end of a jet stream. Members 
may be moved through 90 deg with their 
concave sides facing upstream to a spoil- 
ing position. 
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OF A SERIES 


Oldsmobile Engineering Leadership sets the 
industry pace with a unique electronic wheel 
alignment device that dynamically computes 
**toe-in’” measurements for precision steering 
and handling. 


Handling and steering ease depend upon precise, 
minute measurement and control of front wheel align- 
ment. Because wheels have a tendency to “‘toe- 
out” when in motion, they must be adjusted for a 
slight amount of “‘toe-in” to eliminate ‘wheel fight’’, 
wander and undue tire wear. 


To meet the requirement of rapid, yet extremely 
accurate measurements on the production line, Olds- 
mobile engineers developed an electronic computor— 


OLDSMOBILE “TOES THE MARK’’ 
ELECTRONICALLY! 


a linear-differential-variable transformer—that dynami- 
cally and accurately measures the average amount of 
toe-in within .030 inches. As the car is brought into 
position, the wheels are rotated by rollers to simulate 
actual driving conditions and to eliminate errors 
caused by variations in tire run-out. By watching the 
visual gauges, an operator can quickly make the 
necessary adjustments to the steering linkage. 


By using the most up-to-date electronic measuring 
techniques in engineering and manufacturing, Olds- 
mobile is able to offer safe, accurate steering and 
handling . . . a controlled, comfortable ride. Visit your 
local Oldsmobile Quality Dealer, take a ride in a °59 
Oldsmobile and see why it’s the value leader of its class! 
OLDSMOBILE DIVISION e GENERAL MOTORS CORPORATION 


OLDS MOBILE > Where Proven Quality is Standard! 
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Book Reviews 


The Theory and Design of Magnetic 
Amplifiers, by E. H. Frost-Smith, John 
Wiley and Sons, Inc., New York, 1958, 
xix + 481 pp. 

Reviewed by G. J. Murpuy 
Northwestern University 


In spite of a statement to the contrary 
in the Foreword to this book, it does seem 
a peculiar choice for the first volume in a 
series on automation and control engineer- 
ing. Although the book contains much 
information of value to the designer of 
magnetic amplifiers, it appears to have 
been written neither by nor for a control 
systems engineer. 

An elementary introduction is followed 
by five chapters devoted to the transducer 
(i.e., the basic amplifying element in a 
magnetic amplifier). The topics dealt 
with in these chapters include the basic 
operation and the characteristics of trans- 
ductors, the response of a transductor with 
resistive load to a step input, the steady- 
state response to a sinusoidal input, the 
effects of self-excitation and of auto-exci- 
tation, and the effect of introducing reac- 
tive load components. 

Chapter 7, entitled ‘Closed Loop Con- 
trol Systems,” is simply a very brief (12 
pages) treatment of the relation between 
output and input of single-loop systems 
with loop transfer functions that have a 
constant numerator and a denominator of 
no more than second degree. This chap- 
ter is extremely weak and again raises the 
question of the purpose of including this 
book in a series of volumes on automation 
and control engineering. 

Chapter 8 is devoted to the transient 
response and the frequency response of 
transductors with self-excitation. 

The subject of the following three chap- 
ters is the use of transductors in various 
magnetic amplifier circuits. Magnetic 
modulators are discussed in chapter 12. 
In chapter 13, the author discusses the con- 
struction of magnetic amplifiers, suggests a 
design procedure, and illustrates the use of 
the latter in some practical examples. 

The fourteenth chapter is devoted to a 
brief, qualitative treatment of a few sys- 
tems in which magnetic amplifiers are used. 
Although a detailed study of some of these 
systems would be appropriate and of con- 
siderable interest, no such study is in- 
cluded. 

A short list of references is included, 
following the appendixes. The references 
seem to be well-chosen, but the most re- 
cent was published in 1955. 

The text is supplemented with approxi- 
mately 20 photographs of equipment, in 
addition to numerous line drawings. 

The impression of this reviewer is that 
this book may be of considerable interest 
to the designer of magnetic amplifiers and 
possibly to an applications engineer, but 
that it is likely to be very unsatisfying to 
the theoretically inclined control systems 
engineer. 


Ali Bulent Cambel, Northwestern University, Associate Editor 


Aerodynamics of Supersonic Flight, by 
Alan Pope, Pitman Publishing Corp., 
New York, 1958, ix + 243 pp. $6.00. 

Reviewed by P. Roy CuoupHurRY 
University of Southern California 


This is a second edition of “‘Aerodynam- 
ics of Supersonic Flight’’ which appeared 
first in 1950. These two editions are very 
much similar except for an additional chap- 
ter on “Performance Calculations” in the 
second edition. Some new materials are 
also added in the chapter on ‘Wind Tun- 
nel Design.” 

The author has tried “‘to provide a text 
at the lowest possible level for a one term 
course in supersonics,’’ and for that rea- 
son, perhaps, he has introduced a great 
number of topics in a short volume such as 
this. The topics include fundamental 
relations, flow in a duct, two-dimensional 
flow, the oblique shock, approximate 
theories, supersonic wind tunnels, three- 
dimensional flow and performance calcula- 
tions. Like the first edition, the coverage 
is rather sketchy. At a first glance, the 
book seems to be a collection of formulas, 
and many of the explanations seem rather 
shallow. This is because the author fol- 
lows a philosophy which educators are 
currently trying to do away with. He 
states this philosophy in the preface: 
“The course treatment is deliberately as 
non-thermodynamic as possible.’’ In most 
engineering schools emphasis is placed on 
the conservation principles (conservation 
of mass energy and momentum) which 
form the foundation of most engineering 
subjects. To try to teach a course in 
compressible fluids without thermodynam- 
ics is not wise and, needless to say, will 
result in confusion. This is perhaps the 
weakest point of the text. 


The Airglow and the Aurorae, edited by 
E. B. Armstrong and A. Dalgarno, 
Pergamon Press, New York, 420 pp. 
$22.50. 

Reviewed by S. F. Sincer 
University of Maryland 


This volume is the fifth one published 
as a supplement to the Journal of Atmos- 
pheric and Terrestrial Physics. Volume 1 
in this series will be remembered as the 
highly successful ‘Rocket Exploration 
of the Upper Atmosphere.’”’ The present 
volume is very much in the same style. 
It consists of 53 papers presented at the 
Conference of Airglow and Aurorae held 
at the Queen’s University of Belfast on 
September 6-7, 1955, and presents out- 
standing contributions from workers in 
all those fields of research which lead to a 
better understanding of the upper atmos- 
phere. 

Following an introductory paper by 
Prof. Sydney Chapman, the contributions 


are in seven groups. The first deals 
with the observations of the airglow, the 
luminous emission from the upper atmos- 
phere which is produced mainly by photo- 
chemical activity. The discussion centers 
around the emission from excited oxygen 
atoms, the OH bands and the sodium in the 
upper atmosphere, particularly the twi- 
light sodium emission. 

The next section deals with observations 
of aurorae and includes a discussion of 
radio investigations, in addition to the 
optical measurements. 

The next two sections deal with theo- 
retical considerations on how the airglow 
and aurorae may be excited. Some of the 
papers deal with the spectroscopic prob- 
lems of the light emission; others concern 
the ultimate cause of the aurora. There 
is, however, no good review of current 
auroral theories among these contribu- 
tions. 

The next three sections deal with 
laboratory studies in which luminosities 
similar to those in the atmosphere are 
excited, quantal calculations which relate 
the observed amounts of light to the num- 
ber of reactions going on, and finally an 
excellent section on instrumentation, an 
area in which, as can be imagined, the 
problems are severe. High speed in- 
struments collect as much of the faint 
light as possible, and the problem of high 
resolution is important for measuring the 
structure of spectral lines and line profiles. 

The volume is well-produced, especially 
the figures, and has subject and name 
indexes. 


A Guide To Nuclear Energy, by R. F. K. 
Belchem, Philosophical Library, New 
York, 1958, vii + 77 pp. 

Reviewed by J. P. FRANKEL 
University of California 


This is a beginners’ primer with sur- 
prisingly detailed information about re- 
actor principles, process and_ practice. 
Its contents range from basic principles 
of nuclear fission through types of reactors, 
fuels, fuel elements, materials, reprocessing 
(no details here) and applications of 
radiation. 

This is the kind of book that can be 
read by the average high school student in 
a single sitting; he would benefit greatly 
from the experience. The rest of us will 
have to be content with the knowledge 
that such a book is available to the scientifi- 
cally uninformed in readable form. 

The diction is British; it is a little 
startling, for instance, to read about 
fission products thrown off in the 
fission process. ..’’ when one remembers 
that the mean free path of fission frag- 
ments is of the order of millionths of a 
meter. 

The book is recommended for the young 
scientist-to-be. 
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Technical Literature Digest 


Jet Propulsion Engines 


Theoretical and Experimental Study 
of a Propulsion System, by R. Corelli, 
Fusées et Recherche Aéron., vol. 3, no. 2, 
July 1958, pp. 69-73. (In French.) 

Contribution to the Study of Normal 
Combustion in Engines Having Controlled 
Ignition, by B. Sale, Inst. Franc. Pétrole, 
Rév. et Ann. des Combustibles Liquides, vol. 
13, July-Aug. 1958, pp. 1157-1196. (In 
French.) 

A Method for Determining Turbine De- 
sign Characteristics for Rocket Turbodrive 
Applications, by Warner L. Stewart, David 
G. Evans and Warren J. Whitney, NACA 
Fes. Mem. E57K25a, May 1958, 32 pp. 
(Declassified from Confidential by author- 
ity of NASA Pub. Announcement 2, p. 17, 
12/10/58.) 

Rocket Research in Canada, by R. F. 
Wilkinson, Inst. Aeron. Sci., Preprint, 
853, Oct. 1958, 22 pp. 


Eprror’s Note: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged. 
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M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


The Rocket Engine Test Facility Sets 
the Pace, by Benjamin F. Rose Jr., 
Inst. Aeron. Sci., Preprint 840, July 1958, 

9 pp. 

Design of Ramjet Engine Ground Test 
Facilities, by S. b. Fasken, Inst. Aeron. 
Sci., Preprint 841, July 1958, 30 pp. 

Some Operating Experience and Prob- 
lems Encountered During Operation of 
a Free-jet Facility, by John E. McAulay 
and William R. Prince, NACA Res. Mem. 
E56J02, Feb. 1957, 21 pp. (Declassified 
from Confidential by authority of NASA 
Pub. Announcement 2, p. 26, 12/10/58.) 


Experimental Study of Shock-position- 
ing Method of Ram-jet-engine Control, 
by Herbert G. Hurrell, George Vasu and 
William R. Dunbar, NACA Res. Mem. 
E55F21, Aug. 1955, 71 pp. (Declassified 
from Confidential by authority of NASA 
Pub. Announcement 2, p. 24, 12/10/58.) 


Evaluation of Operating Characteristics 
of a Supersonic Free-jet Facility for Full- 
scale Ram-jet Investigations, by Carl B. 
Wentworth, Herbert G. Hurrell and 
Shigeo Nakanishi, NACA Res. Mem. 
E52108, Oct. 1952, 27 pp. (Declassified 
from Confidential by authority of NASA 
Pub. Announcement 2, p. 20, 12/10/58.) 


Some Control Considerations for Ram- 


jet Engines, by Seymour G. Himmel, 
NACA Res. Mem. E52F10, Aug. 1952, 
40 pp. (Declassified from Confidential 
by authority of NASA Pub. Announce- 
ment 2, p. 19, 12/10/58.) 

Design and Applications of Mach 4 
Turbine Engines, by M. A. Zipkin and 
R. E. Neitzel, Inst. Aeron. Sci., Preprint 
859, Oct. 1958, 14 pp. 

High-altitude Performance of an Ex- 
perimental Tubular Prevaporizing Com- 
bustor, by Helmut F. Butze, NACA Res. 
Mem. H54110, Dec. 1954, 29 pp. (De- 
classified from Confidential by authority 
of NASA Pub. Announcement 2, p. 22, 
12/10/58.) 

On the Measurement of High Fre- 
quency Stability Limits in a Liquid Bi- 
propellant Rocket, by L. Crocco and J. 
Grey, Aerotecnica, vol. 38, June 1958, pp. 
135-144. (In Italian.) 

An Examination of Solutions of the 
Underwater Ramjet for Aerial Tor- 
pedoes, by Dino Dini, Aerotecnica, vol. 38, 
June 1958, pp. 145-154. (In Italian.) 

Gas Turbines... Mobile, Flexible Mis- 
sile Support, by P. B. Garner, Space/- 
Aeron., vol. 30, Dec. 1958, pp. 38-41. 

Rocket Throttling: Tricky but Not Im- 
possible, by Kurt R. Stehling, Space/- 


The Applied Research & Development 
Division of Republic Aviation is making 
a multi-disciplined attack on the com- 
plex problems posed by advanced and 
unique flight vehicles. This expanding 
effort affords you the opportunity to 
apply your technical imagination and 


originality in the following areas: 


ASTRODYNAMICS 


e Space Vehicle Design 
e Orbital Mechanics 

e Interplanetary Flight 
e Re-Entry Flight 


AERODYNAMICS 


@ Missiles & Aircraft Aerodynamics 


e Hypersonic Flow 
e Physics of Rarefied Gases 
e Strong Shock Phenomena 


engineers « scientists 


Explore Previously 


Uncharted Areas in 


ASTRODYNAMICS and 


AERODYNAMICS 


There’s a growing atmosphere of technological excitement these days 

at Republic Aviation Corporation. Accelerating space study contracts 

plus a number of original proposals in such areas as ICBM’s, air-to-surface 
ballistic missiles and unique aircraft are responsible for this environment 
that adds incentive for you to use your best thinking. The scope 


of challenge at Republic ranges from design concept to proof of 
feasibility. Working within a flexible framework of theoretical and 
experimental research, you will be able to define your own approach 
‘n attaining the requirements of your assignment. 


For further information on these positions which combine 
scientific pioneering with career growth, address your resume to: 


Mr. George R. Hickman, Engineering Employment Manager, Dept. 10H 


FTE 


Farmingdale, Long Island, New York 
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Aeron., vol. 30, Dec. 1958, pp. 66-71. 

The World’s Aircraft Engines, nteravia, 
vol. 13, Dec. 1958, p. 1279. 

Soviet Development of Nuclear Air- 
craft Engines, by G. N. Nesterenko, A. I. 
Sobolev and Yu. N. Sushkov, Space/- 
Aeron., vol. 31, Jan. 1959, pp. 132-136. 


Aerodynamics of Jet 
Propelled Missiles 


On the Stability of Flow in the Bound- 
ary a Near the Nose of a Blunt Body, 
by A. M. Kuethe, Rand Corp., Res. Mem. 
RM-1972 (ASTIA AD 150687), Aug. 
1957, 13 pp. 

Thin Airfoils in Hypersonic Flow with 
Strong Shocks, by Joseph L. Raymond, 
Rand Corp., Paper P-1189, Dec. 1957, 14 
pp. 

Sweepback Theory for Shock Waves at 
Hypersonic Speeds, by J. D. Cole, Rand 
Corp., Paper P-1337, April 1958, 12 pp. 

Sweepback Theory for Shock Waves at 
Hypersonic Speeds, by J. D. Cole, Rand 
Corp., Res. Mem. RM 1991 (ASTIA AD 
144284), Oct. 1957, 13 pp. 

A New Approach to the Theory of 
Slender Body in Compressible Flow, by 
Haruo Moriguchi, J. Phys. Soc. Japan, 
vol. 13, Nov. 1958, pp. 1384-1389. 

Stability Derivatives for Spinning Bodies 
of Revolution in the Superaerodynamic 
Region, by Thomas W. Finch, Bur. Ordn., 
TN 43, 1958, 20 pp 

Experimental Studies at Mach Numbers 
12 to 19 of Conical and Blunted Bodies at 
Zero Angle of Attack, by 8S. M. Bogdonoff 
and I. E. Vas, Princetun Univ., Dept. 
Aeron. Engng., Rep. 435 (AFOSR-TN-58- 
pa ASTIA AD 203181), Sept. 1958, 

pp. 

Theory of Flight of the Sounding Rocket, 
by V. C. Liu, ASME, Paper 58-A-32, 
Nov. 1958, 3 pp. 

Axisymmetric Stagnant Flow of a Vis- 
cous and Electrically Conducting Fluid 
Near the Blunt Nose of a Spinning Body 
With Presence of Magnetic Field, II: 
Consideration of Realistic Conditions: 
Compressible Viscous Layer and Small 
Magnetic Reynolds Number, by Ching- 
Sheng Wu, Princeton Univ., Dept. Aeron. 
Engng., Rep. 431 (AFOSR-TN-58-711; 
ASTIA AD 162246), Sept. 1958, 36 pp. 

Laboratory Studies of Missile Re-entry 
Aero-thermodynamics, by D. Bershader, 
J. Instr. Soc. America, vol. 5, Nov. 1958, 
pp. 62-71. 

On the Effect of Flow Separation on the 
Lift of Slender Bodies, by Svetopolk 
Pivko, J. Roy. Aeron. Soc., vol. 62, Nov. 
1958, pp. 832-833. 

An Investigation of the Aerodynamic 
Characteristics of a Series of Cone- 
cylinder Configurations at a Mach Num- 
ber of 6.86, by Ralph D. Cooper and 
Raymond A. Robinson, NACA Res. Mem. 
L51J09, Dec. 1951, 43 pp. (Declassified 
from Confidential by authority of NASA 
Pub. Announcement 2, p. 27, 12/10/58.) 


Generation of High-velocity Projec- 
tiles, by Melvin A. Cook and Robert T. 
Keyes, J. Appl. Phys., vol. 29, Dec. 1958, 
pp. 1651-1657. 

Theoretical Investigation of the Flow 
Field about Blunt-nosed Bodies in Super- 
sonic Flight, by Roberto Vaglio-Laurin 
and Antonio Ferri, J. Aero/Space Sci., 
vol. 25, Dec. 1958, pp. 761-768. 


The Effects of Blunt Leading Edges on 
Delta Wings at March 5,8, by Kenneth F. 
Nicholson, J. Aero/Space Sci., vol. 25, 
Dec. 1958, pp. 768-787. 
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Heat Transfer and Fluid 
Flow 


Heat Transfer in the Turbulent In- 
compressible Boundary Layer, I: Con- 
stant Wall Temperature, by W. C. Rey- 
nolds, W. M. Kays and §. J. Kline, NASA 
Mem. 12-1-58W, Dec. 1958, 36 pp. 

Heat Transfer in the Turbulent In- 
compressible Boundary Layer, II: Step 
Wall-temperature Distribution, by W. C. 
Reynolds, W. M. Kays and §S. J. Kline, 
NASA Mem. 12-2-58W, Dec. 1958, 
38 pp. 

Heat Transfer in the Turbulent In- 
compressible Boundary Layer, III: Arbi- 
trary Wall Temperature and Heat Flux, 
by W. C. Reynolds, W. M. Kays and 
S. J. Kline, NASA Mem. 12-3-58W, Dec. 
1958, 52 pp. 

Heat Transfer in the Turbulent In- 
compressible Boundary Layer, IV: Effect 
of Location of Transition and Prediction 
of Heat Transfer in a Known Transition 
Region, by W. C. Reynolds, W. M. Kays 
and S. J. Kline, NASA Mem. 12-4-58W, 
Dec. 1958, 28 pp. 

Air-flow and Thrust Characteristics of 
Several Cylindrical Cooling-air Ejectors 
with a Primary to Secondary Temperature 
Ratio of 1.0, by W. K. Greathouse and 
D. P. Hollister, NACA Res. Mem. 
E52L24, March 1953, 79 pp. (Declassified 
from Confidential by authority of NASA 
Pub. Announcement 2, p. 21, 12/10/58.) 

Use of IBM 704 Computer for Com- 
plicated Thermodynamic Equilibria Cal- 
culations, by D. C. Villars, NAVORD 
Rep. 6382 (NOT'S 2064), July 1958, 12 pp. 

Triple Correlations, by M. E. Rose, 
AEC, ORNL-2516, May 1958, 60 pp. 

Real Gas Effects on Shock Tube Per- 
formance at High Shock Strengths, by 
J. L. Stollery, Gt. Brit., Aeron Res. Council, 
Current Paper 403, 1958, 12 pp. 12 figs. 

The Influence of the Absorption of 
Radiation in Shock Tube Phenomena, by 
Jacob Pomerantz, NAVORD Rep. 6136, 
Aug. 1958, 134 pp. 

A Study of Some Nearly Free Molecular 
Flow Problems, by D. Roger Willis, 
Princeton Univ., Dept. Aeron. Engng., 
Rep., 440, 1958, 115 pp., 33 figs. 

Explosively Induced Nonuniform Ob- 
lique Shocks, by John O. Erkman, Stanford 
Res. Inst., Poulter Labs., Tech. Rep.010-58, 
June 1958, 17 pp. 


The Development of a Variable Mach 
Number Effuser, by J. B. McGarry, Gt. 
Brit., Aeron. Res. Council, Rep. & 
Mem. 3097 (Formerly ARC Tech. Rep. 
19619; Natl. Gas Turbine Estab., 
Rep. R210), 1958, 15 pp. 

Static Crucible Investigation of Corro; 
sion and Mass-transfer of Nickel in Mol- 
ten Sodium Hydroxide with Various Addi- 
tives, by Americo F. Forestieri and 
William F. Zelezny, NACA Res. Mem. 
E53F12, Aug. 1953, 25 pp. 

Preliminary Investigation of Corrosion 
by Molten Sodium Hydroxide Flowing in 
Tubes of AISI 347 Stainless Steel, Inconel, 
and Nickel Having Average Outer-wall 
Temperatures of 1500° F and a Circum- 
ferential Temperature Gradient of 20° F, 
by Don R. Mosher and Leland G. Desmon, 
NACA Res. Mem. F51J16, March 1952, 
20 pp. 

Determination of High Gas Tempera- 
tures at High Speeds, II, by E. Neher, 
Archiv Fiir Tech. Messen, 273, Oct. 1958, 
pp. 205-208. (In German.) 

Solution to the Previous General Motors 
Institute Classroom Problem: Determine 
the Static and Total Temperatures of a 


igh Temperature, High Velocity Gas 
Seki by Kenneth W. Woodfield and 
Rodger Bloomfield, Gen. Motors Engng. J., 
vol. 5, Oct.-Dec. 1958, pp. 47-49. 

Unsteady Laminar Boundary Layers in 
an Incompressible Flow, by Kwang-Tzu 
Yang, J. Appl. Mech., vol. 25, Dec. 
1958, pp. 421-427. 

Thermoelectric Effects and Irreversible 
Thermodynamics, by G. N. Hatsopoulos 
and J. H. Keenan, J. Appl. Mech., vol. 
25, Dec. 1958, pp. 428-432. 

Comparison between the Linear and 
Nonlinear Steady-state Behavior of a 
Heated Tube, by B. L. McFarland, J. 
Appl. Phys., vol. 29, Dec. 1958, pp. 
1682-1684. 

Approximate Three-dimensional Solu- 
tions for Transient Temperature Distribu- 
bution in Shells of Revolution, by Maurice 
A. Brull and Jack R. Vinson, J. Aero/- 
Space Sci., vol. 25, Dec. 1958, pp. 742-750. 

Turbulent Heat Transfer Through a 
Highly Cooled, Partially Dissociated 
Boundary Layer, by Peter H. Rose, 
Ronald F. Probstein and Mac C. Adams, 
J. Aero/Space Sci., vol. 25, Dec. 1958, 
pp. 751-759. 

On Pitot Pressure in an Almost-free- 
molecule Flow—a Physical Theory for 
Rarefied-gas Flows, by V. C. Liu, J. 
Aero/Space Sci., vol. 25, Dec. 1958, pp. 
779-780. 

Experimental Mixing Profiles of a 
Mach 2.6 Free Jet, by Edward T. Pitkin 
and Irwin Glassman, J. Aero/Space Sci., 
vol. 25, Dec. 1958, pp. 791-792. 


On the Formation of Magneto-hydro- 
dynamic Shock Waves, by S. Segre, 
Nuovo Cimento, vol. 9, Sept. 16, 1958, 
pp. 1054-1057. 

Growth of Electric Space-charge and 
Radio Waves in Moving Ion Streams, by 
J. H. Piddington, Phil. Mag., Ser. 8 vol. 3, 
Nov. 1958, pp. 1241-1248. 

The Boundary Layer Equations of 
Compressible Flow. Separation, by D. 
Meksyn, Zeitschr. f. Angewandte Math. wu. 
Mech., vol. 38, Sept.-Oct. 1958, pp. 372- 
379. 

Experimental Investigation of Harmonic 
Oscillations in Rotating Flow, by H. Oser, 
Zeitschr. f. Angewandte Math. u. Mech., 
vol. 38, Sept.-Oct. 1958, pp. 386-391. 
(In German.) 

The Exterior Limiting Conditions of the 
Boundary Layer Differential Equations, 
by K. Nickel, Zeitschr. f. Angewandte 
Math. u. Mech., vol. 38, Sept.-Oct. 
1958, pp. 400-401. (In Ginn) 


Heat Transfer and Reynolds’ Analogy 
in a Turbulent Flow with Heat Release, 
by R. R. Hunziker, Zetischr. f. Angewandte 
Math. u. Phys., vol. 9a, no. 4, Nov. 11, 
1958, pp. 307-315. (In English.) 


Low Prandtl-number Free Convection, 
by J. L. Gregg and E. M. Sparrow, 
Zeitschr. f. Angewandte Math. u. Phys., 
vol. 9a, no. 4, Nov. 11, 1958, pp. 383-387. 
(In English.) 

A Simple Relationship Between the 
Shock and Expansion Pressure Coefficients 
as a Basis for Studying Two Dimensional 
Hypersonic Flow, by J. L. Raymond, 
Rand Corp., Paper P-1011, Jan. 1957, 
45 pp. 

Transport Coefficients of Dissociating 
and Slightly Ionizing Air, by P. S. Frei- 
finger, Rand. Corp., Res. Mem. RM-1794 
(ASTIA AD 133002), April 1957, 25 pp. 


Surface Protection and Cooling Systems 
for High Speed Flight, by D. J. Mason 
and C. Gazley Jr., Rand Corp. Res. Mem. 
RM-1735, March 1956, 42 pp. 


On the Propagation of Thermal Waves 
ARS JourNAL 


‘ 


in Media with Stored Energy, by A. J. 
ik. Foreman, Gt. Brit., Atomic Energy 
Res. Estab. AERE T/M 162, 1958, 10 
pp. 

Investigation of the Transmission of a 
Shock Wave Through an Orifice, by Louis 
L. Monroe, Calif. Inst. Techn., 
heim Aeron. Lab., Hypersonic Res. Proj., 
Mem. 46, Sept. 1958, 54 pp. 


Combustion, Fuels and 
Propellants 


A Modified Analysis for the Determina- 
tion of the Burning Velocity of a Gas 
Mixture in a Spherical Constant Volume 
Combustion Vessel, by K. H. O’Donovan 
snd C. J. Rallis, Witwatersrand Univ., 
Dept. Mech. Engng., Rep 1/58, May 1958, 
pp. 

Ignition of and Flames above Hydro- 
gen Peroxide Solutions, by Charles N. 
“atterfield, Ephraim Kehat and Maria 
\. T. Mendes, MIT, Dept. Chem. Engng., 
ivep. 54, Oct. 1958, 15 pp. 

Effect of Fuel Variables on Carbon For- 
mation in Turbojet-engine Combustors, 
iy Edmund R. Jonash, Jerrold D. Wear 
and William P. Cook, NACA Rep. 1352, 
(958, 18 pp. 

Spontaneous Flammability of Penta- 
borane and Pentaborane-3-methyl- 
pentane Blends, by Edward A. Fletcher, 
VACA Res. Mem. E53117, Feb. 1957, 
il pp. (Declassified from Confidential 
by authority of NASA Pub. Announcement 
i, p. 12,11/14/58.) 

Performance Estimation of the Coolant 
Fuel System, Final Summary Report, 
Odin Associates, Tech. Rep. 100-1 
(\AFOSR-T R-58-43; ASTIA A D 154150), 
Dee. 1957, 23 pp. 

Development of a Double Base Pro- 
pellant Charge for T17E2 JATO, by John 
Berti and Theodore Zimmerman, Picatinny 
Arsenal, Feltman Res. and Engng. Labs., 
Tech. Rep. 2540 (ASTIA AD 201101), 
Sept. 1958, 18 pp. 

A Pressure Bomb for the Photography 
of Burning Strands of Solid Propellant, 
by David R. Margetts and Benjamin V. 
Connor, Calif. Inst. Tech., Jet Prop. 
Lab., Rep. 20-119, July 1958, 21 pp. 

Energy Partition in the Exploding Wire 
Phenomenon, by F. D. Bennett, Aberdeen 
Proving Ground, Ball. Res. Labs., Rep. 
1056, Oct. 1958, 32 pp. 

Development and Evaluation of New 
Nitrometer Reaction Bulb, by Milton 
Roth, Raymond F. Wegman and George 
Albansoder, Picatinny Arsenal, Feltman 
Res. and Engng. Labs., Tech. Rep. 2579 
(ASTIA AD 203400), Nov. 1958, 11 pp. 

On the Stability of a Plane Deflagration 
Wave, by J. Mendes, Calif. Inst. Tech., 
Jet Prop. Lab., Progr. Rep. 20-360, 
Aug. 1958, 18 pp. 

A Contribution to the Theory of Lami- 
nar Flames with Radial Symmetry, by 
J. Menkes, Calif. Inst. Tech., Jet Prop. 
Lab., Progr. Rep. 20-361, Aug. 1958, 29 pp. 

An Experimental Investigation on Flame 
Stabilization in a Heated Turbulent 
Boundary Layer, by Donald Lawson Tur- 
cotte, Calif. Inst. Tech., Jet Prop. Lab., 
Progr. Rep. 20-355, June 1958, 88 pp. 

Experiments Concerning the Occurrence 
and Mechanism of High Frequency Com- 
bustion Instability, by Calvin Larue 
Barker, Calif. Inst. Tech., Jet Prop. 
Lab., Progr. Rep. 20-356, Aug. 1958, 92 pp. 

The Effect of Fuel Types and Admission 
Method upon Combustion Efficiency, by 
H. M. MeManus Jr., W. E. Ibele and 
T. Ek. Murphy, ASME, Paper 58-A-128, 
Nov. 1958, 4 pp. 


Aveust 1959 


Explore 
new areas in 
mathematics 


at IBM 


Creative IBM mathematicians are discovering 
important, new applications of mathematics 
in the electronic computer field. With the aid 
of large-scale computers, for example, one 
group of mathematicians simulated ina matter 
of weeks, eight years of engineering work which 
have yet to begin. Other IBM mathematicians 
are studying vehicular penetration problems 
involving thousands of variables. Rewarding 
careers are now available in such areas as 
computer system design and analysis; human 
factors engineering; mathematical and numeri- 
cal analysis; probability, communication, and 
information theory; reliability; scientific pro- 
gramming; statistics; and switching theory. 


A career with IBM offers advancement oppor- 
tunities and rewards. You will enjoy profes- 
sional freedom, comprehensive education 
programs, and the assistance of skilled special- 
ists. Working independently or as a member of 
a small team, your contributions are quickly 
recognized. This is a unique opportunity for a 
career with a company that has an outstand- 
ing growth record. 


Qualifications: B.S., M.S., or Ph.D. in Mathe- 
matics, Physics, Statistics, Engineering Sci- 
ence, or Electrical Engineering—and proven 
ability to assume important technical respon- 
sibilities in your sphere of interest. 

For details, write, outlining background and 
interests, to: 


Mr. R. E. Rodgers, Dept. 572H 
IBM Corporation 

590 Madison Avenue 

New York 22, N. Y. 


INTERNATIONAL BUSINESS 
MACHINES CORPORATION 
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Senior staff 
openings 

in basic 

research 


: THE BOEING SCIENTIFIC RESEARCH 
o LABORATORIES are engaged in a 
program of fundamental research 
designed to make contributions to 
the progress of the physical 
ae sciences. High-level staff positions 
abs 7 are open now in the fields of 


a Gas Dynamics 
Plasma Physics 
Mathematics 

Solid State Physics 
Electronics 
Physical Chemistry 
Geoastrophysics 


Boeing grants scientists the lati- 
tude and independence needed to 
achieve and maintain leadership in 
their special fields. Scientists inter- 
ested in carrying on their work in 
this kind of stimulating research 
environment are invited to com- 

municate with Mr. G. L. Hollings- 

: worth, Director, Boeing Scientific 
Research Laboratories. 


P. O. Box 3822 -JTB 
Boeing Airplane Company 
| Seattle 24, Washington 


BOEING 


Heat Transfer to Anhydrous Ammonia 
Flowing Turbulently in an Unsymmet- 
rically Heated Rectangular Tube, by 
R. Mascolo, North American Aviation, 
Inc., Rocketdyne Div., Rep. R-1155, Nov. 
1958, 46 pp. 

Cooperative Test Projects of the Joint 
Army Navy Air Force Panel on Analytical 
Chemistry of Solid Propellants, Round 
Robin 4: Ethyl Centralite and Nitro- 
glycerin in Cordite N; Round Robin 5: 
Heat of Explosion; Round Robin 6: 
Triacetin in OGK Propellant, by R. H. 
Pierson, Naval Ordn. Test Sta., NOTS 
1937, part 3, Feb. 1958, 60 pp. 

Halogen Fluorides and Other Covalent 
Fluorides, Physical and Chemical Prop- 
erties, by Howard C. Clark, Chem. Revs., 
vol. 58, Oct. 1958, pp. 869-894. 

A Theory for Estimation of Interfacial 
Energies, II: Application to Surface 
Thermodynamics of Tefion and Graphite, 
by Robert J. Good, Louis A. Girifalco 
and Gerard Kraus, J. Phys. Chem., vol. 
62, Nov. 1958, pp. 1418-1421. 

The High Temperature Heat Contents 
of Sodium Ferrite and A-Aluminum Oxide, 
by John L. Margrave and Robert T. 
Grimley, J. Phys. Chem., vol. 62, Nov. 
1958, pp. 1436-1438. 

Heats of Formation of Diborane and 
Pentaborane, by Edward J. Prosen, 
Walter H. Johnson and Florence Y. 
Pergiel, J. Res. Nat. Bur. Standards, vol. 
61, Oct. 1958, pp. 247-250. 

Infrared Emission Spectra of Flames 
under High Resolution, by Earle K. Plyler 
and Eugene D. Tidwell, J. Res. Nat. 
iol Standards, vol. 61, Oct. 1958, pp. 263- 
68. 


Boron Plus Deuteron Reactions, by 
R. W. Kavanagh and C. A. Barnes, 
Phys. Rev., vol. 112, Oct. 15, 1958, pp. 
503-511. 

Cyclotron Resonance in Flames, by 
Jiirgen Schneider and F. W. Hofmann, 
Phys. Rev. Letters, vol. 1, no. 11, Dee. 
1, 1958, pp. 408-409. 

Photometric Investigations of Alkali 
Metals in Hydrogen Flame Gases, III: 
The Source:of the Alkali Metal Continuum 
by C. G. James and T. M. Sugden, Proc., 
Roy. Soc., London, vol. A248, Nov. 11, 
1958, pp. 238-247. 

Photometric Investigations of Alkali 
Metals in Hydrogen Flame Gases, IV: 
Thermal and Chemiluminescent Effects 
Produced by Free Radicals, by P. J. 
Padley and T. M. Sugden, Proc., Roy. 
Soc., London, vol. A248, Nov. 11, 1958, 
pp. 248-265. 

Analytical Evaluation of Effect of Inlet- 
air Temperature and Combustion Per- 
formance of Boron Slurries and Blends 
of Pentaborane in Octene-1; Supple- 
ment I: Influence of New Boric-oxide 
Vapor-pressure Data on Calculated Per- 
formance of Pentaborane, by Leonard K. 
Tower, NACA Res. Mem. E56D02, 
May 1956, 11 pp. (Declassified from 
Confidential by authority of NASA 
Pub. Announcement 1, p. 16, 11/14/58.) 

Analytical Evaluation of Effect of 
Inlet-air Temperature and Combustion 
Pressure on Combustion Performance of 
Boron Slurries and Blends of Pentaborane 
in Octene-1, by Leonard K. Tower, 
paw Res. Mem. E55A31, June 1955, 

_ (Declassified from Confidential by 
aut ority of NASA Pub. Announcement 
1, p. 15, 11/14/58). 


Materials of Construction 


High-temperature Mechanical Proper- 
ties of Graphite, I: Creep in Compres- 


sion, by Paul Wagner and Allen R. 
Driesner, J. Appl. Phys., vol. 30, Feb. 
1959, pp. 148-151. 

High-temperature Properties of Graph- 
ite, II: Creep in Tension, by Pau! 
Wagner, Allen R. Driesner and Larry A. 
Haskin, J. Appl. Phys., vol. 30, Feb. 
1959, pp. 152-154. 

New Aspects in Ceramic Coatings, by 
Paul A. Huppert, ARS Journat, vol. 29, 
Jan. 1959, pp. 19-21. 

Measurements of Total Hemispherica! 
Emissivity of Several Stably Oxidized 
Metals and Some Refractory Oxide 
Coatings, by William R. Wade, NASA 
Mem. 1-20-59L, Jan. 1959, 3 pp. 

Approximate Equations for Intercon- 
verting the Various Mechanical Propertie: 
of Linear Viscoelastic Materials, by) 
Thor L. Smith, Calif. Inst. Tech., Je’ 
Prop. Lab., Mem. 20-176, Nov. 1958, 
27 pp. 

State of the Art: 
terials, Space-Aeron., vol. 31, 
1959, pp. 40-41. 


Structures and Ma- 
March 


Instrumentation, 
Telemetering, Data 
Recording 


Determination of the Speed of Moving 
Bodies with the Help of the Doppler 
Effect, by Siegmund Wintergerst, Arch. 
f. Tech. Messen, No. 275, Dec. 1958, pp. 
241-242. (In German.) 

New Strain Gages for the Space Age, 
ISA J., vol. 6, Feb. 1959, pp. 50-51. 

Device for Thermal Control of Liquid 
Hydrogen Bubble Chambers, by Pierre 
Amiot, Nuclear Instruments, vol. 3, no. 
5, Nov. 1958, pp. 275-277. 

Electronic Apparatus for the Measure- 
ment of Dynamic Magnetic Fields, by 
I. F. Quercia, Nuclear Instruments, vol. 
3, no. 5, Nov. 1958, pp. 292-297. 

Forced Precession of Symmetrical Gyro- 
scopes under Vertical Shock, by F. 
Weidenhammer, Zeitschr. f. Angewandte 
Math. u. Mech., vol. 38, Nov.-Dec. 
1958, pp. 480-483. (In German.) 

Electronic Flow Control, by J. G. 
Ziegler and N. B. Nichols, JSA J., vol. 
6, Jan. 1959, pp. 64, 69 

Avionics, Aviation Week, vol. 69, Mid- 
December 1958, Section C, 106 pp. 

Instruments: Key to Missile Program, 
by F. Leary, Electronics, vol. 32, Jan. 16, 
1959, pp. 47-51. 

Preliminary Development of Telemetry 
for Aeroballistic Ranges, by M. K. 
Kingery, Arnold Engng. Dev. Center, 
TN-59-7 (ASTIA 208692), Feb. 
1959, 25 pp. 


Electronic Transient Recorder for 
Tracking of High Speed Projectiles, 
by D. Ritke and E. White, Univ. of 
Maryland, Inst. Molecular Phys., IMP- 
ONR-1, Nov. 1958, 17 pp. 

Transducer Frequency Response Evalu- 
ation for Rocket Instability Research, 
by Charles R. Tallman, ARS Journat, 
vol. 29, Feb. 1959, pp. 119-122. 

A New Mass Flow Meter for Liquid 
Streams, by R. F. Wall, Ind. Engng. 
Chem, vol. 51, Feb. 1959, pp. 69A-70A. 

Some Effects of Vibration and Rotation 
on the Drift of Gyroscopic Instruments, 
by R. M. Stewart, ARS JourNaL, vol. 29, 
Jan. 1959, pp. 22- D8. 

A Radioactive Ionization Gage Pressure 
Measurement System, by N. W. Spencer 
and R. L. Boggess, ARS JouRNAL, vol. 29, 
Jan. 1959, pp. 68-71. 
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Guidance Systems and 
Components 


Selection of an Aerodynamic Configura- 
tion for Improved Beam Rider Guidance, 
by Basil Staros, Richard Gretz and Merven 
Mandel, ARS JourNAL, vol. 29, Feb. 1959, 
pp. 115-118. 

Guidance for Sputnick III; the Russian 
Satellite Made Use of Some Very Ingeni- 
ous Devices, by Richard E. Stockwell, 
Ordnance, vol. 48, Jan.-Feb. 1959, pp. 
675-675. 


Servomechanisms and 
Controls 


An Attitude Control System for Space 
Vehicles, by Walter Haeussermann, 
JouRNAL, vol. 29, March 1959, pp. 203-207. 

Spatial Attitude Control of a Spinning 
Rocket Cluster, by Walter Haeussermann, 
Al’S JouRNAL, vol. 29, Jan. 1959, pp. 56-58. 

The Cluster Spin Control System for 

upiter C Missiles, by Niels Jasper and 
ieter Teuber, ARS JouRNAL, vol. 29, Jan. 
1959, pp. 58-61. 


Flight Vehicle Design and 
Testing 


Missile Ground Equipment, by B. J. 
Meldrum, SAE Preprint 5T, Jan. 1959, 
8 pp. 

Precise Electric Power for Ground 
Support of Missiles, by David Ginsberg, 
SAE Preprint 7R, Jan. 1959, 5 pp. 

Design and Fabrication of a Re-entry 
Vehicle, by Leon L. Farnham, AstrRo- 
Nautics, vol. 4, March 1959, pp. 26-27, 
57-58. 

Heat Protection for Re-Entry Vehicles, 
by Norris F. Dow, Astronautics, vol. 
4, March 1959, pp. 28-29, 82, 84, 86. 

Missile Reliability Now; Surveillance 
and Research Programs Must Assure 
Weapon Dependendability, by T. H. 
Braunstein, Ordnance, vol. 43, Jan.-Feb. 
1959, pp. 686-687. 

Optimization of Rockets for Maximum 
Payload Energy, by Dandridge, M. Cole 
and Michael A. Marrese, ARS JouRNAL, 
vol. 29, Jan. 1959, pp. 71-73. 

Interrelationship of Calculus of Varia- 
tions and Ordinary Theory of Maxima 
and Minima for Flight Mechanics Ap- 
plications, by Angelo Miele, ARS Jour- 
NAL, vol. 29, Jan. 1959, pp. 75-76. 


Space Flight 


Centimetric Radar Observations of the 
Russian Satellites, by J. S. Hey and V. A. 
Hughes, Proc. Roy. Soc., vol. 248A, Oct. 
28, 1958, pp. 34-42. 

Method for Computations of Satellite 
Orbits, by L. Jacchia, Proc. Roy. Soc., 
vol. 248A, Oct. 28, 1958, pp. 43-44. 

The Prediction Service of H.M. Nauti- 
cal Almanac Office, by D. H. Sadler, Proc. 
Roy. Soc., vol. 248A, Oct. 28, 1958, pp. 
45-47. 

The Computation of Orbit Parameters 
from Interferometer and Doppler Data, 
by E. G. C. Burt, Proc. Roy. Soc., vol. 
248A, Oct. 28, 1958, pp. 48-54. 

Perturbations of the Orbit of a Satellite 
Near to the Earth, by D. G. King-Hele, 
Proc. Roy. Soc., vol. 248A, Oct. 28, 1958, 
pp. 55-62. 


Aueust 1959 


This difference between systems can make 


a big difference in your career 


IF YOU ARE QUALIFIED and interested in contributing to programs of 
“total” scope, it will be of value to you to investigate current oppor- 
tunities with General Electric’s DEFENSE SYSTEMS DEPT., whose 
work lies primarily in providing total solutions to large scale defense 
problems of the next 5, 10 and 20 years. 


The work here lies almost entirely in the areas of systems engi- 
neering and systems management. 


Inquire about these positions: 


Guidance Equation Engineers Systems Test Evaluation Engineers 

Systems Logistics Engineers Engineering Psychologists 

Electronic Systems Radar Equipment Engineers 
Management Engineers Weapons Analysis Engineers 

Operations Analysis Engineers Weapons Systems Integration 

Systems Program Engineers Engineers 

Data Processing Engineers Engineering Writers 


Forward your confidential resume at an early date. 

Whereas the growth potential is evident — both for DSD and the 

H9d9)))))) engineers who join us — the positions we fill during these early 
months will carry significant “ground-floor” benefits. 


; ‘y Write fully to Mr. E. A. Smith, Room 8-M. 


DEFENSE SYSTEMS DEPARTMENT 


A Department of the Defense Electronics Division 


GENERAL @ ELECTRIC 


300 South Geddes Street, Syracuse, N. Y. 
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ENGINEERS 
and 
SCIENTISTS 


Aeronutronic, a new divi- 
sion of Ford Motor Com- 
pany, has immediate need 
or qualified people to 
staff senior positions at 
its new $22 million Re- 
search Center in Ne rt 
Beach, Southern ali- 
fornia. 


The Space Technol 
Operation offers the 
highly desirable combina- 
tion of new facilities and 
advanced equipment, lo- 
cated in California’s finest 
environment for living 
and raising a family. In- 
vestigate these excep- 
tionally rewarding posi- 
tions now: 
VEHICLE TECHNOLOGY 
Aerodynamic design 
and testing 
Rocket engine 
development 
Rocket nozzle and re- 
entry materials 
High temperature 
chemical kinetics 
Combustion and 
detonation theory 
Combustion 
thermodynamics 
High temperature 
structural plastics 
& ceramics 
Advanced structures 
Rocket vehicle 
systems 
MISSILE DEFENSE 
Supersonic 
aerodynamics 
Aerothermodynamics 
High temperature heat 
transfer 
Space physics 
Re-entry programs 
ASTRO SCIENCES 
Space electronics 
Guidance & control 
Communications 
Instrumentation 
Experimental physics 
magneto- 
y ynamics 
studies 
Visit Aeronutronic’s ex- 
hibit booth 3822-24 at 
the WESCON show. 
Qualified applicants are invited 
to send resumes and inquiries to 
Mr. R. W. Speich, Aeronutronic, 


Dept. 19, jox 451, Newport 
Beach, California. 


AERONUTRONIC 


a division of 
FORD MOTOR COMPANY 
Newport Beach 


/ Santa Ana e Maywood, California 


Direction Finding Observations on the 
20Mc/s Transmissions from the Artificial 
Earth Satellites, by F. A. Kitchen, Proc. 
6 aaa vol. 248A, Oct. 28, 1958, pp. 


Radio Observations on 20 and 40 Mc/s, 
by W. C. Bain, Proc. Roy. Soc., vol. 248A, 
Oct. 28, 1958, pp. 68-71. 

Post Office Observations on the First 
Russian Satellite, by C. F. Booth, Proc. 
Roy. Soc., vol. 248A, Oct. 28, 1958, pp. 
72-75. 

The Ionosphere and the Radio Emission 
from the Satellites, by K. Weekes, Proc. 
Roy. Soc., vol. 248A, Oct. 28, 1958, pp. 
77-79. 

The Rocket’s Dimensions, by R. d’E. 
Atkinson, Proc. Roy. Soc., vol. 248A, 
Oct. 28, 1958, pp. 80-81. 

The British Astronomical Association’s 
and Radio Society of Great Britain’s 
Earth Satellite Observations, by J. Hey- 
wood, Proc. Roy. Soc., vol. 248A, Oct. 
28, 1958, pp. 82-87. 

Electric Power for Space Flight, by 
J. H. Huth, Rand. Corp., Paper P-1244, 
Dec. 1957, 16 pp. 

Power in Space, by J. H. Huth, Rand 
Corp., Paper P-1389, May 1958, 13 pp. 

Fuel Cells—Current Problems and 
Potential Usefulness, by J. H. Huth, 
Rand Corp., Paper P-1187, Sept. 1957, 
22 pp. 

Skin Temperature Variation During Re- 
Entry of Scientific Satellite, by D. J. 
Masson, Rand Corp., Res. Mem. RM- 
1693, March 1956, 8 pp. 

Heat Transfer Aspects of the Atmos- 
- Re-entry of Long Range Ballistic 

issiles, by Carl Gazley Jr., Rand Corp., 
Rep. R-273, Aug. 1954, 105 pp. 

The Penetration of Planetary Atmos- 
ae, by C. Gazley Jr., Rand Corp., 

‘aper P-1322, Feb. 1958, 47 pp. 
Lunar Trajectory Studies, by H. A. 


Lieske, Rand Corp., Paper P-1293, Feb. 
1958, 36 pp. 
Lunar Flight Trajectories, by R. W. 


Buckheim, orp., Paper P-1268, 
Jan. 1958, 49 pp. 

A Discussion of Energy Sources for 
Space Communications, by J. H. Huth, 
Rand Corp., Rep. P-1318, March 1958, 
25 pp. 

Propagation Considerations in Space 
Operations, by Cullen M. Crain, Rand 
Corp., Paper P-1305, March 1958, 17 
pp. 

Zero Gravity Effects Largely Unknown, 
by Brig. Gen. Don Flickinger, Aviation 
Week, vol. 70, Jan. 5, 1959, p. 35. 


Radio Tracks Orbits of Man-made 
Moons, by C. A. Schroeder, C. H. Looney 
Jr. and i. E. Carpenter Jr., Electronics, 
vol. 32, Jan. 2, 1959, pp. 33-37. 

Reliability and Stability Most Critical in 
Satellite Communication, by DuMont 
Labs, Space/Aeron., vol. 31, Jan. 1959, 
pp. 38-45. 

Atmospheric Attenuation of Space Com- 
munication, Space/Aeron., vol. 31, Jan. 
1959, pp. 74-75. 

Retro-rockets and Drag Devices 
Teamed for Satellite Recovery, by F. W. 
Ross, Space/Aeron., vol. 31, Jan. 1959, 
pp. 180-184. 

Three-dimensional Orbits of Earth 
Satellites, Including Effects of Earth 
Oblateness and Atmospheric Rotation, 
by Jack N. Nielsen, Frederick K. Good- 
win and William A. Mersman, NASA 
Memo 12-4-58A, Dec. 1958, 85 pp. 


Practical Limitations on Orbit Deter- 
mination, by Robert M. L. Baker Jr., 


Inst. Aeron. Sci., Preprint 842, July 1958, 
10 pp. 

Some Problems of Space Navigation, by 
Robert E. Roberson, Aero/Space Engng., 
vol. 17, Dec. 1958, pp. 39-43. 

SFI—A Prime Factor in Space Tech- 
nology, by Irwin Kane, Arma Engngq., vol. 
2, Oct.-Nov. 1958, pp. 25-57. 

Satellite Survey: Analysis of Their 
Orbits and Guidance, by Murray State- 
man, Arma Engng., vol. 2, Oct.-Noy. 
1958, pp. 28-31. 


Aerophysics, Astrophysics 


The ARDC Model Atmosphere, 1956, 
by R. A. Minzer and W. S. Ripley, «ir 
Force Cambridge Res. Center, Geophysics 
Res. Directorate, Air Force Surveys in 
Geophys. 86 (AFCRC-T N-56-204; AST/ A 
AD 110233), Dec. 1956, 201 pp. 

Measurements, by Two Methods, of 
the Equivalent Black Body Radiation 
Temperature in the Atmosphere and Infra- 
red Divergence Measurements in the 
Atmosphere with the Double Radiometer, 
by Roger L. Aagard, Minnesota Uni”, 
School of Physics, Dec. 1958, 35 pp. 

Variations in the Cosmic-ray Nucleonic 
Intensity, by J. A. Lockwood, Phi. 
Rev., vol. 112, Dec. 1, 1958, pp. 1750-1758. 

Suprathermal Particle Production by 
Binary Stars, by D. A. Tidman, Phys. 
Rev., vol. 112, Dec. 1, 1958, pp. 1759-1762. 

Acceleration Cosmic Radiation, ly 
G. S. Murty and R. K. Varma, Phys. 
Rev., vol. 112, Dec. 1, 1958, pp. 1789-1792. 

Unusual Cosmic-ray Intensity Fluctua- 
tions Observed at Southern Stations dur- 
ing October 21-24, 1957, by K. G. 
McCracken and N. R. Parsons, Phys. Rev., 
vol. 112, Dec. 1, 1958, pp. 1798-1801. 

Rapid Decrease of Cosmic Ray In- 
tensity, by J. F. Steljes and H. Carmichael, 
Nuovo Cimento, vol. 10, no. 2, Oct. 16, 
1958, pp. 393-394. 

International Conference on Cosmic 
Rays, Varenna, June 1957. Proceedings, 
Nuovo Cimento, Supplement, 2d Quarter, 
no. 2, 1958, pp. 126-804. 

Comparison of High-altitude Rocket and 
Satellite Density Measurements, by H. 
E. LaGow and R. Horowitz, Phys. of 
Fluids, vol. 1, no. 6, Nov.-Dec. 1958, 
pp. 478-479. 

Plasma Motions Induced by Satellites 
in the Ionosphere, by Lester Kraus and 
Kenneth M. Watson, Phys. of Fluids, 
vol. 1, no. 6, Nov.-Dec. 1958, pp. 480- 
488. 

Free-flight Test of a Technique for 
Inflating an NASA 12-foot-diameter 
Sphere at High Altitudes, by Alan B 
Kehlet and Herbert G. Patterson, NASA 
Mem. 2-5-59L, Jan. 1959, 18 pp. 

An Astronomical Telescope in Space, 
by R. J. Davis, F. L. Whipple and C. A. 
Whitney, Astronautical Sciences Rev., 
vol. 1, no. 1, Jan—March 1959, pp. 9-12. 

IGY Solar Flare Program and Ionizing 
Radiation in the Night Sky, by Herbert 
Friedman, ARS Journat, vol. 29, Feb. 
1959, pp. 103-107. 

U. S. Extension to ICAO Standard 
Atmosphere Tables and Data to 300 
Standard Geopotential Kilometers, by 
R. A. Minzner, W. S. Ripley and T. P. 
Condron, Air Force Cambridge Res. 
Center, Geophys. Directorate, 1958, 236 pp. 


Atomic Energy 


Potentialities and Problems of Nuclear 
Rocket Propulsion, by T. P. Cotter, Jnst. 
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Aeron. Sci., Rep. 59-24, Jan. 1959, 14 
pp 

Performance of Nuclear Electric Propul- 
sion Systems, by V. P. Kovacik and D. P. 
Ross, Inst. Aeron. Sci., Rep. 59-25, Jan. 
1959, 12 pp. 

A Method for Determining Core Dimen- 
sions of Heat Exchanger with One Dom- 
inating Film Resistance and Verification 
with Experimental Data, by John N. B. 
Livingood and Anthony J. Diaguila, 
NACA RM E56K26a, Jan. 1957, 16 pp. 

Potentialities and Problems of Nuclear 
Rocket Propulsion by T. P. Cotter, 
Aero/Space Engng., vol. 18, Feb. 1959, 
pp. 50-53. 

Kiwa-A to Test Nuclear Rocket Feas- 
ibility, by J. S. Butz Jr., Aviation Week, 
vol. 70, Feb. 16, 1959, pp. 48, 49, 51, 53. 

Piasma Motors, by Winston H. Bostick, 
J. Astron. Sci., vol. 5, Summer 1958, 
pp. 39-44. 

Performance Calculations for Hybrid 
Nuclear-Chemical Rocket Propulsion Sys- 
tems, by Leon Green Jr. and James M. 
Carier, ARS JourNaL, vol. 29, March 
1959, pp. 180-186. 

NASA Reactor Facility Hazards Sum- 
mary, Vol. II, by Lewis Res. Center, 
NASA Mem. 12-10-58E, Feb. 1959, 194 
pp. 

Lockheed’s Georgia Nuclear Labora- 
tory, by W. R. Rhoads, SAE Preprint 
28, Jan. 1959, 7 pp. 

A List of United Kingdom Atomic 
Energy Research Establishment Unclassi- 
fied Memoranda and Reports on Con- 
trolled Thermonuclear Reactions, by 
C. S. Sabel, Gt. Brit., Atomic Energy 
Res. Estab., AERE Inf/Bib 117, 1958, 
6 pp. 


Ignition of a Thermonuclear Plasma 
by High-energy Injection, by Albert 
Simon, Phys. of Fluids, vol. 1, no. 6, 
Nov.-Dec. 1958, pp. 489-494. 


Russian Technical Articles* 


Streamlined Bodies in Non-ideal Gas 
Flows at High Supersonic Speeds, by 
G. A. Liubimov, Akademii Nauk SSSR, 
Tzvestiia, Oldelenie Tekhnicheskikh Nauk, 
Mashinostroenie 1 Mekhanik, no. 1, 1959, 
pp. 173-176. (In Russian.) 

Calculation of Slip during Flow of a 
Viscous Gas around a Semi-infinite Plate, 
by V. P. Shidlovskii, Akademii Nauk 
SSSR, TIzvestiia, Otdelenie Tekhnicheskikh 
Nauk, no. 9, Sept. 1958, pp. 82-90. 
(In Russian. ) 

Radiation Cooling of Air, I: General 
Description of the Phenomenon and the 
Weak Cooling Wave, by Ia. B. Zel’dovich, 
A.S. Kompaneets and Iu. P. Raizer, Soviet 
Physics-JETP, vol. 34(7), no. 5, Nov. 
1958, pp. 882-892 

Stabilization of Temperature of Work- 
ing Substance of Thermal Resistance, 
by N. P. Udalov, Abtomatika i Telemekh- 


*The Editors of Technical Literature 
Digest are making a systematic search for 

rtinent Russian articles, both in the orig- 
inal Russian and in translated form. For 
the balance of 1959, these will appear in a 
separate section with this heading, mainly 
to draw the readers’ attention to the new 
listing. After that, they will appear item 
by item under the proper subject headings. 


anika, vol. 19, no. 11, 1958, pp. 1070- 
1072. (In Russian.) 

Some Problems of the Thermo-elasticity 
in Convection with Heat Stresses in 
Turbine Motors, by A. D. Kovalenko, 
Akademii Nauk SSSR, Izvestiia, Otdelenie 
Tekhnicheskikh Nauk, Oct. 1958, pp. 68- 
76. (In Russian.) 

Application of the Method of Electric 
Analogy to the Solution of the Problem of 
Radiant Heat Transfer, by V.N. Adrianov, 
Akademii Nauk SSSR, Izvestiia, Otdelenie 
Tekhnicheskikh Nauk, Energetika i Abto- 
matika, no. 1, 1959, pp. 20-25. (In 
Russian. ) 

On the Derivations of Basic Equations 
in the Hydrodynamic Theory of Heat 
Transfer, by P. N. Pinsker, Akademii 
Nauk SSSR, Izvestiia, Otdelenie Tekhni- 
cheskikh Nauk, Energetika 1 Abtomatika, 
no. 1, 1959, pp. 26-32. (In Russian.) 


Oblique Shock Waves in a Plasma with 
Finite Conductivity, by M. I. Kiselev 
and V. I. Tsepliaev, Soviet Physics- 
JETP, vol. 34(7), no. 6, Dec. 1958, pp. 
1104-1106. 

Theory of Excitation of Hydromagnetic 
Waves, by A. I. Akhiezer and A. G. 
Sitenko, Soviet Physics-JETP, vol. 35(8), 
no. 1, Jan. 1959, pp. 82-86. 

Cooling of Air by Radiation II: Strong 
Cooling Wave, by Ia. B. Zel’dovich, A. 8. 
Kompaneets and Iu. P. Raizer, Soviet 
Physics-JETP, vol. 34(7), no. 6, Dee. 
1958, pp. 1001-1006. 

Streamlined Conical Flow in Gases 
Moving at High Supersonic Speeds, by 
A. L. Gonor, Akademii Nauk SSSR, 
Tzvestiia, Otdelenie Tekhnicheskikh Nauk, 
Mashinostroenie i Mekhanik, no. 1, 1959, 
pp. 3-40. (In Russian. ) 
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